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Nonstationary Envelope Process and First 
Excursion Probability 

J.-N. Yang 

Engineering Mechanics Division 


A definition of the envelope of nonstationary random processes is proposed. The 
establishment of the envelope definition makes it possible to simulate the nonstation- 
ary random envelope directly. The envelope statistics, such as the derr'ty function, 
the joint density function, the moment function, and the level crossing rate, relevant 
to the analyses of the catastrophic failure, fatigue, and crack propagation of 
structures, are derived. Applications of the envelope statistics to the prediction of 
structural reliability under random loadings are demonstrated in detail. 


Introduction 

It is well k.iown that loadings to structures, such as earthquakes, 
oceanographical wave forces, strong gust loadings, aircraft impact loadings 
during landing, spacecraft excitations resulting from booster engine ignition 
or shutdown, etc., are nonstationary random processes. Because of the 
relative severity of such random excitations, the prediction of structural 
safety or failure becomes an important task of structural analysis, design, and 
test. 

In stationary random vibration, the envelope process plays an important 
role in the prediction of the catastrophic failure (e.g., References 1 and 2) 
fatigue failure, and crack propagation (e.g., Reference 3). Similarly, the 
nonstationary envelope process is expected to have important application to 
the prediction of structural failure under nonstationary random excitations. 

The purpose of this article is (1) to establish a definition of the envelope of 
nonstationary random processes so that the nonstationary random envelope 
processes can be simulated efficiently with the aid of a digital computer, and 
(2) to derive relevant envelope statistics useful in the prediction of structural 
reliability. 

The definition of stationary random envelope proposed by Rice (Refer- 
ence 4) and Cramer and Leadbetter (Reference 5) is extended herein to the 
envelope of nonstationary random processes possessing evolutionary power 
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spectral densities (e.g., References 6 and 7). The establishment of the 
nonstationary random envelope makes it possible to simulate the envelope, 
peaks, and troughs of a nonstationary random process directly (Reference 8). 
Statistics of envelope processes, such as the density function, the joint 
density function, the moment function, and the upcrossing rate of a barrier 
level, relevant to the analyses of the catastrophic failure probability, fatigue, 
and crack propagation, are derived. On the basis of these envelope statistics, 
approximate solutions to the first excursion probability are obtained. 

Sample functions of the structural response and envelope to strong gust 
loadings or earthquake excitation are simulated, and applications to the 
prediction of structural reliability are demonstrated in detail. 


Nonstationary Envelope Processes 


Consider a real nonstationary gaussian random process X(t) with mean 
zero and an evolutionary power spectral density. Then, X(t) has a spectral 
representation: 


X(t) = /_“ A(t, u)e lu 'dX(u) (1) 

where A(f,co) is a deterministic function of time t and frequency a>, called the 
oscillatory function. dX(o)) is an orthogonal random process with mean zero 
and 


E [dXiuJdX *(« 2 )] =0, w, *gj 2 
E [ldX(w)l 2 ]= S(u)du 


in which the asterisk denotes the complex conjugate; S(a>) is, in general, the 
power spectral density of a suitable stationary random process X(t), 


m = 


iu)t — 

re- dX(u) 


( 3 ) 


where e^dXia)) is the elementary oscillation of X(t). It is further assumed 
that X(t) *0 for t< 0 so that A(f,a>)—0 for K0. 
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Since X(t) is a real process, the real value autocorrelation function R x (ti,t 2 ) 
of X(£j) and X(t 2 ) can be shown as 


*7 



. t 




R x^r t 2^ = ^o’ lA ( t r w )! W(^2’^) |cos 

(4) 

+ 6(t J , oo) - 8(t 2 , co)] G(cj)dco 

where G(o>) is the one-sided spectral density of X(t) (Equation 3), i.e., 
G(co) = 2S((o ) for a>>0, and 

8(t, co) - tan -1 [\mA(t, co)/ReA(t, co)] (5) 

in which ImA(f,a>) and ReA(f,a>) are the imaginary part and the real part of 
A(t, (o), respectively. 

Furthermore, since X(t) is real, Equation 1 can be written as (Reference 5) 
X(t) = ClA(t, w)| {cos [cor + S(t, «)] dU(w) 

( 6 ) 

+ sin [cor + 5(/, co)] c/F(co)J 


where dU(a>) and dV(a)) are real orthogonal processes with mean zero, and 


E [dU(u) 2 ] =E [dV(u>f] = G(co)dco 


2 1 - 


(7) 


and dU( o>) and dV( a>) are mutually orthogonal. 

Following Cramer and Leadbetter, the stationary random process X(t) 
(Equation 3) is passed through a filter with gain g(o>), 


l, to < 0 

*(co) = {0, co = 0 

to > 0 


( 8 ) 


to obtain a stationary process X(t) with elementary .oscillation g(<o)e^dX(ct>). 
Then, one forms the nonstationary random process A (t) as follows: 


= 0>)g(u)e iu, dX(u) 


(9) 
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Since X(t) is real, it follows from Equation 9 that X(t) is real, and it can be 
written as 


)£{t) - J r 0 °° \A(t, co)l{sin [w t + S(t, w)] dU( co) 

-cos [atf + 60, co)] dV{co) ^ 

The envelope V(t) of X(t) is defined as 

W «[**(/) + **(f)] !/l (ii) 

where it is clear that 

V(t)> 1*0)1 

V(t) = 1*0)1. when X(t) = 0 

This definition of envelope (Equation 12) is originally due to Rice (Reference 

4). 

Equation 6 states that the nonstationary random process X(t) is a 
superposition (in the frequency domain) at each time t of elementary 
components 

l40» w)l {cos [cot + 50, w)] dU(co) 

+ sin [wf + 60, cc)] dV(co)\ 

= F(t, co) cos [cot + 60, co) + \p(co)] 

while Equation 10 states that X(t ) is a superposition of elementary 
components 

L40, w)l{sin [cot +50, to)] dU(co) 

- cos M + 8(/,o,)] 

= r{t, co) sin [cot + 60, co) + 4>(co)\ 
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Therefore, the envelope definition in Equation 11 can be appreciated by 
noting that the elementary component of X(t) is 90 deg out of phase with the 
elementary component of X(t) and 


Icos [a )t + 5 (/, co) + t//(co)] I = 1 


whenever 


sin [co/ + 8(t, cj) + i//(co)] = 0 


The autocorrelation function RJJxJq) of X(t) (Equation 10) can easily be 
shown to be equal to the autocorrelation function R^ti,t 2 ) of X(t) (Equation 
4). 

A A 

The correlation function E[X(t 2 )X(t\)] between X(t 2 ) and X(t{) can be 
shown as 


EiXitJcitJ] = /~ U(f,, 0 ))l L4(f ? ,a>)lsin Mf, -f 2 ) 
+ , w) - 5(r 2 , w)] G(a))du 


U3) 


A A 

Therefore, the covariance matrix A of X(t 2 ), X(t 2 ), X(t i) and X(t x ) can be 
written as 


A = 


( o 2 ) 2 

0 

a i<V> 

o,o 2 X 

0 

<**>* 

- a i a 2 X 


<w 

-o,o 2 X 

(ffj) 2 

0 

°l°2 X 


0 



(14) 


where o-j =o-(0) and cr 2 = cr(t 2 ) are the standard deviations of X(t\) and X(t 2 ), 
respectivelv, and p and A are nondimensional quantities as follows: 


/ 

‘’.t, 

% 


P ss (P l T 1 (o 2 T i Iq \A(t v <j)I \A(t r w)lG(w)cos Mr, -t 2 y 
+ 5(/ 1 ,u>)-5(/ 2 , co)] du 

^ = (^ 1 )" 1 ( a 2 )” 1 4“ •*(*!» <*>)• l4(f a ,w)lC(w)sin [w(t l -t 2 ) 
+ 6(/ 1 ,co)-6(r 2 , cj)] dcj 


( 15 ) 
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With the aid of the definition of envelope given in Equation 11 and the 
c ovariance matrix A given in Equation 14, the probability density function 
f t (x;t i) of V(ti) and the joint density function f vc (x,y;UJ i) of V(f 2 ) and V(fj) 
can be shown as 


[, x/(o ,) 2 ] exp (-x 2 /2(a 1 ) 2 ] 


(16) 




xy 


2/1 _ 2 \ 0 


( p . r^n 1-r 2 ) 


xvr 


OjOjO-r 2 ) 


X exp<- 


x 2 (a t ) 2 +y 2 (o 2 ) 2 

2(a 1 ) 2 (c7 2 ) 2 (l-r 2 ) 


(17) 


where I 0 [D] is the zero-order modified Bessel function of the first kind and 


r- [p 2 +X 2 ] 1/2 


(18) 


The moment function of V{t { ) and V(t 2 ) can be shown as 


v/2 

E { r ( tl ) r >( t 2 )) = [2(o t ) 2 ] [2(Oj) 2 ] 

x r (' +?)*'. (‘I 



(19) 


where T(*) and denote the gamma function and the hypergeometric 
function, respectively. 

Furthermore, the envelope crossing rate N(t\,a) of a barrier level a at 
time t\ can be shown as 


N(t v a) = 


aJK pIA + CWllf 

— 1 — «xp { 

vsfo.j* L 2 (°,) a J ( 


X exp 


i \fiiCct 
2V2 A(o,) J 


(20) 


' 2-2 


C‘a 


8(o, ) 4 A 


erf 
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where 


A = V [(X, ) 2 /(o, ) 2 ] - [C 2 /4 (o, ) J ] - c \ 

Xj = /“ U(/|,uJI 2 C(coXw + 6yrfw, /= 1, 2 I 

S' = d&(t l ,w)ldt l > (21) 

C = 24“ [ U((, , cj)I dU(r, , w)l/3f, lG(w)d<o = 2£W , », )] \ 

(2 - [3 L4(r, , cj)I/ 3/, ] 2 C(o))do) / 

When A(t,(o) is a slowly varying function of time t, 

C 2 /(o, ) 2 « A 

e« a 

5 r « co 

A-Xj-KX,) 2 /^,) 2 ] 


Therefore, Equation 20 becomes 

/V(f,,<*)=»— exp [-^(o,) 2 ] (22) 

V^(a, ) 2 

It can easily be shown that the nonstationary envelope process V(t) 
defined in Equation 11 degenerates into the stationary envelope process of 
Cramer and Leadbetter (Reference 5) when X(t) is stationary; in which case, 

A(t, cj) * 1 


50,, w)* 0 

and Equations 16, 17, 19, and 20 degenerate into the formulas associated 
with the stationary random processes. 

Consider a linear time-invariant system characterized by an impulse 
response function h(/), or by a frequency response function //(«), under the 
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excitation \(t) given in Equation 1. The response process Y(t) can also he 
expressed as a spectral representation (e.g. Reference 7): 


iut _ 


Y(t) = 12 M(t, a)e' dX(u) 

(23) 

, # -/U>T 

M(t, co) s f* h(j)A(t - t, co)e dr 

(24) 


where the assumption has been made that the system is at rest initially. 
Therefore, the envelope V y (t) of the response process Y(t) is as follows: 


vjt) = [r 2 (,') + rV)l 1,2 


(25) 


A 

in which Y(t) and Y(t) are in the forms of Equations 6 and 10, respectively, 
with A(t,<o) replaced by M(t,to). 

First Excursion Probability and N^rr.erical Example 

The probability that the structural response exits the safe domain of 
performance for the first time within the time interval [0,T] is called the first 
excursion probability. The first excursion probability of a structure under 
random loadings has received considerable attention among engineers 
because of its direct relation to the structural safe performance and design. 
We shall discuss here the applications of nonstationary envelope statistics 
derived previously to the prediction of the first excursion probability of 
structures under random disturbances. 

Let h(t,a) be the failure rate at time / of a gaussian random structural 
response X(t) with mean zero and threshold levels ±a. The first excursion 
probability F(T,a ) of the structure within the time interval [0,71 can be 
written as 


F(T, a) = 1 - exp [-£ h(t, a )dt ] (26) 


The simplest approximation, called the Poisson approximation, assumes 
that the events of crossings of the threshold level are independent. Hence, 


h(t.a) *2v\t,a) (27) 

in which v+(t,a) is the upcrossing rate of the level a at time t by X(f). The 
expression of v 4 (f,a) has been given, for instance, in References 7 and 9. 
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Another approximation is to replace the upcrossing rate v + (t,a) by the 
envelope crossing rate N(t,a) and assume that the events of envelope 
crossings are independent, i.e.. 


h(t,a)s*N{t,a) 


(28) 


where N(t,a ) is given by Equation 20 or 22. 

The approximation obtained in Reference 10 can be generalized to the 
nonstationary random processes as follows: 


Kt, a) w + (f. a) *- ex P 

1 - lv + (r,a)/v\t, 0)] 


( 29 ) 


Furthermore, the approximations based on the extreme point process 
(References 1 and 2) can also be generalized to the nonstationary random 
processes, since the joint density function (Equation 17) of the envelope 
process has been derived. It is noted that the establishment of nonstationary 
envelope statistics (Equation 17) sheds considerable light on the difficult 
problems of nonstationary first excursion probability. Further applications of 
envelope statistics based on the extreme point process (References 1 and 2) 
will be reported in the near future. 

The structural response to strong gust loadings or earthquake excitations 
discussed in Reference 9 is considered here. Using the definition of envelope 
established previously (Equation 25) and the simulation method proposed in 
Reference 8, the structural response X(t) and its envelope V(t) are simulated 
on a digital computer. Sample functions of X(t) and V(t) are given in Figure 
1. It is observed from Figure 1 that the envelope practically passes through 
all the peak values and the absolute values of all troughs of X(t). 

The first excursion probability based on (1) the independent level crossing 
(Equation 27) of X(t), (2) the independent envelope crossing (Equation 28), 
and (3) the generalization of Reference 10 (Equation 29) are summarized in 
Table 1 and plotted in Figure 2. Previous results obtained in References 11 
and 12 are also given in Figure 2 for the purpose of comparison. 


It is observed from Table 1 that the approximation based on envelope 
crossing does produce a slight improvement over Equation 27 when the 
threshold level oc is not too high. For the high threshold level, however, 
Equation 28 results in a worse approximation. This p\^nomenon is well 
known in the stationary random process when the average clump size is 
smaller than unity. It can be observed from Figure 2 that the approximation 
of Equation 29 gives results close to the simulation results (Reference 11). 
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TIME, $ 

Figure 1. Sample functions of structural response and its envelope to nonstationary 

random excitation 


THRESHOLD LEVEL a, cm 

0 0.25 0.50 0.75 1.00 



0 0.1 0.2 0.3 0.4 

THRESHOLD LEVEL a. In. 


Figure 2. Probability of structural failure vs threshold level 
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DISPLACEMENT, cm 


Table 1. Probability of failure 





Barrier levels 



Approximation 

0.38 cm 
(0.15 in.) 

0.51 cm 
(0.20 in.) 

0.64 cm 
(0.25 in.) 

0.76 cm 
(0.30 in.) 

0.89 cm 
(0.35 in.) 

Equation 27 

0.9998 

0.8958 

0.3776 

0.073 

0.0091 

Equation 28 

0.9961 

0.8551 

0.3900 

0.089 

0.0128 

Equation 29 

0.9931 

0.7496 

0.2678 

0.052 

0.0069 


Concluding Remarks 

A definition of the envelope of nonstationary random processes possessing 
evolutionary power spectra has been presented. The envelope definition 
enables one to simulate the nonstationary envelope process directly. The 
envelope statistics, such as the moment function, the joint density function, 
and the level crossing rate, relevant to the analysis of the first excursion 
probability, fatigue, and crack propagation, have been derived. Based on the 
envelope statistics, approximate solutions to the first excursion probability 
have been obtained. 
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Nonequilibrium Ionization Measurements in 
Hydrogen-Helium Mixtures 

L. P. Leibowitz 

Environmental Sciences Division 


Time-resolved emission measurements of several atomic line and continuum 
radiation channels have been made behind the incident shock wave of the JPL 
electric arc shock tube. Test times and nonequilibrium ionization times were 
obtained for shock velocities up to 2.5 X 10 4 m/s in a 0.2 H 2 - 0.8 He gas mixture. 
The shock-heated test gas was found to be free from driver gas contamination, and 
the test times were adequate to achieve steady-state conditions. An activation energy 
of 4 eV was obtained from the nonequilibrium ionization time measurements. 
Modifications to experimental technique to determine the effect of test gas impurity 
level on ionization time measurements are discussed. 


Introduction 

The existence of nonequilibrium flow during Jupiter entry could have a 
considerable effect on entry vehicle heating. An atmospheric en* / probe 
sent to Jupiter must withstand very large aerodynamic heating upon 
entering the Jupiter atmosphere. Heating rates based on the assumption that 
the shock-heated gases are in chemical equilibrium cannot be relied upon 
until the ionization kinetics of the shock-heated Jupiter atmospheric gases 
are understood. By measuring in the shock tube the time that is required to 
reach equilibrium behind shock waves in hydrogen-helium mixtures, the 
thickness of the nonequilibrium shock layer surrounding a Jupiter entry 
probe can be estimated. An electric arc shock tube built at JPL to simulate 
Jupiter entry conditions has been described by W. A. Menard (Reference 1). 
In the present investigation, the fluid flow produced in a hydrogen-helium 
mixture by the electric arc shock tube was studied and nonequilibrium 
relaxation times were measured. Measurements have been made of emission 
intensities for atom lines and continuum radiation as a function of time 
behind the incident shock wave. Characteristic ionization times and flow 
test times have been obtained as a function of shock velocity and initial 
pressure. These results are discussed and future plans are presented. 
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Experimental Technique 

The emitted radiation from important species present in the H 2 -He test 
gas and in the helium driver gas has been simultaneously measured as a 
function of time behind the incident shock wave of the JPL electric arc 
shock tube (Reference 1). Spectral intervals were selected in order to follow 
the chemical reactions in the test gas, to determine the length of the slug of 
test gas, and to determine whether there was any mixing between the test 
gas and the driver gas. It was necessary to determine spectral intervals for 
continuum measurements that were free from impurity radiation and to 
identify the strongest radiators in the driver gas. Spectrograms of the time- 
integrated emission from a shock-heated 0.2 H 2 - 0.8 He gas mixture were 
obtained using a modified Baush and Lomb grating spectrograph with 2.0 X 
10" 6 (20 A /mm) dispersion. A series of spectrograms were taken over the 
wavelength range 250 to 900 nm (2500 to 9000 A) using Polaroid Type 107 
film for the shorter wavelengths and Polaroid Type 413 film for wavelengths 
larger than 600 nm. Lines from hydrogen and helium were identified on the 
spectrograms; however, the most intense radiation was from metal atom 
lines due to impurities in the driver gas. 

Time-resolved emission measurements were made of five spectral 
intervals. Channels were selected to observe radiation from the line, 
atomic hydrogen continuum, a helium line, a portion of a hydrogen 
molecular band, and an atomic line from chromium, which is an impurity 
introduced into the driver gas by the vaporization of a stainless-steel trigger 
wire. The spectrometer was a Jarrell Ash Model 75-000 with a 0.75-m focal 
length and a dispersion of approximately 1.0 X 10 -6 (10 A /mm). The 
spectral interval of each channel and a schematic drawing of the optical 
system are shown in Figure 1. The five photomultipliers used were RCA 
Type 1P28A, and the photomultipliers and electronic circuits were each 
tested and found to have a response time of less than 0.2 /u.s. The time 
j ^solution of the optical system was less than the photomultiplier response 
time for the shock velocities studied in this investigation. 

The test gas was a 0.2 H 2 - 0.8 He mixture of Prep grade gases prepared 
by Matheson, Inc. The shock tube was pumped down to approximately 3.99 
X 10~ 2 N/m 2 (3 X 1(H torr) before each run. 

Results 

Shock tube runs were made over the range of shock velocities, U s , from 
1.6 X 10 4 to 2.5 X 10 4 m/s and for initial pressure of 33.3 and 133 N/m 2 
(0.25 and 1.0 torr). The output signals from the five photomultipliers were 
recorded on three Tektronix 555 oscilloscopes. Oscillograms obtained from a 
typical shot are shown in Figure 2. The signals from all channels reached a 
maximum slope rapidly and then asymptotically approached an equilib. ii m 
level without an intermediate overshoot. The He line channel and H 2 band 
channel produced intensity profiles similar to the H continuum channel. It 
was concluded that H continuum was the primary source of radiation for 
these channels. The Cr signal profiles were similar to the continuum channel 
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SPECTRAL INTERVALS 


CHANNEL 

X , nm 

AX , nm 

SPECIES 

1 

425.4 

0.5 

Cr 

2 

447.1 

0.5 

He 

3 

462.8 

1 .0 

h 2 

4 

486.1 

0.5 


5 

514.5 

1.0 

H + + e 


Figure 1. Schematic diagram of optical system 


until arrival of the contact surface, at which point there was a rapid and 
large increase in the signal level. The sudden rise in the Cr signal occurs at 
the same time as the sudden decrease in the signal. This marks the 
interface between the driver gas and the test gas. The decay of the signal 
over several microseconds indicates that there is a mixing region between 
the driver and test gases. It was concluded from the Cr and signals that an 
impurity-free slug of test gas with measurable test time was formed in the 
JPL electric arc shock tube. Test times obtained for initial pressures p\ = 
33.3 N/m 2 and p\ * 133 N/m 2 are shown in Figure 3. The curves in the 
figure are for test times calculated assuming constant values for test slug 
length l t ; reasonable agreement between calculated and experimental data 
can be seen. These test slug lengths are considerably smaller than the shock 
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tube diameter of 0.152 m and are several times smaller than the slug lengths 
that can be obtained from the test time measurements of Reference 1. This 
may be due, in part, to the sensitivity of the Cr channel in resolving the 
arrival of the driver gas. 

Both the H continuum and the line signals showed similar rise times. 
The rise time was defined as the time at which a line drawn through the 
maximum slope of the signal intersected the plateau signal level. The 
quantity (H 2 )° T jon is plotted as a function of 1/T 2 in Figure 4, where r ion is 
the rise time of the continuum signal in shock fixed coordinates, (// 2 )° is the 
frozen molar concentration of H 2 immediately behind the shock wave, and 
T 2 is the corresponding temperature. A reasonable fit to the p\ = 33.3 N/m 2 
data can be made by a straight line whose slope gives an activation energy of 
4 eV. The values of r ion measured with p\ = 133 N/m 2 were found to differ 
from the values measured with p ] = 33.3 N/m 2 by a factor approximately 
equal to the ratio of the initial pressures. Thus, r ion was found to be 
independent of the initial pressure. 


SHOCK SPEED PHOTOTUBE 
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Figure 2. Photomultiplier signals as a function of time for typical shot 


I 


16 












Conclusions 


Time-resolved measurements of emission intensity of the line and Cr 
line radiation have shown that a test slug free from driver gas impurities is 
formed in the JPL electric arc shock tube. Test times which allow 
nonequilibrium and steady-state measurements to be made were obtained 
for shock velocities up to 2.5 X 10 4 m/s. 

Both the low value for the activation energy and the insensitivity of the 
ionization time to initial pressure indicate that the collisional ionization of 
hydrogen is not the dominant process for the intensity profiles. Although the 
activation energy for H 2 dissociation (4.5 eV) is close to the experimental 
activation energy, the decomposition time is expected to be 100 times faster 
than the measured ionization times. Comparison with ionization processes 
in other gases leads one to expect activation energies closer to the 13-eV 
ionization energy of hydrogen and r ion to be inversely proportional to initial 
concentration. However, experimental measurements of ionization rates of 
noble gases have been found to be sensitive to impurities in the part-per- 
million range. In the runs of this experiment the impurity level in the initial 
test gas was approximately 880 ppm. Therefore, the domination of 
measured intensity profiles by impurity ionization processes is a definite 
possibility. 

Several modifications to the experiment'd technique are being made in 
order to determine the effect of impurity level on ionization time 
measurements. The pumping system is being modified so that a diffusion 
pump and liquid nitrogen cold trap will be used at each end of the shock 
tube. The test gas filling system will be changed to permit continuous 
flushing of test gas down the shock tube, thus minimizing the effects of 
absorption and desorption of gases on the shock tube wall. Finally, a test gas 
mixture made from ultra high purity grade gases and with an analysis of all 
impurities is being obtained. 

The measurement of r jon as a function of U s and p\ w ill be repeated to 
determine the effect of the improved purity of the system. If these 
measurements indicate that impurities are no longer a problem and that the 
ionization times have not become larger than the available test time, an 
improved spectroscopic technique is planned. By making both absolute and 
restive continuum intensity measurements across the Balmer series limit, 
the electron concentration and electron temperature can be obtained as a 
function of time without the reliance on collisional equilibrium assumptions 
required of other spectroscopic techniques. 
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On the Possibility of Earth Re-entry Simulation 
of Shallow Angle Jupiter Entry 

A. D. McRonald 
Environmental Sciences Division 


Possible Earth re-entry simulation of shallow angle (3-30 deg) Jupiter entry has 
been investigated in terms of four parameters of the bow shock layer ahead of a 
blunt vehicle: peak (equilibrium) temperature, peak pressure, peak inward radiative 
flux, and time-integrated radiative flux. The comparison shows that simulation 
ranging from fair to good can be achieved, generally the easiest (lowest Earth re- 
entry speed) at steep Earth re-entry, in the Earth entry speed range of 15-22 km/s, 
for both the Jupiter “nominal” and “cool” atmospheres. Increasing Earth speed is 
required, generally, for increasing Jupiter entry angle, and for temperature, radiative 
flux, time-integrated flux and pressure, in that order. It appears that a meaningful 
simulation test could be done using a launch vehicle with the speed and payload 
capability of the Titan IIID/Centaur/BII. 


Introduction 

The possibility of simulating the major features of entry from near- 
parabolic inertial entry speed into the atmosphere of Jupiter by means of an 
Earth re-entry test is of considerable interest on two main accounts: 

(1) As a “calibration” point for several aspects of Jupiter entry probe 
design, including: 

(a) Theories and models for the very high enthalpy hypersonic flow 
field with strong coupling of intense radiation, convection, and 
ablation. 

(b) Ablation materials response when subject to high radiative flux 
(of order of 100 kW/cm 2 ), large gasdynamic mass flux over the 
body, and large “inertia” forces due to violent deceleration. 

(2) As a prior-to-launch validation of design of a prototype Jupiter entry 
capsule, including trajectory, aerodynamic stability, structural 
integrity, heat shield, communications, instrument performance, etc. 
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Analysis 

Simplifying Assumptions 


For this first analysis, v irious simplifying assumptions hive been made, 
including: 

(1) The shock layer is in thermochemical equilibrium, is uniform in 
composition, and is an isothermal parallel slab of uniform thickness 
(presupposes a fairly large blunt vehicle and fairly nigh gas density; 
good approximation near peak conditions during entry; facilitates 
temoerature and radiative flux calculations). 

(2) The shock layer thickness is constant during the entry (ignores change 
in body (nose) radius due to ablation and change in gas density ratio 
over the shock; facilitates radiative flux calculations). 

(3) The shock layer pressure is Pi(V rd ) 2 , here p | is the undisturbed 
atmospheric density, and V rc | is the speed of the vehicle relative to 
the atmosphere. 

(4) All other factors relating to the comparison of Jupiter and Earth 
entry, e.g., diffusion, chemical reaction, and ablation, are eithei 
similar in effect or minor. 

Trajectories 

The trajectories considered for Jupiter entry were for a ballistic 
coefficient 1 of 154 kg/m 2 , an inertial entry speed of 61 km/s (slightly 
hyperbolic), entry angles of from 3 (near ballistic skip-out) to 30 deg, and a 
“downwind” equatorial entry (giving V rd at entry in range 48-50 km/s). For 
Earth, similar entry at all entry angles, the same ballistic coefficient, and 
inertial entry speeds of 12, 14, 16, 18, 20, and 22 km/s were evaluated. A 
two-degree-of-freedom (motion in a plane) computer program was used to 
calculate approximate trajectories. 


Radiative Flux 

Data on the emergent radiative flux from the parallel, isothermal slab of 
hot gas of uniform thickness, as a function of slab pressure and temperature, 
have been published (Reference 1 for air. Reference 2 for the Jupiter 
nominal and cool atmospheres). The JPL thermochemistry program 
(Reference 3) was used to calculate shock layer temperature, pressure, and 
gas density ratio across the shock. The slab (shock layer) thickness 8 was 
evaluated from the Lighthill formula: 


1 The value is one for which many prior trajectory calculations have been made by the author, 

and corresponds to a 45-deg cone of base radius 1 m and an average specific gravity of K. 


20 





where R is the vehicle nose radius, € = p\/ p 2 , and p 2 is the shock 
density. For convenience of calculation, the 5 values for Jupiter gas ai 
were taken as 5 cm and 3 cm, respectively; the body is sketched, and i 
of R, 8, and p 2 l p\ are given in Figure 1. 

Results 

The absolute levels of peak shock layer pressure, temperature, ii 
radiative flux, and integrated flux for the entry bodies considered are s 
in Figures 2, 3, 4, and 5, respectively. From these, the Earth entry spec 
angle to simulate a given Jupiter entry angle as regards the abovt 
parameters can be calculated. Figures 6 and 7 show Earth entry cond 



Figure 1. Type of blunt entry body considered in radiative flux calculation 
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to simulate 10-, 20-, and 30-deg Jupiter entry into the nominal and cool 
atmospheres, respectively. The 90-deg Earth entry speeds to simulate the 
various parameters for 10-, 20-, and 30-deg entry into the Jupiter nominal 
and cool atmospheres are shown in Table 1. 

As a prelude to discussion of simulation, it is appropriate to comment on 
performance of available launch vehicles. The feasibility of an Earth entry 
test depends on launch-vehicle capability, which is outlined in Figure 8 (data 
from Reference 4). 



Figure 2. Peak shock layer pressure versus entry angle 
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Figure 4. Peak inward radiative flux versus entry angle 


approach constancy as increasing density brings the flux toward the black 
body limit. 

As regards the Jupiter entry angle to be simulated, it is noted that the 
pressure is very sensitive to angle, the peak radiative flux quite sensitive 
(particularly at low angles), and the integrated flux and temperature not very 
sensitive (low flux at low angle being compensated by longer entry flight 
time). 

The iaunch-vehicle payload capability is likely to impose serious restraints 
on Earth re-entry vehicle mass, and hence on R and m/CpA. Similarly, entry 
speed and weight appear to be limited, for practical reasons of cost, to about 
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Figure 5. Time-integrated inward radiative flux versus entry angle 


19 km/s and 50-kg payload, corresponding to a vehicle like the Titan HID/ 
Centaur/BII. It may then be difficult to construct an Earth re-entry vehicle 
with R = 60 cm and m/CpA * 154 kg/m 2 as calculated here, subject to a 
50-kg mass, and there is clearly a trade-off between mass and speed for best 
simulation. However, the conclusions below will be given in terms of an 
Earth vehicle of comparable size to the Jupiter vehicle. 


Conclusions 

The following conclusions can be drawn from Figures 6, 7, and 8, as 
regards Earth entry simulation of Jupiter entry. 
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Figure 6. Simulation of 10-, 20-, and 30-deg Jupiter entry by Earth re-entry: Jupiter 

nominal atmosphere 


(1) As the Jupiter entry becomes more shallow, less Earth entry speed is 
required for simulation of the four parameters. 

(2) For a given simulation, less Earth entry speed is required at steep 
Earth entrv than at shallow. 

(3) Shock layer temperature is simulated most readily, requiring Earth 
entry speeds in the range 15 to 17 km/s for the Jupiter nominal, and 
18 to 19 km/s for the cool atmosphere, at Jupiter entry angle of 10 to 
30 deg. 

(4) Inward radiative flux simulation can be achieved for 10-30-deg 
Jupiter entry into the nominal atmosphere (Earth entry speeds of 15 
to 18.6 km/s), but only for about a 10-deg entry into the cool 
atmosphere (Earth entry speeds of 18.8 to 21.0 km/s for 10-30-deg 
entry). 

(5) Time-integrated radiative flux can barely be simulated, requiring 
Earth entry speeds of 19.2 to 20.3 km/s for 10-30-deg Jupiter entry 
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Table 1. Earth re-entry speeds to simulate shallow Jupiter entry 


Parameter 


Earth entry speed at 90 deg, 

km/s 


Nominal atmosphere 

Cool atmosphere 

10 a 

20 a 

30 a 

10 a 

20 a 

30 a 

Temperature 

15.0 

15.9 

17.0 

18.0 

18.6 

19.0 

Pressure 

12.4 

19.3 

23.0 

18.0 

24.0 

— 

Flux 

15.1 

17.4 

18.6 

18.8 

20.2 

21.0 

Integrated flux 

19.2 

19.8 

20.3 

22.2 

22.2 

22.2 

a Jupiter entry 

angle, degrees. 









Figure 7. Simulation of 10% 20-, and 30-deg Jupiter entry by Earth re-entry: Jupiter 

cool atmosphere 
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Figure 8. Performance of present and proposed launch vehicles 


into the nominal atmosphere, and about 22.2 km/s for similar entry 
into the cool atmosphere. 

(6) Shock layer pressure can be simulated for 10-20-deg entry into the 
nominal atmosphere (Earth entry speeds of 12.4-19.3 km/s) and for 
only about 10-deg entry into the cool atmosphere (Earth entry of 18.0 
km/s for 10 deg, 24.0 for 20-deg entry). 

(7) Prospects of simulating two, three, or four of these parameters 
(temperature, pressure, flux, and integrated flux) for a specific case, 
by proper choice of ballistic coefficient, nose radius, and Earth entry 
angle and speed, are good for Jupiter entry in the range 10-20 deg for 
the nominal atmosphere, and 10-15 deg for the cool atmosphere. 

(8) A meaningful Earth re-entry simulation of shallow Jupiter entry 
appears to be possible using a launch vehicle about the scale of the 
Titan IIID/Centaur/BII (19 k*n/s at 50 kg, 16.5 km/s at 500 kg). 
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Measurement of the Power Spectral Density of 
Phase of the Hydrogen Maser 

A. Sward 

Telecommunications Division 


Hydrogen masers are being developed by JPL to satisfy the needs of the Deep 
Space Network for improved frequency standards in order to achieve better timing 
accc.-icies, reduced doppler residuals, and lower phase noise in the tracking system. 
Performance data for the JPL hydrogen m isers is given in two forms: tirne-domain 
performance, useful for analyzing time-keeping systems and doppler systems; and 
frequency-domain performance, useful for analyzing tracking systems.. Particular 
attention is given to frequency-domain performance, as obtained from both time- 
domain measurements and direct measurements. 


The power spectral density of phase for the JPL hydrogen masers 
(References 1, 2, and 3) has recently been measured and found to be in close 
correlation with that which was calculated using time-domain stability 
measurements. This has allowed the 3-dB bandwidth of the hydrogen maser 
to be calculated as well as the frequencies at which the various types of 
noise are predominant. 

A plot of the measured standard deviation of the average fractional 
frequency departure versus averaging time for the hydrogen maser is shown 
in Figure 1. The flat region of the plot, which occurs for averaging times of 
50 s and greater, is caused by a random noise called flicker-of-frequency 
noise, and is characterized in the power spectral density of phase by a /" 3 
slope. The portion of Figure 1 that varies approximately as 1/r for 
averaging times of less than 50 s is called flicker-of-phase noise. This noise 
would have a spectral density of phase varying as/" 1 . 

Before describing the equations that are used to convert from the time- 
stability plot to the power density spectrum of phase, it is important to 
define several terms. First, the term <r 2 (N,Tj) is a sample variance of N 
measurements of the time-varying frequency, f{t), each of duration r, spaced 
every T units of time. The units of <r 2 (jV,7» are rad 2 /s 2 . The standard 
deviation of the average frequency departure is divided by a> 0 to obtain the 
average fractional frequency departure. The variance of frequency fluctua- 
tions for flicker of frequency exhibits a logarithmic divergence as the 
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number of measurements increases. It is thus necessary to fix the number of 
measurements in order to obtain meaningful data. A convenient measure, 
called the Allen variance, is obtained by fixing .V at 2, and by modifying the 
usual sample-variance definition by the factor .V/(.V-1). It can be shown that 
the optimum value (Reference 4) for .V is, in fact, very close to 2. Normally, 
the mean sample variance for a large number of sample-variance measure- 
ments is taken in order to reduce the uncertainty for such a low number of 
measurements .V. 

In order to calculate the spectral density of phase, the following 
relationships are needed: For flickcr-of-freqtiencv noise (Reference 5>, 


0 (2tt) 2 / 2 In 2 


0 ) 


where S J[f) = power spectral density of phase in rad 2 'Hz. For fiickcr-of- 
phase noise (Reference 5), 


<7 (f\ - [a 2 (2,T,T>] 

0 {J) (0{3[2 + lr (r W „)] - In 2} 


( 2 ) 


where <o n = noise bandwidth in rad/s. An important point should be noted 
here. In the time-domain plot, frequency stabil ity va ries as 1/r for white 
phase noise, whereas the stability varies as , In r/r for flicker-of-phase 
noise. 

White phase noise (J° behavior in the frequency domain) produces a 1/r 
behavior in the time-domain plot of <r versus r, while flicker-of-phase noise 
(f~ l behavior in the frequency domain) produces a y/U\ t/t dependence 
(Reference 5). These two types of noise are difficult to distinguish in the time 
domain because o r measurement dispersion. A direct measurement in the 
frequency domain is often more useful to determine which type of noise is 
present. Such measurements have, in fact, verified the existence of fiicker-of- 
phase noise. 

The time-domain plot of the average fractional frequency departure, 
(t{2,t,t)/io 0 versus averaging time r for the hydrogen maser, is shown in 
Figure 1. The power spectral phase density calculated from Figure 1 is 
shown in Figure 2, with the measured values superimposed. Equating and 
solving Equations 1 and 2 results in a solution for crossover frequency 
between the two types of noise. This value is 1.59 X 10~ 2 Hz at a center 
frequency of / 0 * 100 MHz. Direct measurements at the device 100-MHz 
output have verified the validity of the spectral plot within 1.4 dB over the 
frequency band 4 to 100 Hz on each side of the carrier. In addition, at the 
10- and 1-MHz outputs, obtained by frequency division of the 100-MHz 
signal, measurements have degradations of only about 1 dB in the power 
density spectrum. It was found that at the 1-MHz output, white phase noise 
became the predominant noise above 20 Hz from the carrier, which, if 
extrapolated back to the 100-MHz output, corresponds to a crossover 
frequency of about 200 kHz. 
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The 3-dB bandwidth of the hydrogen maser at its output frequency of 
1420 MHz was calculated by extrapolation to be 6.4 X 10~ 4 Hz from the 
power spectral density plot. 
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Reduction of Near-In Sidelobes Using Phase 
Reversal Aperture Rings 

W F. Williams 


Telecommunications Division 


Spacecraft antenna;: having high gain and extremely low sidelobes will be needed 
for communications satellites. Low sidelobc performance is needed to isolate beams 
of multibeam systems using the same carrier frequency (frequency re-use). Because 
of a limited allowable spectrum for transmission, frequency re-use may be necessary. 
A technique for reducing near-in lobes by cancellation is described. This technique 
takes a small portion of the radiation from the antenna aperture and generates the 
near-in lobes, which are are then fed out-of-phase relative to the main signal. Results 
of sample cases indicate that the first three lobes can be nearly eliminated at a 40% 
reduction in aperture efficiency. 


Introduction 

A current NASA-sponsored work program has as its objective the 
development of antenna techniques that may be useful in high-gain satellite 
antennas used for Earth communications. One of the most important 
characteristics of such future antennas will probably be extremely low 
sidelobes. This is because of the expected use of the same ecu i iei frequency 
channel simultaneously on several different antenna beams, each carrying 
different information. To avoid interference, such beams should not overlap. 
However, interference also exists when the sidelobe of one such beam 
occurs within the main lobe of another. The power transmitted within this 
sidelobe, as dictated by its gain, will determine whether it is detected by an 
Earth station and cause interference. 

Sidelobe size and structure are the result of the illumination function on 
the antenna aperture. Uniformly illuminated circular apertures result in first 
sidelobes of about -17 dB. When the illumination is reduced at the aperture 
edge (tapering), sidelobes are reduced. Interfering structures in front of an 
aperture (feeds and struts for a paraboloidal reflector) will change this 
distribution and also raise sidelobes. A typical reflector antenna with 10-dB 
taper and average blockage may have first sidelobes of about -20 dB. More 
extreme tapers can be used, and, with careful attention to blockage, sidelobe 
structure can be improved. But it is generally conceded to be doubtful that a 
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-30-dB sk\?lobe level could be reached by edge-taper schemes— probably a 
-26-dB number is more likely. 

A required level of sidelobe reduction has not been established for any 
specific system as yet. However, there are indications (Reference 1) that the 
levels mentioned above may not be adequate and that means other than 
edge-taper might be needed to obtain even lower sidelobes. 

The lobes of importance are only the first, second, and possibly the third, 
since beyond that point the lobes will be looking into regions of no interest 
for frequency re-use. 

The following sections describe a technique for lowering near-in 
sidelobes, bat at the expense of gain and an increase in the lobes that occur 
at greater angles from the main beam. 

Sidelobe Cancelling 

An approach to this problem is to construct a second antenna whose 
radiation pattern is similar to only the sidelobe structure of the main 
antenna. If this second antenna is then excited out-of-phase relative to the 
main unit and of proper magnitude, the sidelobes will be eliminated or 
greatly reduced. The example chosen for the cancelling antenna for 
discussion in this article is a selected ring or rings within a circular aperture. 

Consider a circular aperture as the radiating source of the spacecraft 
antenna. It might be a circular paraboloid, or the face of a phased array, 
horn, or lens antenna. As a parabola or lens, the illumination would normally 
be tapered for 20-dB sidelobes; as an array, it could be uniform for 17-dB 
first sidelobes. Further, consider this aperture to be made up of a complete 
set of narrow adjacent concentric circular rings. Each ring may be assumed 
to be uniformly excited, the intensity being a function of its radius and the 
taper on the overall aperture. Each of these rings is then conside r • is a 
separate radiating aperture. The illumination on this sub-aperture is uniform 
around the periphery and zero throughout the rest of the area, which gives 
rise to a multilobed interference type of pattern— the larger the ring radius, 
the narrower these lobes will be. The sum of the patterns of all rings is the 
total pattern of the aperture. This suggests the way to reduce the inner 
sidelobes. A particular ring region will have its first interference lobe at the 
same (or near) position as the first lobe of the total antenna pattern. If this 
region of the aperture were then excited at 180-deg phase reversal, it would 
tend to cancel this first main lobe. Magnitude would be set by choosing the 
proper ring region and varying its width. 

Sample Calculations 

Figure 1 represents a coordinate system and the quantities that make up a 
circular aperture consisting of rings of equal width. A good approximation to 
a tapered illumination from a horn-type feed for either paraboloid or lens is 
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where r 0 is the aperture radius, and a defines the taper magnitude at the 
i aperture edge. As a goes to zero, 100% efficient uniform illumination would 

be obtained. 
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When the aperture is made up of many rings, each ring can be considered 
as uniformly illuminated, at a value 

: f * t 

v^; ' 

watts/cm 2 (2) 



A 

■a- 


P(r) = exp 




> » 



for the particular ring indexed by £ . The far-field pattern of such a ring can 
be considered as the sum of contributions from its smaller segments, A A, 
identically excited: 
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where k g = integer value of 



/eg is chosen in this manner so that the segments that make up the rings will 
be of nearly the same area throughout the aperture. 


The total pattern is the sum of all rings, 



To demonstrate the performance of such an antenna, a 50-wavelength 
aperture with 10-dB taper was assumed, and ring patterns were calculated 
for the outer half. The inner rings have broad patterns and are not useful in 
this cancelling technique. The total number of rings was set at 36; patterns 
of some .of these are presented in Figure 2, which shows the relative 
spacings of their lobes. Boresight is at 0 deg. The total aperture pattern is 
then shown in Figure 3; the first, second, and third sidelobes are at 2.1, 3.25, 
and 4.4 deg. Observation indicates that ring 24 (% of the distance to the 
edge) generates a lobe in such a position as to cancel the first sidelobe. 
However, it does not affect the second lobe, but increases the third lobe of 
the main pattern. Further, ring 30 cancels the second lobe, while ring 31 
cancels the third, but does not aid in cancelling the first sidelobe. Selective 
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Figure 2. Ring radiation patterns 

combinations must therefore be tried to obtain a proper combination to 
reduce the first three lobes. Such a result is indicated in Figure 4. The first 
three sidelobes of the main pattern are indicated (main lobe is not shown), as 
well as the contributions of four selected rings. The result of reversing the 
phase of the selected rings is near-cancellation of the sidelobes from 2 to 5 
deg off boresight, with a lobe coming up beyond 5 deg. Figure 5 shows the 
resulting radiation pattern, to be compared with that shown in Figure 3. 
Sidelobe radiation lies below -37 dB out to about 5 deg, beyond which point 
the lobes increase well above the normal pattern, but still remain below -28 
dB. Other similar combinations can be chosen with slight variations in the 
result. 

Attempts to achieve significant sidelobe reduction using only a single ring 
did not succeed. Although the first lobe was reduced, other lobes increased. 

The technique described could be easily implemented with a waveguide 
lens antenna by placing 180-deg phase shifters in selected ring areas about 
the lens center. 
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Figure 2. (contd) 

Conclusion 

Phase reversal of selective aperture zones can effectively cancel given 
sidelobes in the radiation pattern of a circular aperture far beyond the 
sidelobe reduction obtainable by edge-taper techniques. The technique 
necessarily reduces the on-axis antenna gain because of the reduced 
illumination efficiency. In the example given, this gain was reduced by about 
2.2 dB in cancelling the first three lobes. 
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Figure 4. Illustration of the subtraction of the contributions from four rings 
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Development of a 20-W Solid-State S-Band Power 

Amplifier 

A. G. van der Cape lien 

Telecommunications Division 


As an alternative to the use of traveling-wave-tube amplifiers in spacecraft with 
long-life mission requirements, JPL is developing a solid-state 20- VV S-band power 
amplifier. Traveling-wave-tube amplifiers have limited reliability because of the 
relative short life of the cathode and the complexity of the power supply. A soiid- 
state amplifier with an output of 20 VV at 2295 MHz, a dc/RF efficiency cf 38%, and 
a gain of 27 dB has been developed. This article describes the physical and electrical 
performance of the solid-state amplifier. 


Introduction 

A solid-state 2295-MHz prototype power amplifier is being developed in 
support of anticipated long-life missions such as the Outer Planet Grand 
Tour missions and the Space Station program. Traveling-wave-tube (TWT) 
amplifiers are not suitable because of their limited reliability and short life. 
Presently available microwave power transistors make solid-state power 
amplification at S-band feasible. The sobd-state power amplifier is to be 
designed to operate continuously for 12 yr. Design goals for this amplifier 
are summarized in Table 1. 

Circuit Description 

AH circuitry is housed in one module. Figure 1 is a block diagram of the 
power amplifiei and Figure 2 shows the module. The amplifier consists of 
three driver stages. The first two stages employ MSC 3000 transistors and 
the third stage uses the MSC 3003. The final stage consists of four summed 
MSC 3005 transistors. All stages except the first operate Class C in common 
base configuration. The first stage is biased in Class A (common base). A 
single-stage amplifier is shown in Figure 3 and a schematic in Figure 4. 
Semilumped circuitry is employed in the amplifiers, consisting of two pairs 
of Johansori Type 6453 capacitors (in a back-to-back configuration) forming 
an L matching network. Interstage coupling is achieved by short pieces of 
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Table 1. Design goals for the 20- W solid-state S-band power amplifier 


Electrical: 

RF power output 

Frequency 

Efficiency 

Bandwidth 

Cain 

Spurious signal output 
Input impedance 
Load impedance 

Environmental: 

Life (combined effect of 
radiation, aging, etc.) 

Temperature 

Baseplate 


20 W 
2295 MHz 

30% minimum (goal 40%> 

>20 MHz (1 dB) 

28 dB 

>70 dB below rated output 

50 a nominal (VSVVR to be determined) 

50 a nominal (open or short 
circuit capability to be 
determined) 


12 yr (confidence = 0.95, 
power degradation < 1 dB) 

-40 to +65°C (power degradation 
to be determined) 

35°C (average temperature 
for reliability calculations) 


Physical (exclusive of thermal 
Daseplate): 

Weight 1.3 kg 

Size 12.5 X 22 X 4.5 cm 


50-0 microstrip lines on alumina substrates to fac; . late testing. The power 
splitter and summer are identical in design and consist of 50-0, 3-dB hybrids, 
which provide a port-to-port isolation ot approximately 30 dB and input-to- 
output insertion loss of 0.3 dB (for one splitter or summer). The final 
summing point is connected to a drop-in type circulator (0.3-dB insertion 
loss at 2.3 GHz), which serves as protection against inadvertent load 
removal. All amplifier stages are voltage-derated and operate off a dc supply 
voltage of 25 V. 

Assuming thermal resistances for the transistors to be as stated by the 
manufacturer, i.e., 

MSC 3005 - 8.5°C/W 
MSC 3003 - 17°C/W 
MSC 3000 - 45°C/W 

the highest calculated junction temperature occurs in the MSC 3003 driver. 
For an assumed baseplate temperature of 35°C, this transistor junction is 
calculated to operate at 96°C. This is within the goal of 100°C maximum 
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4 x MSC 3005 



SPLITTER 
-0.3 dB 


SUMMER 
-0.3 dB 


2295 MHz 
20 W 


Figure 1. Block diagram of 20- W solid-state S-band power amplifier 



Figure 2 . Breadboard of 20-W solid-state S-band power amplifier 
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Figure 4. Single-stage amplifier schematic 
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Table 2. Performance data for single-stage amplifier MSC 3005 a 


Output 
power, W 

Gain, dB 

Efficiency, 

% 

Supply 
voltage, V 

6.8 

9.3 

42 

25 

6.4 

9.0 

50 

25 

5.0 

5.6 

54 

25 

5.0 

4.8 

55 

25 

8.0 

7.3 

45 

28 

a Frequency = 2295 MHz: bandwidth = 100 MHz (3 dB). 


junction temperature set for this development. The MSC 3005 operates at a 
junction temperature of 86°C (under the same assumptions). 


Development Status 


Th 


e power amplifier breadboard shewn in Figure 1 has been developed at 
JPL. An output of 20 W at 38% dc/RF efficiency and 27-dB gain has been 
achieved with a dc supply voltage of 25 V. Bandwidth is 35 MHz (3 dB), 25 
MHz (1 dB). Data for the output stage (MSC 3005) are shown in Table 2. 


Four MSC 3005 transistors were radiation- tested. Table 3 summarizes the 
dose values that the transistors were exposed to (including background 
radiation). The transistors were not operating during irradiation. They were 
mounted in test fixtures where they remained throughout all RF and 
radiation testing. This procedure was followed to avoid a shift in the data as 
a result of displacement of the transistor relative to the circuit. All four MSC 
3005 transistors survived the radiation test without observable (resolution 0.2 
to 0.5 dB) degradation. 

RF breakdown tests have been conducted on the final stage. No 
breakdown of any kind was observed at pressures between 10.1 N/cm 2 and 
4.0 X 10 -7 N/cm- (760 torr and 3 X 10~ 5 torr). (Electron source was present.) 


Table 3. Radiation dose values 


Radiation 

type 

Energy, 

MeV 

Flux, 

particles/cm 2 -s 

Fluence, 

particles/cm 2 

Background 

Electrons 

2 

3 X !0 10 

1.4 X !0 13 

None 

Neutrons 

1.5 

1.4 X 10 10 

1.6 X 1G 13 

4.3 X 10 3 rad 

Protons 

140 

1.2 X 10 8 

9 ±2 X 10 12 

10 9 - 10 10 neutrons/cm 2 
5 X 10 4 rad 
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Four MSC 3005 transistors in final stage configuration are presently 
under life test; 1200 h have been accumulated without failure. Baseplate 
temperature is held at 35°C. 

Work has begun on the development of a microstrip version of the 20-W 
power amplifier with plans to achieve a wider bandwidth (200 MHz, 3 dB 
minimum). This amplifier will employ an RCA TA 7995, which has emitter 
ballasting for better reliability. 

Future plans include a trade study between efficiency and other require- 
ments such as bandwidth, life, and temperature. Radio frequency life testing 
will continue to gather more data. Also, failure mechanisms in microwave 
power transistors will be studied along with a transistor RF circuit model. 
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A Markov Chain Technique for Determining the 
Acquisition Behavior of a Digital Tracking Loop 

H. D. Chadwick 

Telecommunications Division 


Tracking loops have two characteristic modes of operation: acquisition and 
steady-state tracking. The steady-state behavior of such nonlinear tracking loops as 
the phase-locked loop has been the subject of considerable analysis. The acquisition 
behavior of a loop, the transition period between turning the loop on and the steady 
state, has resisted analysis for all but the simplest configurations. An iterative 
procedure is presented for determining the acquisition behavior of discrete or digital 
implementations of the tracking loop. The technique is based on the theory of 
Markov chains and provides the cumulative probability of acquisition in the loop as 
a function of time in the presence of noise and a given set of initial condition 
probabilities. A digital second-order tracking loop to be used in the Viking command 
receiver for continuous tracking of the command subcarrier phase has been analyzed 
using this technique, and the results agree closely with experimental data. Possible 
extensions of the technique include the analysis of continuous loops using discrete 
approximations. 


Introduction 

A digital second-order servo loop is to be used in the Viking command 
receiver for continuous tracking of the command subcarrier phase (Refer- 
ence 1). An initial coarse estimate of the subcarrier phase is provided by a 
separate circuit, and the tracking loop is expected to begin tracking with an 
initial phase error as large as the maximum error of the coarse estimate and 
a phase error rate of change due to the doppler shift. The subcarrier loop 
should be able to acquire lock and eliminate this initial error, thus yielding a 
small steady-state average phase error. The behavior of this loop during this 
initial acquisition is studied in this article. The cumulative probability of 
acquisition as a function of time in the presence of noise has been calculated 
for the loop by a Markov chain approach. 

Tracking Loop Model 

The subcarrier tracking loop being analyzed is illustrated in Figure 1. It is 
a typical digital second-order loop which is complicated by the presence of 



the limiter in the error signal path. Figure 2 shows the same loop reduced to 
its important details. The aceumulate-and-dump circuit for the error signal 
accumulates the values of 2 M phase samples and has the effect of only 
improving the signal-to-noise ratio per 2 M samples. Otherwise it may be 
ignored for analytical purposes. The effect of using the sign bit of this 
accumulator is equivalent to that of a limiter circuit. The subcarrier 
waveform, which is sampled in the phase detector, is basically a square 
wave, with a small linear region in the vic nity of the zero point. This 
waveform is illustrated in Figure 3. The phase quantization level or “bump” 
size is determined by the phase detector and the phase estimate accumula- 
tor. Sampling in the phase detector can be performed twice per cycle with 
appropriate sign reversal circuitry. This circuitry can also be ignored in the 
analysis. The phase detector can thus be considered to provide a signal that 
is linearly related to the phase error over a small region close to zero and 
saturates to a constant level beyond the linear region. A new phase estimate 
is provided once every M subcarrier cycles (one data bit period). 


Discrete State Space 

Using the model described in the previous section, there are two state 
variables inherent to the tracking loop. These variables have been assigned 
and labeled and x 2 in Figure 2, where X\ is the phase error signal and x 2 
the loop accumulator output. The state space f or the loop is thus two- 
dimensional and discrete because movement from one point to another can 
only take place in discrete jumps. The system is also discrete in time, and the 
movement in the state space is describable by difference equations. 

In the absence of noise, because of the nonlinearity presented by the 
limiter, the state space must be analyzed in two regions: the region in which 
Xj > 0 and the limiter output is positive, and the region in which xj < 0 and 
the limiter output is negative. In this noiseless case, the loop will follow a 



Figure 1. Subcarrier tracking loop 
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Figure 2. Equivalent subcarrier tracking loop 



Figure 3. Phase detector characteristic 


deterministic path from a point representing a given set of initial conditions 
to an equilibrium point. 

When noise is present, however, the path becomes a random walk in the 
x h x 2 plane. At each point in the plane there are two possible subsequent 
points. One such point is the same as the noiseless deterministic point and 
will occur if the limiter input signal (phase error) plus noise has the same 
polarity as the signal only. The other point may be considered the wrong 
point and occurs if the noise biases the limiter input to a polarity opposite 
that of the phase error. Thus, for every possible origin point in the plane, ' 
there are two destination points, and the transition probabilities can be 
computed from the noise probability density and the phase error amplitude. 
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Transition Probability Matrix 


The movement of the state of the tracking loop in the state plane can be 
considered as a random walk or Markov chain. If the transition probabilities 
from state to state are written in the transition probability matrix P, then, in 
matrix notation, 


P 


n + 1 



( 1 ) 


where p n is the vector of state probabilities after n iterations of the loop. 
Given the initial state probabilities p 0 , iterative application cf Equation 1 
will provide the state probabilities for any time thereafter. The boundary 
conditions for the state space are crucial to the application of this equation. 
In a steady-state analysis, it does not matter which subcarrier cycle the loop 
is presently tracking. For this purpose, the boundaries of the phase error 
variable at ±Vi cycle from the origin may be made reflecting. This is 
equivalent to considering the phase error modulo 2ir, which is common in 
phase-locked loop analyses. In this case a steady-state solution to the state 
probabilities can be obtained from the eigenvector equati " n 


p-Pp 


( 2 ) 


where p is the vector of steady-state probabilities. For acquisition, however, 
movement from one subcarrier cycle to another is considered failure to 
acquire and this probability is of interest. (This definition of acquisition 
without cycle slipping was a requirement of the Viking mission.) The 
boundary at = ±Vi is, therefore, made absorbing, and the probability of 
reaching this boundary is called the probability of failure to acquire. With 
these boundary conditions, the steady-state probability vector is zero for all 
stales except the boundary. The loop will eventually go out of lock with 
probability one. For the purposes of studying the acquisition behavior, a 
small region around the equilibrium point was defined as the lock region, 
and the cumulative probability of reaching this region was calculated as a 
function of the number of iterations of the Markov chain equation. 


Determination of Transition Probabilities 


Examination of the system diagram in Figure 2 shov's that the difference 
equations for the transitions that can occur are: 
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E(x 2 ) = K 2 sgn | f{x x ) + n] 


( 3 ) 


E(x x ) = E{6 )-x 2 sgn \f{x x ) + n] (4) 


where E is the forward difference operator, 6 is the phase angle input, E(0) is 
the change per time interval of the phase angle, and n is the sum of 2 M 
gaussian noise samples. The function /(•) is the phase detector characteristic 
of Figure 3 multiplied by the accumulator gain of 2 M. When 


f( x x ) + n > 0 


the signum function is + 1 and 


£•(**)= *2 

( 5 ) 

H\x l )°E(e)-x 1 -K l 

(6) 

The probability of this occurrence is 


= | 1 +erf (vft(x,)/2) 

( 7 ) 


where 

R(x x ) = phase accumulator output signal-to-noise ratio 
- f( Xi )l2MN 0 B 

NJ2 s two-sided noise spectral density 
B - noise bandwidth 

The probability of the complementary event that 

/(x, ) +/i <0 
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and 


E(x 2 )--K 2 

( 8 ) 

E(x^E(0)-x 2 *K i 

( 9 ) 

\ -p\x v x 2 )= X -e:ic^JW^n 

00 ) 


These relations depend only on the values of X\ and x 2 and the signal-to-noise 
ratio R. The transition matrix can therefore be calculated exactly for any 
specified signal-to-noise ratio. 


Acquisition Behavior of the Viking Tracking Loop 

A computer program was written to calculate by the above method the 
cumulative probability of acquisition for the Viking subcarrier tracking 
loop. The specific values used in the Viking loop analysis are as follows: 


M 

1 phase increment 

K 2 

R 

m 


— 1G cycles 

*■ 1/512 cycles = tt/ 256 rad 

* 6 increments 
= 1 increment 

= 4 dB in the saturated region of the phase 
characteristic 

* 6.5 increments/ time unit 


The linear region of the phase characteristic extends ±1G increments from 
the zero point. The initial condition state vector p 0 was calculated from the 
statistics of the coarse phase estimator by C. Tegnelia (see Figure 4). 
Because of the large size of the transition matrix (for the above numbers the 
matrix would be 262144 X 262144), the calculations were performed only in 
a reduced area of the state plane. The area used for computation was the 
region -175 < xj <75, - 50 < x 2 <50 (25000 X 25000 transition matrix). The 
degree of approximation induced by using this reduced area is small due to 
the phase error rate and initial conditions that were used. The asymptotic 
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value of the probability of acquisition is biased low, however, because the 
reduced area causes the probability of absorption at the boundary to be 
high. This asymptotic value was determined to be high enough to provide a 
good lower bound for the Viking system. 

The results of the calculation are shown in Figure 5, together with 
experimental values taken from simulation runs. The difference between the 
\ curve and the experimental values is probably due to a slight difference in 

the definition of acquisition between the two cases. The analytical results 
show that the probability of failure to acquire without cycle slipping within 
75 bit periods is less than 10 -5 , a figure which could not be determined 
experimentally in a reasonable amount of time. This result indicated that the 
proposed design would be adequate for the Viking mission. 

Further Applications of the Technique 

In the future th . technique is expected to be used to determine the 
behavior of the tracking loop during the bit synchronisation acquisition 
phase, which follows the subcarrier acquisition phase in the command 



' Figure 4. Experimentally determined initial condition state probabilities 
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0 25 50 75 

TIME, number of v'rto bit pe*iods 


Figure 5. Acpuisition probability of loop 


receiver. In this case the state probability vector can be used to determine 
the probability density function of the phase error at the end of bit 
synchronization phase, during which the effective signal-to-noise ratio in the 
loop is reduced. 

The technique may also be used in the study of the acquisition behavior of 
other tracking loop systems. While the Markov chain transition matrix is 
inherently discrete in nature, it may be applicable to continuous systemj 
with suitable quantized approximations, although the matrix dimensions 
may make this approach unfeasible. 

Summary 

A discrete Markov chain technique has been usei to determine the 
acquisition probability as a function of time for a digital tracking loop. The 
technique has been applied to *be Viking command receiver subcarner 
tracking loop with results that agice very closely lo the experimental results, 
and extend these results beyond those that can be achieved experimentally 


i \ 




% 
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in a reasonable time. The technique is applicable to a wide range of tracking 
loop problems. 


Reference 

1. Holmes, J., Performance of the Second Order All-Digital Command 
System Suhcarrier Timing Loop, Interoffice Memo 331-71-117, Mar. 
24, 1971 (JPL internal document). 
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Bias and Spread in Extreme Value Theory 
Measurements of Probability of Error 


J. G. Smith 

Telecommunications Division 


Performance tests of communication systems characterized by low bit rates and 
high data reliability requirements frequently utilize classical Extreme Value Theory 
(EVT) to avoid the excessive test times required by bit error rate (BER) tests. If the 
underlying noise is gaussian or perturbed gaussian, EVT will produce results either 
biased or with excessive spread, if an insufficient number of test samples are used. 
This article examines EVT to explain the cause of the bias and spread, gives 
experimental verification, and points out procedures that minimize these effects. 
Even under these conditions, however, EVT test results are not particularly more 
significant than BER tests. 


Introduction 

Classically, performance testing of a digital communication system 
consists of counting the numbers of wrong outputs of the system in response 
to a sequence of known inputs. This so-called bit error rate (BER) test 
procedure is known (Reference 1) to produce an estimator of the probability 
of error which converges to the right answer eventually regardless of the 
type of underlying noise; i.e., it is unbiased, consistent, and distribution-free. 
But, the lower the probability of error of the system, the longer the 
sequence length must be to generate a sufficient number of errors to 
construct meaningful confidence intervals; and the lower the data rate of the 
system, the longer the time required to generate this sufficient sequence. 

As a time-saving alternate to BER testing of communication systems 
characterized by low bit rate and low probability of error, classical Extreme 
Value Theory (EVT) (References 1-9) has been utilized. Procedures based 
on classical EVT do not record the number of errors generated by the 
system, but rather examine the noisy data input to the decision element of 
the system, and attempt to measure probability of error from the statistics of 
the noisiest data, i.e., the extremes of the data. Classical EVT is known 
(Reference 2) to be asymptotically unbiased, consistent, and distribution-free 
for a broad class of noise. Unfortunately, the rate of convergence of the EVT 
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estimate to its asymptote is sensitive to the statistics of the underlying 
noise— and particularly slow when that noise is gaussian (References 6-9). 
Further, for a large class of noise, including gaussian noise, the estimator 
suffers either excessive bias or excessive spread for finite sequence length. 
These two features combine to render classical EVT of questionable 
advantage in a gaussian-noise type of environment for these communication 
systems. 

This article documents the bias and spread of EVT in the given test 
environment and indicates how to minimize these effects. Specifically, for a 
given number of test samples, the effects of the number of extremes and the 
number of samples from which each extreme is drawn are examined. 

Classical EVT 

Suppose we have a set of n independent samples (data points) xj,* • *,x„ of a 
random variable X with unknown distribution function 

F x (x) = ?r[X<x] 

and we want to estimate 

Pr[X>x 0 ]=l-F x (x 0 ) 

for fixed threshold % The largest (the extreme) value, x, of this sample has a 
distribution function 

*(*) = Pr[allX.<x] 

= n Pr[Jf.<jc] 

1=1 

= (F x T(x) 

If the random variable X can take on all values of x (that is, it is not limited 
to some finite interval), and its distribution function F x (x ) is asymptotically 
lower bounded by the exponential distribution, then for large enough n, 
<t>„(x) has been shown to be of the asymptotic form (References 1 and 2) 

$>„(*) = exp j -exp [-a n (x -u„)]j 

where a n and u n are unknown parameters. Thus, the procedure followed in 
classical EVT is this: We take several ( k ) independent groups of data, each 
group of size n, and from each group we record the largest (the extreme) 
sample value. Then, using this block of extremes, we estimate a n and u n 

59 




[using, for example, a maximum-likelihood criterion, a minimum distance 
criterion, or possibly others (Reference 2)]. Then our estimate of probability 
of error becomes 


P ~ 1 ~ F x (* 0 ) 

= 1 -($ n ) 1/w (x 0 ) 


= 1 - 


exp -exp l-cx n (x 0 -u n )]\ 


i In 



Note that the total number (m) of samples required for the test is m = n • k. 
Unfortunately, if n is too small, the estimate turns out to be biased. We 
consider why for gaussian noise. 

The Bias of EVT 

We define a new random variable v n (the so-called reduced variate) as a 
function of the underlying random variable X by 


1 

u„(*) = -log n(l -F x (x)^j 


This can be rewritten to obtain 


F x (x) = 1 - exp^-t> n (.x)j In 

4>„M = (F x f(x) 

= |l - exp (-e„(*))/n\” 

t exp (-exp [-e„(x)]) 


Thus, the extremes of X have a distribution function <J> n (x), which asymptoti- 
cally in n has the double-negative exponential form characteristic of classical 
EVT. Convergence to this asymptotic form will be fairly rapid, in terms of 
increasing n, whenever v n (x) is relatively insensitive to changes in n. 
Further, the development so far has left F x unspecified (within the broad 
class mentioned earlier). 

However, the dist-ibution function is in terms of the reduced variate 
values v n (x), rather than the extreme values x themselves. Measured or 
recorded extreme values x must be converted to values v„(x) before statistics 


' * 


■t 
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can be computed. Because the relationship of v n (x) to x depends on the 
unknown functional nature of the characteristics of <J> n (x) need to be 
further understood. 


For example, if F x were exponential, i.e., 

F x (x) - 1 - exp (-ax) 

then we could immediately solve for v n (x), 

v n (x)=a(x-u n ) 


where the only n-dependence appears in the parameter 


= Oog n)/a 


In such a case there would be a nice linear relationship between the 
extremes x and the values v n (x) of the reduced variate, so that 


$„(*) = exp|-exp [ -a(x - « n )] j 


would hold for any reasonably large value of n. 

But if the noise were gaussian (it suffices here to consider a zero-mean, 
unit-variance process), then 


1 - F x (x ) = -j= jf exp (-x 2 12) dx 


which can be expanded asymptotically as (Reference 1) 




where R(x) is positive, small, and falls off as 1/r 2 . Defining 
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e = ~log 1 1 -/?(*)) 
P n = \/2log7z 

log\/27r(p n ) 2 


and manipulating, we obtain a nonlinear relationship between x and v n (x): 

v n (x ) = + ( x ~P n ) 2 l2 + log C*c/P„) + + e 

However, since log y ® y - 1 for t/ * 1, then for values of x near p n , i.e., 
when | x - p n | << p n , we can neglect € and the quadratic term and obtain 

V n (x) * O'/* - u n ) 

where 

cr ^ p + 1/p 
n '*n 

w - C /or 

n n' u n 

Thus, for all x near p n , we have 

$„(*) exp ^-exp [-o^(x - w„)]) 


Classical EVT thus draws a “best” linear fit to the log-log empirical 
distribution function. If we assume we have a sufficient number of extremes 
so that the data “clusters” well around the empirical distribution function, 
and at least a moderate sample size per extreme, then classical EVT also 
draws a “best” linear fit through v n (x). The classical EVT estimate of 
probability of error is a function only of the ordinate of the linear curve at 
threshold (i.e., at Xq). Thus, for classical EVT to yield accurate results, we 
require a large enough n that v n (x ) is reasonably linear near threshold; i.e., 
we require n such that 


l*-P. l«o. 
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for all x near threshold x 0 , or, in particular, we might require 


1*0 


«p 


n 


For example, if x 0 = 4 (i.e., Xq = 4<r with cr = 1), and we require x 0 = p n , 
then 


n = exp [(* 0 ) 2 /2] « 3000 


Thus, we see the need for some minimal sample size n per extreme to 
guarantee dominance of the linear portion of v n (x). For gaussian noise of 
mean m and variance cr 2 , a n is scaled by 1/cr, u n is scaled by cr and shifted 
up by m, and the criterion for dominance of the linear term near threshold 
becomes 



a 



«P 


n 


Furthermore, since d 2 v n (x)/ dx 2 is positive for all x greater than 1, classical 
EVT places a linear curve through an upward-curved function [assuming we 
have enough extremes that the data clusters around the curve o n (x)]. We 
should point out that if the threshold x 0 exceeds u n , then the linear curve 
must cross threshold at a lower value than v n (x), because most of the 
recorded extremes used in making the best linear fit lie near u n , and 
probably none exceed % If we have a large number of extremes exceeding 
threshold, our value of n is probably sufficiently large to permit classical bit 
error testing. Thus, we can assume u n < x 0 . 

Recall now that the true probability of enor p is 


p - 1 - exp 


(-exp [-»„(*„)]) 


l/n 


while the classical EVT probability of error p ’ is 

p' = 1 - (exp j -exp [-o„(* 0 -«„)] J J 


Thus, since 
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then p' y p\ that is, classical EVT estimates an excess probability of error, 
and thus the bias b k p - p' is negative. Clearly, the magnitude of the bias 
will decrease as n is increased, since v n (x) becomes increasingly linear as n is 
increased. 


The Spread of EVT 

Given a set of independent data 1 xj,* • *,x*) , ordered so that 

Y C X ^ M ^ Y 
**1 2 * 

with some distribution function 4> n (x), we note that the fi action of points 
with value no greater than x, (for any i - 1,* • -,k) is i/k. Thus, we could define 
an empirical distribution function e* k (x) by 



for then, in the limit, as k gets large, e* k (x) would converge to <I> n (x), as 
demanded by the law of large numbers. 

To avoid problems with log-log functions of the e* k (x) defined above, 
however, it is conventional to modify the definition somewhat to 

TTk^ (x - */> 

1=1 

Again, there is ultimate convergence of e^x) to <I> n (x). 

We may then define the reduced empirical variate u^x) by 

Hog e nk^l 


so that if 


*,W n CX P ( _eX P ) 


(which it does for moderate n), tf m u n ^x) also converges to c n (x). This 
provides a qualitative measure of the number of extremes needed. We 
require a sufficient number that the reduced empirical variate u^x) clusters 
well around the reduced variate o n (x). If k is too low, then the “best” linear 
fit of the data will not necessarily “fit” well, no matter how linear v n (x) may 
be. That is, an insufficient number of extremes will cause excessive spread of 
the parameter estimates in any sequence of tests. 
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Figure 1 shows the reduced variate v n (x), the reduced empirical variate 
u n k(x), and the linear curve a n (x - u n ) for a sample of 30 extremes, each 
drawn from a block size of 100, with threshold x 0 = 4. 

Experimental Evidence 

We have pointed out that the number k of extremes used must be large 
enough so that the reduced empirical variate clusters near the reduced 
variate, and that the sample size n from which each extreme is drawn must 
be large enough to minimize the bias (effect of nonlinearity). For a given 
total test size m, we must ask what apportioning values of n and k give a best 
tradeoff between bias (n too low) and spread (k too low). Figures 2 and 3 
compare the effects of various values of n and k as a function of test size m. 
These results were obtained from examination of independent samples from 
a gaussian source with known statistics, and a threshold x 0 of 4 sigma points. 
Figure 3 summarizes the results at m — 3000, 12000, 30000, and 60000, for k 
= 10, 20, 30, 50, and 100. As k increases, the spread decreases but moves 
upward, reflecting the increased bias due to decreasing n. 

As a reasonable compromise between the effects of bias and spread, i k of 
30 (30 extremes/ total test) was selected for this threshold of 4o\ and EVT 
was compared with BER using the same gaussian source. The performance 
curves are shown in Figure 4. The five solid EVT performance curves 
include the 90% quantiles, the first sigma points, and the average. The two 
dashed BER performance curves are the 90% confidence curve and the 
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Figure 4. Comparison of BET and EVT for gaussian noise 


actual bit error count. Clearly, EVT holds no distinct edge over BER for this 
data source. Comparable tests at other threshold values also bear out this 
conclusion. 

For noise with density function of the form 

f(x) = a • exp 

for r >1, analysis again reveals a nonlinear form of v n (x) and a negative bias. 
Experiments have shown the same tradeoffs between spread and bias, and 
comparable performance. 


Summary 

We have shown that bias and spread impose conflicting requirements in 
the use of EVT for testing digital communication systems, have indicated 
how to minimize the effects of these factors, and have compared EVT with 
BER testing of a gaussian source. It seems apparent that under the best of 
conditions EVT does not perform significantly better than BER. 
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On the Selection of an Optimum Design Point 
for Phase-Coherent Receivers Employing Band- 

Pass Limiters 

M. K. Simon 

Telecommunications Division 


In the design of phase-coherent receivers employing band-pass limiters, it is 
customary to specify system performance relative to its value at a fixed design point. 
For a given design point, an optimum tradeoff can be found between the power 
allocated to the carrier and the sideband signals. This article describes an attempt to 
further improve the performance of such coherent carrier systems by optimizing the 
design point based upon a given practical optimizing ci iterion. The single-channel 
system is treated in detail and a brief discussion is given on how to extend the 
technique to a two-channel system. 


introduction 

Over the years a great deal of interest has been demonstrated relative to 
the problem of optimum power allocation in single-channel command, and 
one- and two-channel telemetry coherent carrier systems for spacecraft 
applications (References 1-4). One fact in common to all of the solutions 
that have been suggested is that the system design point [i.e., the carrier 
tracking loop signal-to-noise ratio (SNR) at “threshold”] has been arbitrarily 
fixed while the other system parameters are varied to achieve the optimum 
tradeoffs. To date, very little consideration has been given to the question of 
what is the optimum design point based upon a given practical optimization 
criterion. 

To answer this question, the problem must be posed in such a way that 
the relation between an actual operating point and the design point is 
clearly placed in evidence, thereby resulting in a solution of significant 
importance to the practicing engineer. This article suggests a method for 
doing so. The technique involved makes use of the up-to-date phase-locked 
loop and band-pass limiter theories recently contributed in References 5 and 
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Statement of the Problem 


There are two parts to the optimum powe* allocation problem associated 
with the design of a phase-coherent receiver employing a band-pass limiter. 
The first part (which is that most commonly considered in the literature) is 
to fix either the system error probability or total transmitter energy-to-noise 
ratio and minimize the other quantity by varying the system modulation 
indices. This approach can be applied to either single- or two-channel 
systems and implicitly assumes that the system design point is held fixed 
during the optimization. 

The second part of the problem considers, in effect, the locus of these 
minima as the design point is varied, and this information is used to select an 
optimum design point. This part of the optimum power allocation problem 
has, to the author’s knowledge, not been considered and is the principal 
motivation behind this article. 

System Model 

A functional diagram of a phase-coherent receiver preceded by a band- 
pass limiter (BPL) is illustrated in Figure 1. The input signal (point 1) is 
characterized by a carrier biphase modulated by a unit power digital data 
modulation. The loop phase detector, loop filter, and voltage-controlled 
oscillator (VCO) comprise a phase-locked loop (PLL) whose purpose it is to 
establish from the amplified, band-pass limited input signal a coherent 
reference for demodulating the data. In computing the error probability 
performance of the data detector, several possibilities exist depending on 
how the data signal is demodulated off the carrier. Specifically, the following 
three cases are of interest. 

(1) The input IF filter passes the total modulation which is subsequently 
transferred to the data detector directly from the PLL phase detector 



Figure 1. A phase-coherent receiver employing a band-pass limiter 
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output (point 4 in Figure 1). This situation is typical of present and 
past command transponders used in deep-space applications. 

(2) The input IF filter passes the total modulation; however, the 
modulated carrier is tapped off prior to entering the band-pass filter 
and is demodulated by the PLL reference in a separate phase 
detector (point 5). The output of this phase detector represents the 
modulation to be recovered by the data detector. This case is typical 
of present-day methods for recovery of ranging data and is of interest 
in future command applications. 

(3) The input IF filter completely filters out the modulation, w* h the 
modulation being recovered as in case (2). This situation would 
correspond to the case where the data modulation is placed on a 
subcarrier whose frequency is outside the IF bandwidth. This case is 
typical of phase-coherent reception of telemetry data. 

In cases (1) and (2), the BPL suppression factor, and hence the effective 
loop SNR, is a function of the total power-to-noise ratio in the input 
bandwidth. In case (3), the BPL loss is increased [relative to case (2)] because 
it is now only a function of the carrier component of the total power-to- 
noise ratio in the input bandwidth. Hence, with all other parameters 
unchanged, the effective loop SNR is reduced and the noisy reference loss 
increased. This effect produces a minor degradation in case (3) relative to 
case (2). 


Regardless of which case is applicable, an important parameter to 
consider is the design point loop bandwidth-symbol time product. Ordinar- 
ily, for command applications 5 0 = l/iW^T) is sufficiently small (e.g., < 
0.5) such that the phase error can be assumed to vary rapidly over a symbol 
interval T. In many low-to-medium rate telemetry systems, this assumption 
still remains valid. Hence, to a good approximation, the data detector error 
probability P £ is for all three cases described by 


P E = ^ erfc (P d ) 1/2 cos 0 


( 1 ) 


where 



V 

N oo 


COS 0 * E {cos 0} 


( 2 ) 


The parameters S d and N od are the data signal and noise powers, respec- 
tively; both depend upon which of the three cases is being considered. The 
validity of the gaussian noise model is justified as follows. In cases (2) and (3), 
the noise affecting the data detection process is clearly gaussian because it 
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has not been transmitted through the BPL. In case (1), the noise entering the 
matched filter detector (i.e., the phase detector output at point 4 in Figure 
1) is not gaussian, but, for small input signal-to-noise ratio, is approximately 
zero mean. However, because the effective bandwidth of this noise is much 
wider than l/T, the integrate and dump action of the matched filter acts as a 
narrow-band filter and the output statistic is approximately gaussian. 

Optimization of Performance as a Function of Design 

Point 

In designing phase-coherent receivers preceded by band-pass limiters, it is 
convenient to characterize loop performance at a given operating point 
relative to its value at a fixed design point. In the past, it has been customary 
to choose this design point to correspond to a carrier SNR in the design 
point loop bandwidth Y 0 = 2P c0 /(NqW 1j0 ) of 3 dB, i.e., so-called “loop 
threshold.” More recently, it has become apparent that the threshold effect 
exhibited by the loop is a more complicated phenomenon than can be 
expressed entirely in terms of a given carrier SNR. Hence, the original 
motivation behind choosing Y 0 = 2 is, strictly speaking, no ionger 'alid. 
Rather, one should allow Y 0 to be chosen as that value that results in 
optimum receiver performance. Our first step, then, is to divide the power 
between carrier and sideband to minimize the total energy-to-noise ratio, 
PT/N 0 , at the transmitter for a fixed error probability, P E , in the data channel 
and a fixed ratio of data rate, l/T, to design point loop bandwidth, W uq. This 
procedure enables the design engineer to constrain, a priori, the minimum 
acceptable signal-to-noise ratio in the design point bandwidth, x min ; this 
minimum is based upon tracking performance considerations and, once 
established, allows the engineer to seek the optimum power split subject to 
that constraint. If the minimum total transmitted energy-to-noise ratio 
occurs at a value of x — 2P c /NqW w less than x min , the system designer 
selects the value of PT/Nq at x min . As we shall see shortly, this additional 
constraint determines an optimum design point. 

For a given P E , a fixed design point Y 0 , and specified values of the ratio of 
loop filter time constants, t 2 /tj, loop damping factor at the design point, r 0 . 
input IF filter bandwidth, W,, and 8 0 , one can find an infinite set of paired 
coordinate values (x, Rj). For each pair of values, the ratio of the total 
energy per bit to noise density may be found from 1 


K R « 


28. 


( 3 ) 


For any of the three cases being considered, a plot of PT/Nq versus x exhibits 
a minimum. As an example, consider a command application wherein the 


1 For case (1), flj = whereas for cases (2) and (3), R'j = Rj. The parameter Tp is the 

limiter performance factor defined in Reference 6. Also, any symbol or subcarrier 
synchronization and demodulation losses that might be present have been neglected. 
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command data is obtained by demodulating the BPL input with the PLL 
reference signal [case (2)]. Plotted in Figure 2 is PT/N 0 in dB versus x in dB 
for P F = 10 -5 and parameter values r 2 /r j = 0.002, r 0 = 2, Y 0 = 2, = 9 

kHz, W lj0 — 18 Hz, i/T = 4 bps. The minimum value of PT/N 0 is observed 
to occur at x = 5.8 dB, which is greater than Y 0 . If Y 0 is now increased in 
value (all other parameters held fixed), one finds that the minimum PT/N 0 
together with the value of x at which it occurs both continue to decrease. 
Thus, it appears, at least at first glance, that continued improvement in 
performance (in the sense of minimum PT/N 0 ) can be had simply by raising 
the design point indefinitely. Of course, this cannot be true in practice, and 
hence some additional constraint must be placed on the problem to 
counteract this anomaly. 

One practical consideration is that the value of x chosen for the final 
system design must be greater than some minimum value x min . Where x min is 
related to the design point by x min = K 0 Y 0 , and where Y 0 increases, a poini 
is eventually reached at which the value of x at the minimum of PT/N 0 
becomes equal to x min . From that point on, the value of PT/N 0 at x — x min is 
selected as the best operating condition. Figure 3 illustrates a plot of (P17 
N 0 ) min versus Y 0 for case (2) and the same parameters as those used in 
arriving at Figure 2. 2 Also assumed in Figure 3 is that x min is chosen at the 
design point itself (i.e., K 0 — 1). The physical significance of such a choice is 
that the operating point loop bandwidth can never fall below its value at the 
design point, i.e., W^. One notes from Figure 3 that an improvement of 
0.41 dB can be obtained by allowing the design point, Y 0 , to be increased 
from its previously chosen value of 3 dB (i.e., 0 dB in 2 B[jq) to 5.8 dB, the 
point at which the minimum value of ( PT/N 0 ) min is achieved. 


2 The results for case (1), where the data signal is derived directly from th. phase detector 
output, are negligibly different from those of case (2). The difference is strictly due to the 
effect of the BPL loss on flj because the noisy reference loss (approximately 0.2 dB) is the 
same for both. Hence, Figure 3 can also be thought to apply to case (1). 



x, dB 

Figure 2. Total power-to-noise ratio versus carrier signal-to-noise ratio in the design 
point loop bandwidth for a fixed design point: parameters are typical of command 
application 
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Figure 3. Minimum total power-to-noise ratio as a function of design point: 
parameters are typical of command application— case (2) 


Assuming case (3), i.e., the single-data channel is on a subcarrier that is 
outside the bandwidth of the input IF amplifier, Figures 4 and 5 illustrate 
the choice of an optimum design point for two sets of parameters typical of 
telemetry systems for P E = 10~ 3 . In both Figures 4 and 5, t 2 /tj = 0.002, r 0 
= 2, W, = 4 kHz, 1/T = 816 bps; in Figure 4, = 12 Hz; in Figure 5, 

W w = 48 Hz. 

Indicated by dashed lines on Figures 3, 4, and 5 is the asymptotic behavior 
of Equation 3 as x becomes large, corresponding to the case where the noisy 
reference due to the RF carrier becomes negligible. This asymptote satisfies 
the equation 



erfc’ 1 (2/> £ )] 


* *o T , 
28. 


( 4 ) 
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Application to Two-Channel Systems 

This article is concluded with a brief discussion on how the optimization 
technique presented here might be applied in a two-channel system. Of 
particular interest in deep-space applications is the situation where both 
channels are used for the transmission of telemetry information, e.g., the 
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Figure 4. Minimum total power-to-noise ratio as a function of design point: 
parameters are typical of telemetry application— case (1) 
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Figure 5. Minimum total power-to-noise ratio as a function of design point: 
parameters are typical of telemetry application— case (2) 
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science and engineering channels in the Mariner Mars 1969 mission. 
Assuming that the two-channel telemetry model is classified as case (3), the 
optimization equation analogous to Equation 3 is given by 
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where the subscripts 1 and 2 now refer to channels 1 and 2, respectively. For 
given error probabilities in the two channels, i.e., P Ei and P E2 , and fixed 
parameter values for T 2 /T b r 0 , T\, T 2 , W„ W lJ0 , H 0 , and Y 0 , a plot of PT\/N 0 
versus x once again exhibits a minimum. If, as before, Y 0 alone is now 
allowed to vary subject to the constraint x > x min , then, as in the single- 
channel case, there exists an optimum design point in the sense of minimum 
PTi/Nq. The practical significance of this optimum design point depends 
upon the specific parameter values that characterize the system. The 
asymptotic behavior of Equation 5 for large x satisfies an equation similar to 
Equation 4, namely. 
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Jupiter Trapped Radiation Belts 

A. J. Beck 

Project Engineering Division 


Previously developed models for trapped electron radiation belts and trapped 
proton radiation belts in the Jovian magnetosphere are described. The spatial 
distribution of flux and the L-shell dependence of the characteristic energy are 
displayed for both models. Based on these models, the fluence accumulated by a 
Jupiter flyby spacecraft is given in terms of the equivalent 3-MeV fluence for 
electrons and the equivalent 20-MeV proton fluence for protons. Finally, some 
impacts of these fluences on outer-planet missions are described. 


Introduction 

The observed synchrotron radiation from Jupiter implies the existence of 
relativistic electrons trapped in the Jupiter magnetosphere (Reference 1). 
Although no direct observational data exists to confirm the existence of a 
trapped proton belt, both trapping theory and analogy with the Earth’s 
radiation belts suggest that Jupiter should have an energetic trapped proton 
belt as well. Based on the existing models for the trapped radiation 
(References 1, 2, and 3), the number density flux and energy of the particles 
are expected to be quite large compared to the Earth’s radiation belt. 
Consequently, the radiation belts may have a strong impact on the design of 
spacecraft which fly by or orbit Jupiter, and on the missions which are 
designed to investigate Jupiter or investigate the outer planets using Jupiter 
for gravity assist. 

On July 13-15, 1971, a Jupiter Radiation Belt Workshop was conducted at 
JPL. The workshop brought together scientists and engineers from JPL, 
NASA, various universities, and industry with the goal of establishing models 
from which particle flux and fluence distributions of both electrons and 
protons could be derived in forms suitable for use as space vehicle design 
criteria. The workshop produced four models of Jupiter trapped radiation 
belt: a nominal electron model, an upper limit electron model, a nominal 
proton model, and an upper limit proton model. These models are referred 
to as the workshop models. Divine (Reference 2) has described models from 
which the early design requirements for the Thermoelectric Outer-Planets 
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Spacecraft were derived. These models are referred to as the preworkshop 
models. 


Trapped Radiation Belt Electron Models 

Electron models are based on the observed synchrotron radiation from 
Jupiter. Understanding the synchrotron radiation requires a knowledge of 
the magnetic field strength, which is taken as 10 gauss in the equatorial 
plane at a distance of one planetary radius (Rj) from the center of the ■ lanet. 
In addition, in order to structure a spatial distribution in tY? model, radial 
diffusion of solar wind electrons from the magnetosheath and conservation of 
the first invariant (the magnetic moment fi) arc assumed. As a consequence 
of these assumptions, the characteristic energy £ 0 and the differential flux <f> 
of the electrons both vary as L~ 3 , where L is the magnetic shell parameter. 

Assuming a value for the magnetic field strength, the characteristic energy 
of the electrons in the synchrotron region at L = 2 can ue estimated from 
the conservation of the magnetic moment. To make tHs estimate, the value 
of the magnetic moment at L — 2 is set equal to the magnetic moment of 
solar wind particles outside of the magnetosheath. This process leads to a 
value for the electron characteristic energy which is of the same order of 
magnitude as that required for electrons producing the observed synchro- 
tron emission (Reference 1). Although this general agreement provides some 
support for the assumed conservation of the magnetic moment, the electron 
flux (2 X 10 7 electrons/ cm 2 -s) and the characteristic energy (6 MeV) for the 
nominal electron model have been derived at L — 2 from the observed 
synchrotron emission. 

Because of the estimated uncertainties in the magnetic field strength, 
synchrotron beaming, and unfolding of the synchrotron data, the upper limit 
mode! has both the electron flux and the electron energy larger than the 
corresponding nominal model values by a factor of three. 

Figure 1 shows models for the distribution of the electron flux in the 
equatorial plane of Jupiter’s magnetosphere. The dashed lines are prework- 
shop models (Reference 2), and the solid lines are workshop flux models. The 
horizontal axis is the L-shell parameter, while the vertical axis i c the flux of 
omnidirectional electrons having energy greater than zero for the model. 
However, fluxes of electrons having energy much less than the characteristic 
energy may also be present, as for example, in a thermal plasma. All the 
models are constant for L < 2, but for L > 2, the fluxes are proportional to 
L ~ 4 , L~ 3 , E -3 , and L~ 2 for the four lines, starting from the bottom. Although 
this figure does not show the latitude dependence, the flux models are 
assumed to be latitude dependent, having an e-folding value of about 30 deg 
(Reference 3). 

The variation of the characteristic energy is shown in Figure 2 for the 
electron models. These electrons have been assumed to have a differential 
energy spectrum, given by 
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where Eq is the local characteristic energy. The hoiizontal axis represents 
the magnetic shell parameter, as evaluated either in or away from the 
magnetic equatorial plane, and the vertical axis represents the local 
characteristic kinetic energy of the electrons. The dashed lines represent the 
preworkshop upper limit model, and the solid lines represent the workshop 
models. All models’ energies are flat for L < 2, but for L > 2, the energies 
are proportional to L~ 3 , L -3 , and L~ l for the three lines, starting from the 
bottom. 


Trapped Radiation Belt Proton Model 

The workshop proton models are theoretical in nature. Furthermore, no 
observational data exists to confirm the existence of a Jupiter trapped proton 
radiation belt nor to limit the concentration of protons below the limit 
imposed by the maximum concentration that the magnetic field can contain. 

In the absence of any observational data, the same physical processes 
which trap and accelerate electrons are assumed in formulating the nominal 
proton model. Thus, the nominal proton model assumes many of the features 
of the nominal electron model. It is assumed that the equatorial magnetic 
field strength at 1 Rj is 10 gauss, that radial diffusion transports protons from 
the magnetosheath, that the first adiabatic invariant is conserved, and that 
there are no losses, except those which occur at small L values. The number 
density of protons at L = 2 is set equal to the corresponding electron 
number density, and the proton energy at L = 2 is set equal to 10 times the 
corresponding electron energy. It is assumed that the protons are able to 
diffuse past the satellites of Jupiter without loss. 

An upper limit model was also proposed. The proton flux is limited at 
intermediate L values (L = 10) by ion-cyclotron wave instabilities. In 
addition, a 10-gauss equatorial magnetic field, radial diffusion, conservation 
of the first invariant, and losses only at small L values are assumed. The 
assumption of the conservation of the first invariant again leads to the 
characteristic energy varying as L -3 . 

Figure 3 shows models for the distribution of the proton flux in the 
equatorial plane of Jupiter’s magnetosphere. The dashed lines are prework- 
shop models (Reference 2), while the solid lines are workshop proton flux 
models. The horizontal axis is the L-shell parameter, while the vertical axis is 
the flux of omnidirectional protons having energy greater than zero for the 
model. However, fluxes of protons having energy much less than the 
characteristic energy may also be present, as, for example, in a thermal 
plasma. The preworkshop flux models (Reference 2) are flat out to L — 2 
(nominal) and to L - 12.6 (upper limit); beyond these limits they are 
proportional to L~ 55 and L ~ 6 , respectively. The workshop fluxer are 
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Figure 1. Jupiter electron models in the equatorial plane 


proportional to L -3 for all L> 1. The nearly three orders of magnitude 
difference between the workshop models reflects the uncertainty in this 
environmental description. Although this figure does not show the latitude 
dependence, the flux models are assumed to be latitude dependent, having 
an e-folding value of about 30 deg (Reference 3). 

The variation of the characteristic energy £ 0 with L is shown in Figure 4 
for the proton models. The protons have been assumed to have a differential 
energy spectrum, given by 


d<t> 

dE 



[EI(E 0 ) 2 ] e- E ' E o 


where E$ is the local characteristic energy. The horizontal axis represents 
the magnetic shell parameter, as evaluated either in or away from the 
magnetic equatorial plane, and the vertical axis represents the local 
characteristic kinetic energy of the protons. The dashed lines represent the 
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Figure 2. Jupiter trapped electron characteristic energy 


preworkshop upper limit model; note that the upper limit preworkshop 
model energy is independent of L. The solid lines represent the workshop 
models. All diagonal lines represent energies proportional to L 3 , as required 
by any theory (e.g., L-shell diffusion) in which the particle magnetic moment 
is conserved during acceleration. 


Electron and Proton Fluences for Jupiter Flyby Missions 

The electron fluences and the proton fluences have been computed, using 
Jupiter radiation workshop models, for Jupiter flyby missions. In presenting 
the results, the fluence is given in terms of the equivalent 3-MeV electron 
fluence for electrons and the equivalent 20-MeV fluence for protons. The 
equivalent fluence $ c (Reference 4) is defined by 
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Figure 3. Jupiter proton models in the equatorial plane 


where <I> C is the equivalent fluence of particles at the reference energy E n 
d$(E)/dE is the differential energy fluence spectra, and D(E) is the function 
which gives the radiation damage produced by particles of energy E. 

The electron fluence accumulated during flyby missions past Jupiter is 
shown in Figure 5. The horizontal axis represents the altitude of closest 
approach of an equatorial flyby spacecraft, and the vertical axis represents 
the fluence (particles per unit area) of relativistic electrons, equivalent to 3 
MeV according to the displacement damage they produce in silicon 
semiconductors. The dashed lines are based on the upper limit workshop 
electron model, for zero shielding (uppermost curve) and for 1.0-g/cm 2 
shielding (e.g., 150 mils aluminum). The solid curves are based on the 
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Figure 4. Jupiter trapped proton characteristic energy 


nominal workshop electron model for the same two cases. The variation of 
the flyby turning angle (between 60 and 140 deg) and inclination (within ±10 
deg) yields reductions in the fluence less than a factor of two. A considerable 
range of electron fluence exists for each of the outer-planet flyby missions 
corresponding to the possible missions between the minimum and maximum 
closest approach distances. The minimum and maximum altitudes of closest 
approach for the Jupiter-Saturn-Pluto 1977 (JSP 77) and Jupiter-Uranus- 
Neptune 1979 (JUN 79) missions are shown by vertical lines. For these 
missions the fluence is less than 2 X 10 11 electrons/cm 2 . 

The proton fluence accumulated during the flyby of Jupiter is a sensitive 
function of both the closest approach distance and the amount of shielding 
available. Figure 6 shows the proton fluence for a range of closest approach 
altitudes. The horizontal axis represents the altitude of closest approach of 
an equatorial flyby spacecraft, and the vertical axis represents the fluence 
(particles per unit area) of energetic protons, equivalent to 20 MeV 
according to the displacement damage they produce in silicon semiconduc- 
tors. The dashed lines are based on the upper limit workshop proton model, 
for zero shielding (uppermost curve) and for 1.0-g/cm 2 shielding (e.g., 150 
mils aluminum). The solid lines are based on the nominal workshop proton 
model for the same two cases. The variation of the flyby turning angle 
(between 60 and 140 deg) and inclination (within ±10 deg) yields reductions 
in the fluence less than a factor of two. The radiation effects data for some 
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Figure 5. Equatorial flyby electron fluences (shielded and unshielded) 


piece-parts show that the threshold for damage is a fluence approximately 3 
X 10 9 protons/cm 2 (Reference 4). By selecting radiation-resistant parts and 
employing design techniques to reduce the sensitivity of electronic 
assemblies, it may be possible to increase the threshold fluence by two to 
three orders of magnitude. 


Conclusions 

Models for the electron radiation belt and the postulated proton radiation 
belt at Jupiter have been described. These models resulted from the Jupiter 
Radiation Belt Workshop that was held at JPL during July 13-15, 1971. 
Proceedings for the workshop are expected to be published in the near 
future. In the interim, the models described here represent what the group 
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Figure 6. Equatorial flyby proton fluences (shielded and unshielded) 


considered to be the best available description of the Jupiter radiation belts 
available at that time for the purpose of determining space vehicle radiation 
design requirements. The group, however, recognized the need to examine 
specific features of the models to insure internal consistency with the basic 
assumptions. For example, the effects of relativistic corrections on the 
assumed constancy of the first adiabatic invariant should be carefully 
examined. In addition, the spatial distribution for the flux models will 
require modification. Consequently, the models of Jupiter’s radiation belts 
are still evolving. Yet, any revisions which may be required are not expected 
to affect the basic conclusions that are drawn here. 
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The Jupiter radiation belts, as described by these models, generally 
present a greater hazard to spacecraft than the Earth’s radiation belts do. 
The most important reasons are that the Jupiter belts are much more 
extensive and that the particle energy at L values between 1 and 2 are 
generally larger than the energy of particles in the Earth’s radiation belts. In 
addition, the proton flux is larger at the small L values. However, a 
considerable uncertainty exists in the electron flux as well as in the proton 
flux. A measure of this uncertainty for the proton flux can be readily seen by 
the nearly three-orders-of-magnitude difference between the nominal and 
upper limit flux models in Figure 3. 

The electron and proton fluences are shown in Figures 5 and 6 for flyby 
missions in the equatorial plane having various closest approach altitudes. 
Comparing the fluences received for closest approach distances near the 
planet with the fluences that damage electronic parts (as given, for example, 
in Reference 4) shows that the proton fluence is generally more damaging 
than the electron fluence. Consequently, from a spacecraft design view- 
point, the proton model is more significant. 

Figures 5 and 6 also show the relative effectiveness of radiation shielding. 
Radiation shielding clearly becomes more effective as the altitude of closest 
approach increases. This increased effectiveness results from the softening of 
the radiation because of the L -3 spatial dependence of the characteristic 
energy. A spacecraft probab'y has a minimum inherent self-shielding of the 
electronic parts, ranging from 0.66 g/cm 2 (100 mils of aluminum) to 1 g/cm 2 
(150 mils of aluminum). When this inherent self-shielding is considered, the 
fluence reaching the electronic parts is so small for closest approach 
altitudes greater than 4 Rj for the nominal proton model and about 6 Rj for 
the upper limit proton model that the protons are no longer a significant 
design problem. However, there could still be some effects on science 
instruments. 

For the electrons, radiation effects data (Reference 4) generally indicate, 
with the possible exclusion *f some part types, that parts can be selected 
which can tolerate electron fluences somewhat in excess of 10 11 electrons/ 
cm 2 . Consequently, considering the spacecraft self-shielding and by selecting 
and screening parts, the spacecraft can be made insensitive to the radiation 
fluences described by the nominal models for closest approach distances 
well inside of 4 Rj. The corresponding distance for the upper limit models is 
well within 6 Rj. 
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The flow field of gases and scattered particulates in the forward direction from a 
rocket nozzle is not presently predictable on a quantitative basis. Qualitative tests 
have been made with water injected through the nozzle wall and the droplet 
trajectories observed photographically. These tests were conducted for nozzles 
where the boundary layer is a significant portion of the flow. The tests indicated that 
both gases and particulate matter will be found in the region outside of the plume 
boundaries that can be calculated using current analytical techniques. A test 
program is being conducted in the JPL Molsink facility to measure the distribution 
of exhaust gases from small nozzles using an electron beam/photomultiplier system 
and a matrix of quartz crystal microbalances. Analysis of data from these tests, using 
nitrogen, will result in a better understanding of the backscattering flow field. 
Calibration tests of the Molsink, using simulated hydrazine exhaust product gas 
mixtures, have also been conducted, and future tests with small hydrazine thrusters 
are being considered. The test program and results of the pumping calibration tests 
with hydrazine (simulated hydrazine exhaust products) are reported. 


Introduction 

The interaction of plumes with spacecraft components can cause 
problems in such diverse areas as flight mechanics, materials and structures, 
communications, and scientific experiments. The possible problems are 
aggravated if a large percentage of the plume is composed of nozzle 
boundary layer gases. Because these have a lower velocity than the nozzle 
core gases, they have a higher potential of expansion and will cause a 
backscattering, or forward flow, of the exhaust gases. Although analytical 
methods, such as the method of characteristics, are available to predict 
plume characteristics for nozzles with negligible boundary layers, no 
equivalent methods are available for those with large boundary layers. 
Because the latter type of nozzles are typical of those used aboard 
unmanned spacecraft today, an experimental and analytical program to 
predict the nature of these plumes has been undertaken. 
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The JPL Molsink is being used to study the plume characteristics of small 
thruster nozzles. To obtain accurate measurements of the plume gas flow 
distribution, a number of special systems have been built and evaluated for 
use in the Molsink. An electron beam /photomultiplier system is used to map 
plume densities, and quartz crystal microbalance systems are used to detect 
the mass flux at a number of locations along the wall of the chamber. In 
addition, a gas supply system was designed, which can be adapted to a 
variety of test nozzles. These systems and some interim results of the test 
program are described in the following sections. To obtain a qualitative 
measure of the dispersal of small liquid droplets from a nozzle, a simple 
water injection test was conducted. Water was injected through a static 
pressure hole near the nozzle exit plane. Figure 1 is a photograph taken at a 
vacuum of 133.3 N/'in 2 (1 torr). 1 The nozzle exit plane is in the center of the 
photo and the flow is downward. It is significant to note that even under the 
conditions of a finite back pressure, the photograph indicates that some of 
the water droplets are “backscattered,” that is, they flow upward, forward of 
the exit plane of the nozzle. 

Facility and Test Apparatus 

The Molsink is an ultra-high-vacuum facility which consists of three 
concentric chambers. The inner-most chamber is a sphere approximately 
3.048 m (10 ft) in diameter, maintained at a temperature between 15 and 20 
K with gaseous helium. The walls are wedge-shaped, resembling an anechcic 
chamber, with a total surface area of approximately 185.8 m 2 (2000 ft 2 ). The 
chamber walls are also coated with titanium, which acts as “getter” material 
to trap helium and hydrogen that are not cryopumped by the 15-K surfaces. 
The titanium coating can be renewed while the chamber is evacuated. 
Recent tests, described in a later section, have shown that by coating the 
walls with carbon dioxide, the amount of hydrogen that can be pumped is 
greatly increased over that which can be pumped by the gettering action of 
the titanium alone. 

The test gases are supplied from gas cylinders, located on a platform 
above the Molsink, into the nozzle plenum chamber, which is a 3.8-cm (1.5- 
in.) diameter tube penetrating through the top of the Molsink. It can be 
raised or lowered to change the location of the nozzle in the Molsink by 
about 15.24 cm (6 in.). The plenum chamber is instrumented to record both 
plenum temperature and pressure immediately above the nozzle throat. 
Heating elements control the temperature of the test gases. 

A variety of small nozzles is available, with area ratios as low as 15:1 and 
as high as 240:1. Wall half-angles of these conical nozzles are 15, 25, and 40 
deg. The throat diameter of most of the nozzles is 0.254 cm (0.100 in.) (see 
Reference 1). Pressure taps are located at the throat and exit plane a« well as 
at the midpoint. 


1 Values in customary units are included in parentheses after values in SI (International System) 

units if the customary units were used in the measurements or calculations. 
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An electron beam /photomultiplier system measures the densities using 
the molecular fluorescence technique. A thin collimated beam of elections 
at 10 to 20 kV is projected into a gas, resulting in an illuminated path, the 
intensity and spectral characteristics of which can be related to gas density. 
The fluorescence is measured with a photomultiplier. Its output, normalized 
bs the electron beam current, is calibrated against known gas densities. For 
nitrogen a pressure range of 6.67 X KH to 1.33 X 10* N/m 2 (5 X 10 -fi to 
1(H torr> c;<n be measured. At the low density limit of operation, there is 
insufficient light intensity, while at the high density limit quenching 
collisions suppress light intensity. In terms of Mach number, the lower limit 
means that, for a nozzle using nitrogen and operating at a nozzle plenum 
pressure of 6.89 X 10 3 N/m 2 ( 1 psia), densities in the region where the Mach 
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number is approximately 18 can be measured. Orientation of the electron 
beam /photomultiplier system is schematically shown in Figure 2. Motors, 
which are not shown in the figure, drive the system: one moves the electron 
beam and photomultiplier assembly relative to the test nozzle, and another 
moves the photomultiplier relative to the electron beam gun so that the 
density along the beam path can be measured. Potentiometers continuously 
record the position of these parts of the system. Mounted along the wall of 
the Molsink chamber are a series of quartz crystal microbalances (QCM), 
which are capable of detecting mass fluxes on the order of a monolayer or 
lens. The key element of the QCM systems is a piezoelectric crystal which 
changes resonant frequency when mass is added to its surface. For a given 
temperature, the relationship between frequency and mass deposition is 
linear. By pairing quartz crystals the temperature effect is eliminated. The 
QCM is capable of detecting a cumulative buildup of mass, and it can detect 
discernible particles impinging and adhering on the crystals. Figure 3 is a 
schematic of the chamber, showing the relative location of the test nozzle, 
the electron beam and photomultiplier, and the quartz crystals mounted 
along the chamber wall. Mounted on the bottom door of the chamber is the 
electron beam apparatus used to sublimate the titanium which is deposited 
on the chamber wall surfaces. More detailed information about the chamber 
can be found in References 2, 3, and 4. 



Figure 2. Electron beam/photomultiplier angles 


92 


. V *• 


20 K CRYOGENIC 
QUARTZ CRYSTAL 
MICROBALANCES 



Figure 3. Molsink 







* 


Hydrogen Pumping Capacity 

Molsink pumping capacity tests were conducted to determine its 
capability for plume tests of hydrazine thrusters in the 0.222- to 0.444-N 
(0.05- to 0.10-lb) range. The exhaust products from these engines contain 
hydrogen, ammonia, and nitrogen. The nitrogen and ammonia are easily 
cryopumped by the 15-K walls of the Molsink. Previous calibration tests, 
using nitrogen, have indicated that the Molsink pressure could be main- 
tained at 1.33 X 10 -3 N/m 2 (10 -5 torr) with flow rates as high as 3 g/s. 
Hydrogen is more difficult to cryopump because 5-K walls are required to 
condense the hydrogen. To maintain the walls at this temperature, liquid 
helium would have to be circulated through the wall, which is expensive. 
Some hydrogen will be absorbed by the titanium, which is on the Molsink 
wall, but calculations indicate that the amount of hydrogen that can be 
pumped by the walls is inadequate. 

Pumping of hydrogen by carbon dioxide cryodeposits has been reported in 
Reference 5. This method consists of injecting carbon dioxide gas into a 
vacuum chamber whose walls have been cooled to about 20 K or less. The 
carbon dioxide forms a porous frost on the walls and removes the hydrogen 
by sorption pumping. Therefore, calibration tests of the Molsink using 
hydrogen/nitrogen mixtures were conducted. Hydrazine exhaust products 
contain hydrogen, nitrogen, and ammonia, but it was decided to substitute 
nitrogen r oi the ammonia to simplify the tests. A mixture of 53% hydrogen 
and 47% nitrogen, by volume, was used. This percentage of hydrogen is 
typical in the exhaust products of small hydrazine monopropellant thrusters. 


93 


V ■**>!**£ 



r 


t 


/‘ 




Carbon dioxide gas was injected through the test nozzle at a flow rate of 
0.2 g/s for 24 min to coat the walls of the Molsink. After the frost coating 
had been allowed to stabilize on the walls, the hydrogen/nitrogen mixture 
was injected through the test nozzle at a flow rate of approximately 0.1 g/s 
as shown in Figure 4. The results of the tests are shown in Figure 5, which 
indicates that a pressure in the 1.33 X 10 _2 -N/m 2 (10 _4 -torr) range was 
maintained for approximately 4Vi min. The pressure in the Molsink prior to 
the test gas injection was 6.67 X 10~ 4 N/m 2 (5 X 10~ fi torr). The pressure 
rose to 2.67 X 10 -2 N/m 2 (2 X 10 -4 torr) within 0.3 min after the test gas was 
injected. This pressure was maintained for 4Vi min, after which the pressure 
rose rapidly (an order of magnitude per minute) and the gas flow into the 
chamber was terminated. The tests have indicated that the use of carbon 
dioxide frost is an inexpensive and convenient method for pumping 
hydrazine exhaust products in the Molsink. The hydrogen/nitrogen gas flow 
rates used during the tests were those equivalent to a 0.222-N (0.05-lbf) 
thruster. Future Molsink calibration tests are planned to determine the 
pumping capacity, using carbon dioxide frost, and test gas flow rates 
equivalent to 0.444-N (0.10-lbf) thrusters. In addition, cyclic test sequences 
will be used to see if by renewing the carbon dioxide frost coating before it 
is saturated the Molsink pressure can be prevented from rising above the 
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Figure 4. Molsink hydrogen/nitrogen gas mixture injection parameters 
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Figure 5. Moisink pressure history of hydrogen/nitrogen mixture injection 


value established when the test gas is initially injected. The amount of 
carbon dioxide injected will also be varied to determine the relationship 
between the frost coating thickness and mass, the chamber pressure, and the 
amount of hydrogen pumped. 

Future Test Plans 


Tests are continuing in the Moisink using the electron beam and quartz 
crystal systems. The data from these tests will be used to help establish an 
analytical model that can be used to predict the plume characteristics of 
small thrusters where the boundary layer is an important influencing factor. 
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New highly stable polymers are required for materials applications on future long- 
term planetary missions. This article presents the results of a systematic investiga- 
tion on the use of anhydrides to prepare stable elastomeric materials using mild 
reaction conditions. The three anhydrides investigated were found to provide 
effectivt chain extension of hydroxy-terminated poly(alkvlene oxides) and poly(buta- 
dienes). These were tetrahydrofuran tetracarboxylic dianhydride, pyromellitic 
dianhydride, and benzophenone tetracarboxylic dianhydride. The most effective 
catalyst investigated was ferric acetylacetonate, which resulted in chain extension at 
333 K (60°C). A novel featur: of these anhydride reactants is that they are 
difunctional as anhydrides, but tetrafunctional if conditions are selected that lead to 
reaction of all carboxyl groups. Therefore, chain extension can be effected and then 
followed by crosslinking via the residual carboxyl groups. 




Introduction 

The objective of this work is the development of stable elastomeric 
systems for solid-propellant binders applicable to long-term planetary 
missions. Achieving this objective requires dimensional and mechanical 
stability as well as insensitivity to radiation. Such elastomers should evolve a 
minimum of volatile products. These new elastomers should be useful in 
other applications. Previous work has shown that anhydrides can be used to 
chain extend hydroxy-terminated prepolymers (Reference 1 and related 
unpublished JPL work); however, chain extension was limited and reaction 
conditions were relatively severe. The work reported here was conducted to 
obtain higher molecular weights and determine the effects of different 
catalysts, diols, and dianhydrides. 

Diols included polypropylene oxide) (PPO), polyethylene oxide) (PEO), f 

and poly(butadienes). For a polymer of increased stability, saturated ; 

hydrocarbon prepolymers will be more extensively tested in future work. 

Anhydrides included tetrahydrofuran tetracarboxylic dianhydride 


97 


.*1 








(THFTDA), pyromellitic dianhydride (PMDA), and extension can be 
illustrated by 


HO(R) m OH 



(PMDA) 



where R is a monomeric unit in the prepolymer, extension can be illustrated 
by where R is a monomeric unit in the prepolymer, m is the prepolymer 
degree of polymerization, and n is the degree of polymerization of the 
extended prepolymer. 


Results and Discussion 

For most of the experiments, poly(propylene oxide) cf moi wt 2000 and 
hydroxyl of 1.0 eq/kg (Union Carbide, Lot 439430) was used as the diol. 
Approximate intrinsic viscosities [rj] were calculated from the expression 

V c= M +a39f c 


from measurements at one concentration C in benzene c ^ylformamide 
at 303 K (30°C). This relationship is strictly valid only for PPO-toiuene 
diisocyanate (TDI) polymers (Reference 2), but is used here for comparison 
of anhydride polymers with the values obtained for PPO-TDI, as an 
approximation of the degree of chain extension. For the extension of 
prepolymers, the dianhydride and low molecular weight polymer were 
stirred at 448-468 K (175-195°C) to dissolve the anhydride; the catalyst (if 
any) was added, and the solution was cured in a vacuum oven at the desired 
temperature. 

Since imidazole had been found to be a catalyst for the reaction of 
alcohols with PMDA for hydroxyl analytical methods (Reference 3), it was 
tested as a catalyst for the reaction of PPO with PMDA. The desired 
catalytic effect was found to be negligible. At a mole ratio of PMDA to PPO 
(AN /OH) of 1.0 and a curing temperature of 373 K (100°C), the values of [17] 
were 1.6 and 1.5 m 3 /kg in the presence (21-h cure) and absence (24-h cure) 
of imidazole, respectively. The catalyst concentration was 0.2 mole per mole 
of PMDA. There are two functions of a catalyst in systems of this type. Not 
only should it increase the rate of the desired chain-extension reaction, but it 
is also desirable that it increase this rate relative to side-reaction rates. 
Therefore, although imidazole was found to increase the chain-extension 
rate at times shorter than 24 h, the final molecular weight for longer cure 
times was no higher when it was present, indicating that side-reaction rates 
were also catalyzed. When short-reaction times are desi»ed, as for analytical 



methods, imidazole is an effective catalyst for anhydride-hydroxyl reaction. 
Further work showed that imidazole, pyrole, and pyridine did not give 
higher values of [ 17 ] than no catalyst for cure times from 14 to 134 h at 333 K 
(60°C) with a mole ratio of PMDA to PPO of 1.05. The intrinsic viscosities 
were always 1.1 (±0.3) kg/m 3 under these conditions. 

However, the presence of ferric acetylacetonate (FeAA) resulted in an 
intrinsic viscosity ~2.4 m 3 /kg in 66 h at 333 K (60°C). The ratio of FeAA/ 
anhydride was 0 . 2 ; thus, the catalyst concentration was higher and the effect 
less than for isocyanate-hydroxyl reactions. In an experiment in which the 
FeAA to anhydride or diisocyanate ratio was 0.002, the extended anhydride 
and isocyanate polymers had [ 17 ] = 1.8 and 8.8 m 3 /kg, respectively. 
However, the reactant ratio was probably more nearly optimum for the 
urethane polymer, accounting for part of the larger increase in molecular 
weight. 

Figure 1 indicates the effect on [ 77 ] when PPO was chain-extended with 
several different anhydrides in the absence of catalyst at 333 K (60°C). 
THFTDA resulted in the highest molecular weights and was as effective as 
TDI under these conditions. Inspection of the results in Figure 1 indicates 
that the rate of chain extension may depend on the solubility parameters of 
anhydride and diol. 

From Figure 2 , it can be seen that the anhydride-hydroxyl mole ratio is in 
the vicinity of 1.05 for maximum chain extension. The value appears to be 
relatively critical and somewhat less than the value of 1.10 generally used for 
NCO/OH ratios. Figure 2 also shows that FeAA is a much more effective 
catalyst for anhydride extension than CuAA or NiAA, which is also true for 
the isocyanate reaction. 

Some experiments were also carried out at higher temperatures, for 
which gelation times were measured. Results are indicated in Figure 3. 
Because of esterification following anhydride ring opening, it appears that a 
crosslinked network can be formed and gelation may occur in 150 h at 373 K 
( 100 °C). The apparent upward curvature at lower temperatures is probably 
real, but, even if it is not, the relative rate of anhydride opening would be 
approximately six times faster than esterification at 333 K (60°C), indicating 
that negligible esterification would take place under normal low-tempera- 
ture curing conditions for this system. 

Table 1 shows some results obtained for chain extension and crosslinking 
of polyethylene oxide), poly(butadienes), and a reduced saturated poly(buta- 
diene). The higher temperature, 463 K (190°C), was used to dissolve the 
anhydride in the prepolymer. For these prepolymers, considerable chain 
extension and some crosslinking took place during the dissolution, as 
indicated by high initial intrinsic viscosities and gel formation within 10 min 
at 463 K in some cases. 
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TEMPERATURE = 333 K (60°C) 
AN/OH = 1 .05 



Figure 1. Solution viscosity of anhydride chain-extended poly (propylene oxide) 


Conclusions and Recommendations 

It has been found that THFTDA readily chain-extends PPO to high 
molecular weights in the presence of FeAA catalyst at 343 K (70°C) or 
lower. The reaction with hydrocarbon prepolymers appears to take place at 
higher rates, even in the absence of a catalyst. 

Further work, particularly with saturated hydrocarbon prepolymers, 
should be conducted to define the minimum dissolution temperature and to 
determine which dianhydride will produce maximum chain extension and 
minimum esterification at moderate temperatures. Crosslinking of the 
carboxyl groups formed on ring opening should be investigated further, 
especially to form crosslinks at temperatures below 373 K (100°C). 
Trifunctional hydroxyl compounds, as well as metallic salts or oxides, should 
be tested. The former would be analogous to triols in urethane system, 
whereas the latter should provide salt crosslinks that would react with some 
of the free carboxyl groups and may tend to increase chemical stability. 
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Figure 3. Gelation time versus temperature for anhydr ide-poly(propylene oxide) 

polymers 
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Table 1. Prepolymer chair, extension with THPTDA 


Prepolymer 

Time, 

min 

Temperature, 

K ( C> 

Mole ratio 
AN/OH 

Initial [tj], a 

m 3 /kg 

Final [^], a 
m 3 /kg 

Polyethylene oxide) h 

3960 
(66 h) 

345 (72) 

1.05 c 

0.8 

1.6 

Telagen-S^ 

10 

463 (190) 

1.05 

C.8 

2.3 

Telagen-S^ 

60 

463 (190) 

1.1 6 C 

— 

Gel 

Telagen-S^ 

1080 
(18 h) 

345 (72) 

1.05 

— 

Gel 

Hycar HTB e 

10 

463 (190) 

1.00 c 

0.6 

6.4 

Butarez HTS f 

10 

463 (190) 

1.05 

1.5 

3.0 

Butarez HTS f 

18 

345 (72) 

1.05 

— 

Gel 

Poly-BD, CN-15? 

10 

463 (190) 

1.05 

— 

Gel 

Poly-BD, R-45M? 

10 

463 (190) 

1.05 

- 

Gel 

Poly-BD, R-15M? 

10 

463 (190) 

1.05 

— 

Gel 

Poly-BD, R-45HT? 

10 

463 (190) 

1.05 

— 

Gel 

Poly-BD, CS-15? 

10 

463 (190) 

1.05 

— 

Gel 


a Viscosities obtained by single-point measurements using 7j S p/C =* [tj] + 0.4 [tj ] 2 C\ 

^ A. G. Fluka, mol wt 2000. 
c Contained FeAA/anhydride mole ratio = 0.0S. 

^ Saturated poiy(butadiene) a,<u-diol, GT and R, Lot 242 AM 148AH, hydroxyl 0.91 eq/kg. 
e Poly(butadiene) a,<w-diol, B. F. Goodrich, Lot V-74, hydroxyl 0.38 eq/kg. 

* Poly(butadiene) a.oi-diol, Phillips Petroleum, Lot 4760, hydroxyl 0.51 eq/kg. 

? Hydroxy-terminated poly(butadienes), ARCO Chemical (all hydr- xyl values from 
supplier's data): 

CN-15, Lot 003061, hydroxyl 0.58 eq/kg 
R-45M, Lot 008281, hydroxyl 0.70 eq/kg 
R-15M, Lot 710291, hydroxyl 0.65 er i kg 
R-45HT, Lot 006041, hydroxyl 0.84 eq/kg 
CS-15, Lot 912211, hydroxyl 0.62 eq/kg 


2. Moacanin, J., “Diisocyanate-Linked Polymers. I. Dilate Solution 
Properties of Toluene Diisocyanate-Extended Polypropylene Glycol,” 
/. Appl. Polym. Sc\, Vol. I, No. 3, pp. 272-282, 1959. 

3. Kingston, H. M., Garey, J. J., and Hellwig, W. B., “Imidazole Catalysis 
of Organic Hydroxyl Group Determination by Reaction with Pyroinel- 
litic Dianhydride in Dimethylformamide,” Anal. Chem ., Vol. 41, pp. 
86-89,1969. 
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Mariner Venus-Mercury 1973 Midcourse Velocity 
Requirements and Delivery Accuracy 

E. L. McKinley 

Mission Analysis Division 


The Mariner Venus-Mercury 1973 Mission represents the first attempt to navigate 
a single spacecraft to more than one planet. The primary mission consists of 
encounters with Venus and Mercury (a second encounter with Mercury is also 
possible). In this study, the expected navigation sequences arc simulated with a 
Monte Carlo computer program for the purpose of determining midcourse 
correction velocity requirements and delivery accuracies. These simulations provide 
the sensitivity of the velocity requirements and delivery accuracies to the error 
sources affecting the navigation process. The orbit determination capability at the 
final pre-Venus maneuver is shown to be the dominant contributor to the velocity 
requirements for the primary mission. Similarly, the orbit determination capability 
at the final pre-Mercury maneuver is shown to be the dominant contributor to the 
delivery accuracy at Mercury. 


Introduction 

In the 1970’s, several opportunities exist for using the gravitational field 
during a planetary encounter to assist a spacecraft in attaining still another 
planetary encounter. (Missions of *his type are made possible by the 
particular geometrical alignments of the planets.) One such opportunity 
occurs in the 1973-1974 time period and involves the planets Earth, Venus, 
and Mercury. Known as the Mariner Venus-Mercury 1973 mission, the 
primary mission consists of encounters with Venus and Mercury (with 
scientific emphasis on Mercury) where the encounter with Venus is designed 
to provide departure conditions necessary to achieve the Mercury encoun- 
ter. (For a given Earth launch date and Mercury arrival date, arrival 
conditions at Venus are fixed.) Although not presently a part of the Mariner 
Venus-Mercury 1973 mission, a second Mercury encounter may be possible. 
This is made possible by a second fortuitous circumstance, namely, that the 
period of the post-Mercury heliocentric trajectory is approximately twice 
the orbital period of Mercury. Therefore, for appropriately selected aim 
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points at the first Mercury encounter (Mercury I), a second encounter 
(Mercury II) will occur approximately 176 days later. 1 

The essential elements of the mission described above are given by Figure 
1, which shows a heliocentric plane view for a typical Mariner Venus- 
Mercury 1973 trajectory. Also shown is an edge view of the ecliptic plane 
illustrating the motion of the spacecraft out of the ecliptic plane. Additional 
detail on the mission design is given by References 1, 2, and 3. 

This article describes results of a Mariner Venus-Mercury 1973 navigation 
study directed toward answering the following questions: 

(1) What midcourse maneuvers are required to accomplish the primary 
mission, and to accomplish an extended mission consisting of a second 
encounter with Mercury? 

(2) What are the corresponding maneuver magnitudes? 

(3) How well can the delivery to the target planets be accomplished? 

(4) What is the relative importance to the midcourse maneuver 
magnitudes of the various error sources which contribute to the 
overall navigation errors? 

The answers to these questions depend upon the actual spacecraft and 
mission designs, and certain of the results described here are subject to 
minor changes as the design continues. However, unless there are considera- 
ble changes in certain areas affecting the navigation process, the characteris- 
tics of the navigation and maneuver strategy given below are a good 
representation of what can be expected during the Mariner Venus-Mercury 
1973 mission. 2 


Maneuvers Required 

The Mariner Venus-Mercury 1973 mission is one of the more ambitious 
interplanetary space flights which has been undertaken. This is particularly 
true of the navigation phase. Table 1 presents a typical Mariner Venus- 
Mercury 1973 mission sequence for which the maneuver magnitudes, 
delivery capability, and error sensitivities have been determined. The 
rationale, described below, for selecting this maneuver sequence is based 
upon experience from previous Mariner missions. 

The mission begins with the injection of the spacecraft by an Atlas/ 
Centaur launch vehicle. The predicted trajectory dispersions resulting from 


1 Aim points selected for Mercury I based on obtaining good return from the scientific 
experiments do not preclude Mercury II; indeed, little compromise is required for certain 
launch/arrival dates (Reference 1). 

2 The Mission Analysis Division at the Jet Propulsion Laboratory is currently studying many o' 
the error sources contributing to the navigation accuracy. Some of this effort is reported in 
Reference 4. 
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V r NUS, 2/5/74 



OUT-OF-ECl IPTIC DISTANCES 
ARE EXAGGERATED 2 TIMES 
ACTUAL SIZE. 


Figure 1. Heliocentric plane view for Mariner Venus-Mercury 1973 trajectory (dots 

are at approximately two-day intervals) 


this Injection are fairly large (0.99 probability of falling within an ellipse in 
the B-plane, Figure 2, of Venus with a semi-major axis slightly in excess of 6 
X 10 4 km), resulting in a virtually certain requirement for at least one 
midcourse correction. In order to minimize the magnitude of the first 
maneuver, it is performed as soon after injection as is practicable, the only 
limitation normally being the time required to obtain a good estimate of the 
post-injection orbit from tracking data. 

The orbit determination dispersion at the time of the first maneuver and 
the maneuver execution errors resulting from the first maneuver typically 
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Table 1. Typical event sequence for the Mariner Venus-Mercury 1973 mission 3 


Event 


Time (nominal) 


Primary mission: 


Launch 

First midcourse (launch +10 days) 

Second midcourse (Venus encounter - 3 days) 

Venus encounter 

Third midcourse (Venus encounter + 3 days) 

Fourth midcourse (Venus encounter + 20 days) 

Mercury I encounter 

End of primar mission (Mercury I encounter + 15 days) 
Extended mission: 

Fifth midcourse (Mercury I encounter + 37 days) 

Sixth midcourse (Mercury I encounter + 149 days) 
Mercury II encounter (Mercury I encounter + 176 days) 


Oct 24, 1973 
Nov 03, 1973 
Feb 02, 1974 
Feb 05, 1974 
Feb 08, 1974 
Feb 24, 1974 
Mar 30, 1974 
Apr 14, 1974 


May 06, 1974 
Aug 26, 1974 
Sep 22, 1974 


a For this study, an Earth launch /Mercury I arrival of Oct 24, 1973/Mar 30, 1974 was 
used. This trajectory is quite representative, and other launch/ arrival combinations 
will not appreciably alter the findings of the study. 
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Figure 2. Encounter R-S-T coordinate system 


combine to give a three-sigma dispersion ellipse of several hundred 
kilometers. On past single-planet Mariner missions, potential errors of this 
size have been acceptable; however, Figure 3 shows the effect on delivery to 
Mercury of a 1-km spherical error at Venus. Consequently, a second 
midcourse correction is virtually certain to be required. 


To minimize the contribution of the maneuver execution errors on the 
delivery accuracy at Venus, the second maneuver is scheduled just a few 
days before Venus Encounter. Consequently, the uncertainty in the orbit 
estimate used to target the maneuver determines the delivery capability to 
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Figure 3. Contour of the B-plane miss at Mercury resulting from a 1-km spherical 
miss at Venus for the Mariner Venus-Mercury 1973 mission 

Venus. Because of the error magnification at Mercury shown in Figure 3, the 
size of the second maneuver orbit determination uncertainty (while small by 
current standards) will make still another maneuver necessary following 
Venus encounter. Indeed, the Venus-to-Mercury navigation problem is 
practically a duplication of the Earth-to-Venus navigation problem. 
Likewise, if the mission is extended to include a second Mercury encounter, 
the above arguments would apply to the Mercury-to-Mercury navigation 
problem with the following exception: because of a thermal constraint, the 
spacecraft is precluded from performing midcourse corrections at distances 
closer to the Sun than 0.58 AU. The timing of maneuvers 4, 5, and 6 satisfies 
this constraint. 


Maneuver Magnitudes 

A Monte Carlo computer simulation of the mission navigation sequence 
described above was used to determine the maneuver magnitudes. The 
program assumes a linear relationship between perturbations in velocity and 
in encounter parameters, the appropriate matrices being obtained by 
integrating the variational forms of the equations of motion about the 
nominal trajectory. Orbit determination characteristics are described by 
covariance matrices of encounter parameters. Other inputs to the program 
include pointing accuracies, propulsion system impulse accuracies, and a 
covariance matrix of Atlas/Centaur injection errors. 

Table 2 gives the results of the Monte Carlo simulation. Included are the 
average velocity requirements for the primary and extended missions, the 
standard deviations, and, most importantly, the velocity needed, in the 
presence of expected navigation errors, to give a 0.995 probability that the 
velocity required is less than the indicated value. The velocity requirements 
for the 0.995 probability level were derived from the cumulative distribu- 
tions calculated by the Monte Carlo program. The cumulative distributions 
for the primary and extended missions are shown in Figures 4 and 5, 
respectively. 

The velocity requirements given by Table 2 for the primary mission are 
well within the total capability of the spacecraft. This capability will vary 
between approximately 90 and 140 m/s depending upon the final spacecraft 
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Figure 5. Cumulative distribution for the Mariner Venus-Mercury 1973 midcourse 
velocity requirements of the extended mission 



110 


t 








i 


ft 




w_ 




mass. 3 On the other hand, again depending on the final spacecraft mass, the 
capability to accomplish the extended mission will range from marginal to 
relatively certain (subject, as always, to selecting a Mercury return aim point 
at Mercury I). 

Finally, as the discussion of “Maneuvers Required’’ suggested, the 
accurate navigation to the intermediate encounter (Venus for the primary 
mission and Mercury I for the extended mission) is crucial if reasonable 
velocity requirements are to be maintained. This is clearly evident from 
comparison of the post-encounter velocity requirements, maneuvers 3 and 5, 
with the primary and extended mission velocity requirements. This 
conclusion will be reinforced in the section below dealing with navigation 
sensitivities to the contributing error sources. 

Delivery Accuracy 

The purpose of the Mariner Venus-Mercury 1973 mission is to perform 
certain scientific investigations near the target planet. This intent gives rise 
to the definition of a guidance success zone at the taiget planet(s), the 
implication being that if the spacecraft passes through this zone the 
scientific investigations can be performed satisfactorily. In the case of 
Mariner Venus-Mercury 1973, the guidance success zone at Venus is a small 
region about that point for which Mercury encounter would be achieved 
without further corrections. 

The guidance success zone for Mercury I is shown by Figure 6. The 
boundary is labeled with ihe iuuhcs ot tne instrument(s) determining the 
particular segment of the boundary. Two zones are shown with relative 
“worth” values as indicated. Also shown is a typical three-sigma delivery 
dispersion ellipse, describing the delivery capability after performing four 
maneuvers. Clearly, there is a significant probability of passing within the 
90% value guidance success zone. This is confirmed by Table 3, which shows 
the probability of guidance success following the indicated number of 
maneuvers. It should also be noted that the probability is nearly unity for 
passing through the 50% science value zone. 

Currently, no guidance success zone has been defined for Mercury II. 
Consequently, a probability of guidance success at Mercury II cannot be 
calculated at this time. 


Sensitivity Analysis 

In this section the sensitivities of the velocity requirements and delivery 
capability to the various error sources are presented. Information of this 


3 Variations in spacecraft mass alone are not responsible for the relatively large range of 90- 
140 m/s in maneuver capability. Rather, because of costs i. volved in calibrating the 
blowdown propulsion subsystem, the calibrations will be made with propellant tanks 
approximately two-thirds full and full. Then, for a given spacecraft dry mass, the tanks (at 
launch) will be loaded to full or two-thirds capacity in order to stay within the launch vehicle 
payload capability. 
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Table 3. Nominal Mariner Venus-Mer'ury 1973 probabilities of guidance success 


Number 

of 

Maneuvers 

Probability of guidance 

success 

At Venus 

At Mercury I 
(90% value zone) 

1 

0.28 

-0 

2 


<3 X 10" 4 

3 

Not 

0.05 


applicable 


4 

Not 

0.91 


applicable 
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Figure 6. Mariner Venus-Mercury 1973 Mercury aiming zones 


type will help in developing an understanding of the interplay among the 
various navigation elements by showing the effects of changes in element 
capability on the overall navigation capability. 
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To establish these sensitivities, the Monte Carlo simulation was per- 
formed for separate increases of 20% in the expected error of each of the 
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follov lajor cor-ributors to the total navigation capability: orbit 
''V’oi r ' < at er x maneuver time, maneuver thrust vector pointing, and 

ancuvcr imp * ontrol. The pointing and impulse control describe the 
«xe» t ''i. o r >r sources. 

The results of the sensitivity study are given by Table 4, where the change 
in t’ e mean value of the velocity requirement for each maneuver and for 
both the primary and extended mission are shown for each of the indicated 
changes. Inspection of this table shows thr "-‘flowing: 

(1) As expected, changes in the orbit determination capability at or 
following a maneuver epoch do not affect the velocity requirements 
at or preceding that maneuver. 

(2) The velocity requirements are most sensitive to the orbit determina- 
tion capability just preceding a planetary encounter (orbit determina- 
tion at maneuvers 2 and 4). This reflects the fact that the final 
delivery accuracy is dominated by the orbit determination capabil- 
ity. 4 Consequently, it is principally the orbit determination errors 
that experience the magnification effects shown by Figure 2. 

(3) For virtually the same reason as just stated, the delivery accuracy to 
Mercury I is most sensitive to the orbit determination capability at 
the time of the fourth maneuver. 

(4) The velocity requirements and delivery accuracy are relatively 
insensitive to pointing and impulse errors. This, too, reflects the 
comments given for (2) above, since the pointing error and most of 
the impulse error are proportional to the maneuver velocity 
increment. 


Conclusions 

The Mariner Venus-Mercury 1973 spacecraft design provides adequate 
midcourse velocity correction capability for the primary mission. Depend- 
ing upon the eventual propellant loading (see Footnote 3), the midcourse 
velocity correction capability for the extended mission will be either 
marginal or adequate. 

The mission velocity requirements are most sensitive to the delivery 
accuracy to encounters with the intermediate targets, Venus in the case of 
the primary mission, and Mercury I in the case of the extended mission. 
Furthermore, these delivery accuracies are currently limited by the present 
orbit determination capability; however, the orbit determination capability 
is adequate for mission success. 


4 This is the result of (1) performing multiple maneuver j between encounters, thereby reducing 
the contribution of execution errors that are proportional to the maneuver size; and (2) 
performing the final pre encounter maneuver as close as possible to that encounter in order to 
reduce the contribution of execution errors that are independent of the maneuver size. 
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For the primary mission, the probability of requiring at least four 
maneuvers in > /der to accomplish an accurate Mercury encounter is 
approximately 0.95. With four maneuvers, the delivery accuracy to Mercury 
I can be rated as good, since the probability of passing through the guidance 
success zone for the 90% value level is approximately 0.91. 
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Application of New Radio Tracking Data Types to 
Critical Spacecraft Navigation Problems 

V. J. Ondrasik and K. H. Rourke 

Mission Analysis Division 


This article is concerned with newly en/isioned Earth-based radio tracking data 
types involving simultaneous or near-simultaneous spacecraft tracking from widely 
separated tracking stations. These data types are conventional tracking instrumenta- 
tion analogs of the very long baseline interferometry (VLBI) of radio astronomy— 
hence the name quasi- VLBI. Quasi-VLBI promises to help alleviate two particularly 
troublesome problems encountered in interplanetary orbit determination: the zero 
declination and process noise problems. This article motivates quasi-VLBI with a 
preliminary analysis using simplified tracking data models. The results of accuracy 
analysis studies are then presented for a representative mission, Viking 1975. The 
results indicate that, contingent on projected, not overly optimistic future tracking 
system accuracy, quasi-VLBI can be expected to significantly improve navigation 
performance over that expected of conventional tracking data types. 


Introduction 

Two particularly troublesome problems that occur in determining the 
orbit of an interplanetary spacecraft are the zero declination and process 
noise problems. The process noise problem refers to the difficulties 
encountered in determining the orbit of a spacecraft that is subject to 
random nongravitational acceleration uncertainties. The acceleration 
uncertainties, although often negligible in their direct effect on the physical 
orbit of a spacecraft, can severely limit the capability of actually solving for 
the orbit on the basis of conventional tracking data types. The zero 
declination problem is a more familiar difficulty, i.e., the magnification of 
declination errors at low spacecraft declinations. 

In pursuit of new techniques for alleviating the low declination and 
process noise problems, the advantages offered by very long baseline 
interferometry (VLBI) data are being investigated. It was determined that 
the VLBI data does, in concept, provide at least a partial solution to both the 
low declination and process noise problems. Furthermore, it was discovered 
that the information needed for the solution of these two problems may be 



obtained by using sufficiently high quality doppler and range data taken 
simultaneously, or near simultaneously, from widely separated tracking 
stations, without having to resort to the special data taking and reduction 
procedures squired by VLBI. The data quality requirements are well 
within anticipated capabilities which will be associated with the 1-m 
ranging machines, hydrogen masers, and improved transmission media (i.e., 
troposphere, ionosphere, space plasma) calibrations. 

As will be shown later the use of the simultaneous and near-simultaneous 
data can be expected to reduce declination errors that are produced by low 
declination geometries by a factor of 2-4, and to reduce errors in the 
estimate of the spacecraft state (i.e., position and velocity) produced by 
unmodelable accelerations by two orders of magnitude. Because of the 
similarities between VLBI data and simultaneous or near-simultaneous 
doppler and range data, the latter type of data is often called quasi very long 
baseline interferometry data (QVLBI). 

Simultaneous and Near-Simultaneous Data 

Conventional methods of radio interplanetary navigation rely almost 
entirely on the use of the “two-way” coherent doppler and range data types. 
Two-way doppler and range are obtained by a tracking station that 
continuously transmits to a spacecraft which, in turn, receives and 
retransmits to the transmitting station. Transmitted and received signals are 
compared to obtain measurements of the spacecraft’s range rate and range 
with respect to the tracking station. Practical considerations prohibit two- 
way measurements being made by separate stations. Nontransmitting 
stations are, however, free to receive the signal that is being transmitted by 
the spacecraft. This procedure, in principle, provides a measurement of the 
round-trip distance, and its rate of change, between the two tracking stations 
and the spacecraft. Measurements obtained in this way via the separate 
stations are referred to as “three-way” data types. The schematic configura- 
tions associated with two-way and three-way tracking are illustrated in 
Figure 1. (Three-way range has not as yet been used as an explicit data type, 
yet is equivalent to interstation timing techniques that have been used for 
lunar spacecraft tracking. See Reference 1.) 

The three-way tracking data suffers from practical difficulties that have 
heretofore restricted their use from interplanetary applications. The three- 
way measurements are by necessity based on frequency and clock standards 
remote from the transmitting station. The conventional rubidium frequency 
standards, although suitably stable for two-way measurements over moder- 
ate round-trip communication distances, introduce doppler measurement 
errors that preclude the usefulness of the three-way tracking data. The new 
generation hydrogen frequency standards to be installed in the Deep Space 
Network do, however, provide the degree of improvement required to make 
three-way tracking data an attractive interplanetary navigation data type. 

The effective use of three-way range is still hampered by the requirement 
of precise interstation clock synchronization to enable the accurate 
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Figure 1. Transmitting and receiving configuration for two-way and three-way data 
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measurement of the round-trip signal delay. Current microsecond-level 
timing accuracies would introduce excessive three-way range measurement 
errors on the order of 300 m. The clock synchronization problem can be 
circumvented by replacing the simultaneous two-way and three-way range 
measurements with near-simultaneous two-way range measurements be- 
tween the two tracking stations, one just before and one just after an 
interstation handover. In this case the clock discrepancy only introduces 
errors into each measurement epoch specification. This synchronization 
error produces negligible errors, given today’s timing accuracies. This 
technique is limited somewhat by uncertainties present in the motion of the 
spacecraft between the near-simultaneous measurements. For example, if 
the topocentric range rate of spacecraft is known to only 1 mm/s and the 
near-simultaneous measurements are separated by 15 min, then the motion 
uncertainty introduces an error of 

1 mm/s X 900 s = 0.9 m 


Analytical Representations of the Doppler and Range 

Observables 

An insight into how QVLBI techniques help to resolve the low declination 
and process noise problems may be obtained by performing a simple 
analysis based upon analytical representations of the doppler (or, equiva- 
lently, range rate) and range observables. These representations developed 
by Hamilton and Melbourne (Reference 2) are based upon the following 
approximations of the topocentric range, from a tracking station to a distant 
spacecraft: 


p(t) = r(t ) + z sin 6 (Y) - r cos 8(t) cos [co/ + 0 + X -o^r)] 

S S 


( 1 ) 



f 

f 

i 
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Figure 2. Spacecraft and station coordinates 


where (see Figure 2) 

p = topocentric range 
r = geocentric range 
5 = declination 
a = right ascension 

r s = distance of tracking station from the Earth’s spin axis 
z s = distance of tracking station from the Earth’s equator 
X = longitude of tracking station 
<o = Earth’s rotation rate 

<f> = a phase angle that depends on the choice of epoch 
t — time past epoch 

As schematically illustrated in Figure 3, this equation shows that the 
observables are composed of a diurnal sinusoidal motion, provided by the 
Earth’s spin, superimposed upon the slowly varying geocentric motion. The 
term produced by the Earth’s spin is a strong source of angular data and 
usually provides a good estimate of 8 and a within a few days. The slowly 
varying term contains information on all the spacecraft state variables, but 
by itself can seldom provide complete estimates. This is because of 
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Figure 3. Schematic illustration of the doppler observable 


difficulties associated with separating the geocentric motion effects pro- 
duced by specific variations in the state variables. 

A convenient approximation to the doppler observable, which may be 
used to examine the declination and right ascension information contained in 
a single pass of doppler data, is shown in Equation 2, and is obtained by 
assuming that 8 and a are constant: 


p(t)~ 'r (t) + cor s cos 5 Q sin {cot + <f> + X - <* 0 ) (2) 


where 

a da/ dt 


oo = a(t=0) 


The Low Declination Problem 

As mentioned above, when conventional data is used, the determination of 
the declination is made principally through the amplitude of the sinusoidal 
term. The amplitude b may be written as 


b s r„ cos 6 

3 


( 3 ) 


120 




STATION 1 



STATION 2 


Figure 4. Simultaneous ranging from two stations 


Since the data taken during a tracking pass may be thought of as providing a 
measure of b , the orbit determination filter will attempt to leave it 
unchanged. Therefore, an error in r s , A r s , will produce a compensating error 
in the declination, A 5, which for data arcs of a few days may be 
approximated by the following equation: 


A 6 = 


cos 5 
sin 5 



(4) 


The most notable feature of this equation is that there is a singularity at 6 = 
0. Although this singularity does not occur in the actual orbit determination 
procedure, there is almost always a substantial magnification of declination 
errors, as 8 approaches zero. 

An insight into how QVLBI-type data avoids the low declination problem 
may be easily obtained by considering the situation depicted in Figure 4, 
where two stations are making simultaneous range measurements at local 
spacecraft meridian crossing. Assume that the stations share the same 
longitude and have equal but opposite latitudes. From this figure it is easily 
seen that the difference between the two range measurements is related to 
the baseline D and the spacecraft declination through the following 
equation: 



p 2 - pj = D sin <5 


( 5 ) 


Clearly, range-type errors, A p, whether they are the result of actual 
measurement errors or station location errors, will give declination errors of 
the approximate amount 


AS = — — ^ 
cos 6 D 


( 6 ) 





This equation, unlike Equation 4, is well behaved at low declination and 
illustrates that, given data of sufficient quality, simultaneous (differenced or 
nondifferenced) range data from widely separated tracking stations provides 
a superior determination of the declination than does the conventional range 
and doppler data. 


The Process Noise Problem 

This section describes how random accelerations actually affect orbit 
determination solutions, and provides the motivation for the use of QVLBI 
data types as a countermeasure to the process noise problem. 

This discussion is concerned with the short-term process noise effects and 
their influence on a single day’s pass of conventional doppler data. The 
short-term effects of process noise on the orbit determination of a distant 
spacecraft can be illustrated with the aid of the Hamilton/ Melbourne 
doppler data representation as shown in Equation 2. 

Consider particular corrections A r, Aa, and AS to the spacecraft 
geocentric range rate, right ascension, and declination, respectively. They 
produce the following range-rate (doppler) residual signature: 

Ap - Ar - (r to sin 5 sin co/)AS 

S 

+ (r to cos 5 cos cot)A8, - n/2 < to/ < n/2 (7) 

a 

over the day’s pass of data. Normally, Ar, AS , and Aa are considered to be 
constant over one day (as well as S, the nominal spacecraft declination). The 
effect of particular error sources can be expressed in terms of their 
corresponding doppler residuals. Assume that Ar is not constant but 
perturbed by random accelerations. The perturbation produces the stron- 
gest effect on Acr and AS determinations if 


Ar ~ Ar 0 + a sin cot + b cos ajt 


1 


* 


.? 

j 

* 

t 

* 
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i.e., if the error residual is sinusoidal with a one-day period. For example, a 
radially acting acceleration of 10 -12 km/s 2 (typical for a ballistic spacecraft) 
varying as 


(1CT 12 km/s 2 ) X cos c ot 


produces an effect indistinguishable from a A 8 perturbation on the order of 
(10~ ,2 )/r s ca 2 , nearly 0.04 second of arc. The resulting position error for a 
spacecraft at 2 AU from the Earth is 60 km. The actual net movement of the 
spacecraft is zero; its maximum perturbed displacement is 0.2 m. Thus, 
conventional doppler tracking data can be expected to be very sensitive to 
uncertainties in the geocentric radial motion of the spacecraft. Uncertainties 
in the Aa and A8 directions caused by process noise can be shown to be 
negligible for a single pass of doppler data. 

The influence of process noise on Ar estimation errors is comparable in 
magnitude to the corresponding errors in r s coAa and r s coAS. The direct 
effects of the Ar errors are not significant; e.g., a large (with respect to 
typical acceleration uncertainties) 10-mm/s Ar error produces only a 10-km 
range error in 10 days. However, given purely ballistic motion, considerable 
information concerning the dynamical state of the spacecraft is available 
from the variation of Ar over several days. Thus, process noise sensitivities 
similar to those occurring for a single pass of data can be expected over the 
longer data arcs. This phenomenon is particulary evident for orbit 
determination data processing methods that ignore the effects of random 
acceleration uncertainties. This point is developed in considerable detail in 
Reference 3. 

Differenced Tracking Data as a Countermeasure for the 

Process Noise Problem 

The preceding section serves as an introduction to the use of differenced 
simultaneous or near-simultaneous tracking data from two widely separated 
tracking stations as a countermeasure to the process noise problem. The 
differencing of the data types is effective because it eliminates the 
troublesome geocentric corruptions common to both tracking stations, yet 
preserves the “topocentric” information arising from the two stations’ 
rotation. Consider, then, the range-rate residuals observed from two 
separated tracking stations: 

Apj(0 = Ar(t) + b x sin cot +Cj cos cot ) 

} ( 8 ) 

Ap 2 (/) = Ar (/) + b 2 sin cot + c 2 cos cot) 


The parameters b lt b 2 , Cj, c 2 are linear expressions in the Aa and AS 
corrections, their explicit form depending on the particular time reference 
used in the above representations. Two-way doppler residuals obtained at 
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Station 1 can be expressed as 2Apj. Three-way doppler residuals available at 
Station 3 are of the form Api + Ap 2 - cAf/f where the cAf/f term arises 
from the frequency standard discrepancy Af between Stations 1 and 2. The 
difference of two-way doppler from Station 1 and three-way doppler from 
Station 2 is represented as 


Ap i - Ap 2 + 


cAf _ cAf 

T~ = 7~ 


+ (b j -ft 2 )sin c of + (Cj -c 2 )cos cjt 


(9) 


over the stations’ view period overlap tft 2 < o)t < The geocentric range- 
rate terms subtract out and are replaced by a “velocity bias” cAf/f arising 
from the relative station-to-station frequency standard bias Af If ( c = speed 
of light). Thus, in the presence of large unmodelable random accelerations, 
the differenced doppler allows separation of A 8 and Aa determination, 
through b\ - b 2 and C\ - c 2 , from a corrupted Ar determination. The 
technique is hindered, however, by the introduction of a velocity bias 
uncertainty in the place of the geocentric range-rate uncertainty. Clearly, 
the differenced doppler data can be effective in circumventing process noise 
effects only as long as the uncertainties arising from frequency standard 
instability are significantly less than the process noise uncertainties expected 
in the conventional doppler data. 

In the section on The Low Declination Problem , the use of simultaneous 
or near-simultaneous range measurements was discussed with regard to 
determining spacecraft right ascension and declination angles. Since these 
measurements allow short arc solutions of equatorial angles, they exhibit the 
same advantages of reduced sensitivity to short-term process noise effects as 
does the differenced doppler data. In contrast to doppler measurements, 
he -'er, the range measurements permit declination solutions at low 
Spacecraft declination. 

This discussion is not intended to imply that the use of differenced two- 
way and three-way doppler or differenced range data is necessarily an 
efficient use of the data received at both stations from the spacecraft. The 
differenced data is effective in allowing separation of topocentric and 
geocentric tracking information— even in the case of a spacecraft experienc- 
ing large random accelerations. Ultimately, maximum information is 
extracted if concurrent two-way and three-way data are processed together 
with a suitably designed orbit determination filter that takes advantage of 
the known random acceleration characteristics. The differenced data 
provides, nevertheless, a simple conceptualization for preliminary analysis, 
as well as a straightforward and probably adequate approximation to an 
“optimal” treatment of concurrent two-way and three-way data. 

Example Mission 

In the preceding section and in the section on The Low Declination 
Problem, simplified representations of the range and doppler data were used 
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to illustrate how QVLBI data may help resolve the low declination and 
process noise problems. However, to quantitatively describe the utility of 
the QVLBI ;lata, detailed accuracy analyses are required. Such a study was 
performed for the Viking Mission B trajectory described in Table 1, using 
data taken during the tracking pattern shown in Figure 5. The study 
consisted primarily of performing conventional weighted least-squares batch 
filter solutions to examine the effects that data noise, station location errors, 
and process noise have on the estimates of the spacecraft state at epoch. For 
comparative purposes, parallel solutions were made with the data-taking 
strategies described in Table 2. 

These data strategies include: (1) a conventional single-station strategy 
using nominal accuracy range-rate data, with ranging data weighted at 100 
m; (2) single-station range and range-rate tracking with tight, 3-m range 
weight specification; (3) single-station range rate with 3-m simultaneous 
range measurements; and (4) differenced simultaneous range rate and range, 
QVLBI data, augmented with loosely weighted conventional range and 
range-rate measurements. 

The above-mentioned data weights are data accuracy specifications 
provided to the orbit determination filter, and should be distinguished from 
expected measurement precision. For instance, range data of a given quality 
may be assigned different weights, depending on how well the spacecraft/ 
tracking system environment can be modeled. 

The results of this study are summarized in Figure 6, which shows the 
standard deviations in 5, a, 5, and a produced by data noise, and errors in 
* tese same quantities produced by the indicated station location errors and 
process noise. The standard deviation and errors of each component are ihen 
root-sum-squared to obtain the consider standard deviation <r c , which is also 


Table 1. Description of Viking Mission B trajectory 


Coordinate 

Value at 

encounter* - 30 days 

Value at 

encounter* - 9 days 

r, km 

0.3300 X 10 9 

0.3489 X 10 9 

6, deg 

11.44 

6.853 

a. deg 

153.8 

165.8 

r, km/s 

0.1142 X 10- 2 

0.8390 X 10 

• 

6, rad/s 

-0.3952 X 10- 7 

-0.4441 X 10' 7 

a, rad/s 

0.1109 X 10' 6 

0.1112 X 10' 8 

a Encounter 

- 23:45:00, August 6, 1976. 
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Figure 5. Tracking pattern 


shown. The corresponding r 0 and r {) errors are not shown because they are 
always less than errors in the 8 and a directions. 

Figure 6 vividly illustrates some of the problems associated with 
determining the orbit of a spacecraft at low declinations using convention d 
doppler and range (<r p =* 100 m) data. Of particular note is the inability to 
solve for 8, and the sensitivity of 8 to errors in station locations as shown by 
the 500-km declination direction error. The solution is improved by 
approximately 30% if the weight on the range data taken at a single station is 
changed to 3 m. A dramatic improvement in the solution is provided by the 
addition of simultaneous 3-m weighted range points from Australia. 
Simultaneous range data reduces the station location error sensitivity by 
factors of three or five over the two single-station data sets and increases the 
ability to solve for 8 by a factor of two over the conventional solution. 

Although the simultaneous range data dramatically reduces the solution 
errors produced by station location errors, the requirement of tight range 
weights substantially increases the importance of process noise as an error 
source. Frocess noise sensitivities for the various data sets are shown in 
Figure 6. The analysis technique used to evaluate the process noise effects is 
that developed by Curkendall in Reference 3. The technique enables 
evaluation of suboptimai, weighted least-souares data filter performance f< r 
particular process noise environments. The process noise is modeled with 
adequate generality as a three-axis, exponentially time-correlated stationary 
stochastic process. The effect of a range of noise correlation times is shown 
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Figure 6. Position ard velocity errors resulting from the use of conventional, simultaneous range, and QVLBI data 





with the bar shading— the lower level corresponding to a relatively 
optimistic bias noise assumption (infinite correlation time), the upper level 
corresponding to a pessimistic, rapidly varying noise assumption (5-day 
correlation time). Although the noise magnitude is given as 10~ 12 km/s 2 , the 
results are completely general inasmuch as the process noise errors are 
scalable with respect to the noise magnitude. 

Data set 4 includes, in place of the precise simultaneous range (and range 
rate) measurements, precise differenced simultaneous range and range 
measurements. Thus, although some available geocentric information is 
deleted, the process noise sensitivity should be reduced. The actual results 
are remarkable. Process noise declination errors are reduced from 1000-km 
levels to 3 and 2 km for t = oo and r — 5, respectively. Similar reductions 
are observed for declination velocity as well as for the right ascension 
direction errors. The geocentric range and range-rate accuracies are 
commensurate with the data set 4 conventional data type accuracies. 

A criticism of the above analysis is that the nondifferenced data is treated 
unfairly, in that an attempt to account for process noise through more 
modern filtering techniques can be expected to result in significant 
reductions in process noise sensitivity. Nevertheless, the importance of 
simultaneous range measurement for improving declination solutions cannot 
be minimized, given that proper precautions are taken against process noise 
sensitivity— either in the form of explicit data differencing or “optimal” 
processing. “Optimal” processing, by recovering some geocentric informa- 
tion, can be expected to outperform suboptimally processed differenced 
data, although probably marginally, and at the expense of some sensitivity to 
the process noise assumptions used in the filter design. 

QVLBI Error Sources 

Although QVLBI data in certain circumstances is inherently superior to 
conventional data, its practicable superiority depends heavily upon the 
quality of the data. 

The crucial accuracies affecting the feasibility of effective simultaneous 
doppler and range measurements lie in the media and instrumentation error 
categories. The primary limitation to simultaneous doppler measurements 
lies in the three-way doppler sensitivity to local oscillator stability. As 
mentioned previously, feasible application of three-way doppler measure- 
ments requires hydrogen maser frequency standards. Adequate estimates of 
simultaneous range accuracies are difficult to obtain at this time since meter- 
level ranging accuracies have heretofore been considered unnecessary. Thus, 
current specifications are expected to be overly pessimistic with regard to 
differenced range applications. Table 3 presents media calibration and 
instrumentation accuracies for both simultaneous and near-simultaneous 
QVLBI techniques. In light of the uncertainty regarding the actual possible 
accuracies, several values are quoted for each error source, including 
expected present capability, and upper and lower values for projected 
future capability. The future quotations include earliest availability dates. 
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The projected accuracies are illustrated in Figure 4. (The future capabilities 
are drawn from Reference 4 that addresses several specific quasi-VLBI 
configurations.) These values indicate that meaningful demonstrations of 
differenced range techniques can be conducted at present, and that the 
future goal of 1-m level differenced range measurements is indeed a 
plausible one. Nevertheless, the advantages of differenced range techniques 
will become clearer with a more detailed analysis and with experimental 
verification of (1) range instrument capabilities in the 3-m domain, and (2) 
sub-meter charged particle and troposphere calibration capability. 
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For abstract, see Mulholland, J. D. 


AUSMAN, N. E., JR. 

A02 Simulation of Mariner Mars 1971 Spacecraft 

N. E. Ausman, Jr., N. K. Simon, and C. F. Rodriguez 

JPL Quarterly Technical Review, Vol. 1, No. 3, pp. 67-78, 
October 1971 

In preparation for the Mariner Mars 1971 mission, operations 
personnel took part in an extensive training program during which 
the primary source of spacecraft data was a computer program 
simulating the spacecraft. The objectives of a simulation model for 
training purposes differ from objectives appropriate to a design or 
analysis model. Model subsystems were designed to provide realis- 
tic telemetry data reflecting changes due both to commands and 
environmental parameters affecting the spacecraft at various times 
during the mission. 
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The spacecraft is modeled along two separate functional lines. 
Boolean operations are concentrated in the spacecraft logic model, 
which determines the spacecraft state or mode, while mathemati- 
cal operations or algorithms are executed in computational subsys- 
tem models. Although logic parameters are interrogated as a part 
of each computational pass, actual logic model processing occurs 
only when a change-of-state input is generated by the operations 
organization. This article discusses the program design, some of 
the special characteristics of each of the modeled subsystems, and 
how the model was used in support of mission operations training. 


BASSET, R. S. 

B01 DSN Progress Report for September-October 1971: DSN 
Telemetry System Tests 

W. J. Kinder and R. S. Basset 

Technical Report 32-1526, Vol. VI, pp. 10-12, 
December 15, 1971 

For abstract, see Kinder, W. J. 


BATELAAN, P. D. 

B02 DSN Progress Report for September-October 1971: Waveguide 
Voltage Reflection Calibrations of the MXK Cone 
(Modification 1) 

P. D. Batelaan 

Technical Report 32-1526, Vol. VI, pp. 123-124, 

December 15, 1971 

A listing of reflection coefficients for the multiple-frequency X- 
and K-band (MXK) cone (Modification 1), at both X- and K-band 
frequencies, is presented. Included is a short discussion of the 
measurement technique and results. These measurements were 
made on the ground at JPL and the Goldstone Deep Space Com- 
munication Complex during the final assembly stages of the up- 
dated cone. 


BEER, R. 

B03 Absorption by Venus in the 3-4-Micron Region 

R. Beer, R. H. Norton, and J. V. Martonchik (University of 
Texas) 
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Astrophys. J., Vol. 168, No. 3, Pt. 2, pp. L121-L124, 
September 15, 1971 


The continuum absorption by Venus in the 3-4- (i region has been 
determined from a detailed comparison of high-resolution spectra 
of Venus and the Sun. The resultant ratio shows a major depression 
centered near 2580 cm -1 and a less well-defined dip near 3050 
cm -1 . The possibility that these may be due to bicarbonates in the 
upper clouds of Venus is discussed. 


B04 Astronomical Infrared Spectroscopy With a Connes-Type 
Interferometer: II. Mars, 2500-3500 cm 1 

R. Beer, R. H. Norton, and J. V. Martonchik (University of 
Texas) 

Icarus: Int. J. Sol. Sys., Vol. 15, No. 1, pp. 1-10, 

August 1971 

New spectra of Mars in the 3- to 4-/im region at significantly 
higher resolution than previously available were obtained near the 
1969 opposition. No features positively identifiable as being due to 
the Martian atmosphere could be detected. The existence of an 
albedo drop, probably due to surface water of hydration, is 
confirmed. 


BENDER, P. L. 

B05 Preliminary Results of Laser Ranging to a Reflector on the 
Lunar Surface 

J. D. Mulholland, C. 0. Alley (University of Maryland), 

P. L. Bender (Joint Institute for Laboratory Astrophysics), 
D. G. Currie (University of Maryland), 

R. N. Dicke (Princeton University), J. E. Faller (Wesleyan 
University), W. M. Kaula (University of California, Los 
Angeles), G. J. F. MacDonald (University of California, 
Santa Barbara), H. H. Plotkin (Goddard Space Flight 
Center), and D. T. Wilkinson (Princeton University) 

Space Research XI, pp. 97-104, Akademie-Verlag, Berlin, 
1971 

For abstract, see Mulholland, J. D. 
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BERGSTRALH, J. T. 


B06 Recomputation of the Absorption Strengths of the Methane 
3v 3 J-Manifolds at 9050 cnr 1 

J. T. Bergstralh (McDonald Observatory) and J. S. Margolis 

J. Quant Spectrosc. Radiat. Transfer, Vol. 11, No. 8, 
pp. 1285-1287, August 1971 

Absorption intensities of the R-branch /-manifolds (0 < / < 7) of 
the 3v 3 methane band have been computed using half-widths 
measured for the individual components of the /-manifolds. The 
effects of small pressure shifts of the lines are included. The results 
yield a value for the rotational temperature in agreement with the 
laboratory temperature and are very close to the results obtained 
assuming no / dependence of the half-widths. 


BERMAN, P. A. 

B07 Effects of Storage Temperatures on Silicon Solar Cell Contacts 

P. A. Berman and R. K. Vasui 

Technical Report 32-1541, October 15, 1971 

Solar cells having various contact systems and configurations were 
investigated to determine the effect of storage temperature and 
duration on cell electrical and mechanical characteristics. Cells 
having n diffused into p-base silicon, p diffused into n-base silicon, 
and p diffused into n-base lithium-containing silicon were studied. 
Contact systems of silver-titanium, silver-titanium with solder 
coating, silver-titanium with palladium, and electroless nickel 
were investigated. Also included in the study were submodules 
similar to those used in past JPL flight projects. 

The results showed that electrical and mechanical degradations 
previously observed in high-humidity, high-temperature environ- 
ments were the result of the combined environments and not due 
to the temperature component alone. Solder-coated silver-tita- 
nium contacts can be adversely affected by prolonged exposure to 
150°C temperature soak and to -196°C liquid nitrogen exposure. 
The lithium-containing cells exhibited good electrical and mechan- 
ical contact stability with respect to the other cells tested at 
150°C. The addition of palladium to the silver-titanium system 
presented no apparent advantage for the 150°C temperature soak, 
except in the case of relative top contact strength, where stability 
(but not absolute pull strength) was superior to that for the other 
cell types. 
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BOS Supporting Data Package for TR 321541, Effects of Storage 
Temperatures on S;:xon Solar Ceil Contacts 

P. A. Berman and R. K. Yasui 

Technical Memorandum 33-497, October 15, 1971 

This data package is a companion document to JPL Technical 
Report 32-1541 and contains a series of summary curves for cells 
having various contact systems, depicting the effects of tempera- 
ture exposure as a function of time on the electrical and mechani- 
cal characteristics. Each curve represents the results obtained from 
a minimum of 32 and a maximum of 1104 individual data inputs 
and includes the 95% confidence limits associated with each time- 
temperature combination for which measurements were obtained. 
The curves are useful for detailed analysis and comparison of the 
behavior of the contact systems as a function of the environmental 
conditions studied. 


BERWIN, R. 

B09 DSN Progress Report for July-August 1971: Superconducting 
Magnet for a Ku-Band Maser 

R. Berwin, E. Wiebe, and P. Dachel 

Technical Report 32-1526, Vol. V, pp. 109-114, 

October 15, 1971 

A superconducting magnet to provide a uniform magnetic field of 
up to 8000 G in a 1.14-cm gap for the 15.3-GHz (Ku band) 
traveling wave maser is described. The magnet operates in a 
persistent mode in the vacuum environment of a closed-cycle 
helium refrigerator (4.5 K). The features of a superconducting 
switch, which has both leads connected to 4.5 K heat stations and 
thereby does not receive heat generated by the magnet charging 
leads, is described. 


BREMNER, D. S. 

B10 DSN Progress Report for July-August 1971: DSS 

Communications Equipment Subsystem Simulation Center 
High-Speed Data Assembly 

D. S. Bremner 

Technical Report 32-1526, Vol. V, pp. 132-135, 

October 15, 1971 

The 1971-1972 era required expansion of the Simulation Center 
High-Speed Data Assembly involving extensive modifications. An 
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increase in the number of channels was necessary to provide 
simultaneous data handling configurations. New data sets were 
installed to double the data rate as required by the DSN-GCF 
High-Speed Data System. 


CAMERON, R. E. 

C01 Farthest South Soil Microbial and Ecological Investigations 

R. E. Cameron, G. H. Lacy, F. A. Morelli, and 
J. B. Marsh (University of California, Davis) 

Antarctic J. U.S., Vol. VI, No. 4, pp. 105-106, 

July-August 1971 

In austral summer 1970-1971, soil samples were collected asepti- 
cally from the surface to the depth of hard permafrost or bedrock 
at sites near Mount Howe, Antarctica. Environmental measure- 
ments were made for soil and microclimatic characteristics. Air 
samples were also taken with liquid, semisolid, and high volume 
dry membrane-type air samplers. Soil properties determined by 
methods used for desert soil and other geologic material are given. 
Soil and air temperatures taken at one site are also given. 


CAMPBELL, J. K. 

C02 A Method of Orbit Determination Using Overlapping Television 
Pictures 

J. K. Campbell 

J. Spacecraft Rockets, Vol. 8, No. 8, pp. 867-872, 

August 1971 

A method that utilizes the Mariner Mars 1969 television pictures 
has been developed to aid in the post-encounter determination of 
the TV viewing directions. A weighted least-squares estimation 
scheme is used, based on observation of common Mars surface 
features in overlapping TV pictures. Estimates were also obtained 
for the encounter orbits of Mariners 6 and 7, and also for the 
orientation of the Mars spin axis; however, the estimated solutions 
for the spacecraft orbits and for the spin axis did not result in a 
reduction of the a priori uncertainties in these parameters. Good 
results were obtained for the estimated corrections and the a 
posteriori uncertainties of the nominal TV viewing directions 
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CANNON, W. A. 


C03 Lunar Fines and Terrestrial Rock Powders: Relative Surface 
Areas and Heats of Adsorption 

F. P. Fanale, D. B. Nash, and W. A. Cannon 

J. Geophys. Res., Vol. 76, No. 26, pp. 6459-6461, 
September 10, 1971 

For abstract, see Fanale, F. P. 





CARTWRIGHT, D. C. 

C04 Differential and Integral Cross Sections for the Electron-Impact 
Excitation of the a^g and b^g* States of 0 2 

S. Trajmar, D. C. Cartwright (The Aerospace Corporation), 
and W. Williams 

Phys. Rev., Pt. A: Gen. Phys., Vol. 4, No. 4, pp. 1482 -1492, 
October 1971 

For abstract, see Trajmar, S. 


CHAO, C. C. 

C05 DSN Progress Report for September-October 1971: 

Tropospheric Range Effect Due to Simulated Inhomogeneities 
by Ray Tracing 

C. C. Chao 

Technical Report 32-1526, Vol. VI, pp. 57-66, 

December 15, 1971 

A simple ray trace method is developed to study the effect in 
range correction of a radio wave passing through the troposphere 
with inhomogeneities. Inhomogeneities were simulated from previ- 
ous observations. The uncertainties in range correction due to the 
unmodeled horizontal gradient in refractivity are mostly less than 
1% for elevation higher than 5 deg. The average uncertainty due to 
local inhomogeneities, based on nine simulated cases, is below 1-m 
in one-way range correction for elevation angles greater than 5 
deg. 

C06 DSN Progress Report for September-October 1971: New 

Tropospheric Range Corrections With Seasonal Adjustment 

C. C. Chao 
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Technical Report 32-1526, Vol. VI, pp. 67-82, 
December 15, 1971 


A study of two years’ radiosonde balloon measurements shows that 
most of the significant seasonal fluctuations in refractivity profiles 
are due to the variation of the water content in the troposphere. 
At each of the six weather stations, the long-term seasonal varia- 
tion in refractivity profile repeats quite well through the two-year 
span. Based on the two years’ data, a new tropospheric range 
calibration method using monthly mean parameters at each station 
was developed. 

The uncertainty of range correction with the new model over one 
pass is estimated to be 0.30-0.35 m for a 10-deg minimum eleva- 
tion angle and 0.40-0.50 m for a 6-deg minimum elevation angle. 
This is below the required accuracy for tropospheric calibration of 
both the Mariner Mars 1971 and Viking 1975 missions. 


CHELSON, P. 0. 

C07 Reliability Computation Using Fault Tree Analysis 
P. 0. Chelson 

Technical Report 32-1542, December 1, 1971 

A method is presented for calculating event probabilities from an 
arbitrary fault tree. The method includes an analytical derivation 
of the system equation and is not a simulation program. Systems 
that incorporate standby redundancy can be handled with the 
method. Conditional probabilities are used for computing fault 
trees where the same basic failure appears in more than one fault 
path. 

C08 Reliability Computation From Reliability Block Diagrams 
P. O. Chelson 

Technical Report 32-1543, December 1, 1971 

A method and a computer program are presented to calculate 
probability of system success from an arbitrary reliability block 
diagram. The class of reliability block diagrams that can be han- 
dled include any active /standby combination of redundancy, and 
the computations include the effects of dormancy and switching in 
any standby redundancy. The mechanics of the program are based 
on an extension of the probability tree method for computing 
system probabilities. 
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Program Listing for Fault Tree Analysis of JPL Technical 
Report 32*1542 

P. 0. Chelson 

Technical Memorandum 33-512, December 1, 1971 

This memorandum presents the computer program listing for the 
fault tree analysis of Technical Report 32-1542, Reliability Com- 
putation Using Fault Tree Analysis, Jet Propulsion Laboratory, 
Pasadena, Calif., Dec. 1, 1971. The program is written in FOR- 
TRAN V and is currently running on a UNIVAC 1108. 

CIO Program Listing for the Reliability Block Diagram Computation 
Program of JPL Technical Report 32*1543 

P. 0. Chelson and R. E. Eckstein 

Technical Memorandum 33-513, December 1, 1971 

This memorandum presents the computer program listing for the 
reliability block diagram computation program described in Tech- 
nical Report 32-1543, Reliability Computation From Reliability 
Block Diagrams, Jet Propulsion Laboratory, Pasadena, Calif., Dec. 

I, 1971. The program is written in FORTRAN V and is currently 
running on a UNIVAC 1108. Each subroutine contains a descrip- 
tion of its function. 

CHEN, C. J. 

Cll Transition Probabilities for Ar I 

C. J. Chen 

J. Opt. Soc. Am., Vol. 61, No. 9, pp. 1267-1268, 

September 1971 

Interest in the spectra of the rare gases has been renewed in 
connection with laser kinetic studies and other endeavors. Because 
the rare gases show large departures from L-S and jl coupling, 
intermediate-coupling schemes are required to calculate the transi- 
tion probabilities of those gases. This article reports the results of 
calculations for some Ar I arrays of special current interest. The 
parameters adopted for the various configurations, together with 
the final transition probabilities, are included. 


CHOI, S. D. 

C12 High-Power Microstrip RF Switch 

S. D. Choi 
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JPL Quarterly Technical Review, Vol. 1, No. 3, pp. 110-124, 
October 1971 

A microstrip-type single-pole double-throw (SPDT) switch whose 
RF and bias portions contain only a metallized alumina substrate 
and two PIN diodes has been developed. It is superior to electro- 
mechanical and currently used circulator-type switches in many 
aspects of flight-qualified switch characteristics, such as power 
drain, weight, volume, magnetic cleanliness, cost, and reliability. A 
technique developed to eliminate the do blocking capacitors 
needed for biasing the diodes is described. These capacitors are 
extra components and could lower the reliability significantly. 

An SPDT switch fabricated on a 5.08- X 5.08- X 0.127-cm (2- X 2- 
X 0.050-in.) substrate has demonstrated an RF power-handling 
capability greater than 50 W at S-band. The insertion loss is less 
than 0.25 dB and the input-to-off port isolation is greater than 36 
dB over a bandwidth larger than 30 MHz. The input voltage 
standing- wave ratio is lower than 1.07 over the same bandwidth. 
Theoretical development of the switch characteristics and experi- 
mental results, which are in good agreement with theory, are 
presented in this article. 

CLAUSS. R. C. 

C13 DSN Progress Report for July-August 1971: Tracking and 
Data Acquisition Elements Research: Low Noise Receivers: 
Microwave Maser Development 

R. C. Clauss and R. Quinn 

Technical Report 32-1526, Vol. V, pp. 102-108, 

October 15, 1971 

A traveling-wave maser, tunable from 14.3 to 16.3 GHz, has been 
completed and is ready for installation on the 64-m antenna at 
Goldstone Deep Space Communication Complex. The maser can 
provide more than 30 dB net gain at any frequency within its 
tuning range; an equivalent input noise temperature of 8.5 K has 
been measured in the laboratory. The maser is a ruby-loaded comb 
structure (C-axis orientation 90 deg) which operates in a closed- 
cycle helium refrigerator. The 8000-G magnetic field required for 
maser operation is supplied by a superconducting magnet. The 
entire package weight is 70 kg, and the unit is capable of opera- 
tion in any position. 

C14 DSN Progress Report for September-October 1971: Microwave 
Maser Development 

R. C. Clauss and H. F. Reilly, Jr. 
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Technical Report 32-1526, Vol. VI, pp. 118-122, 
December 15, 1971 


A traveling- wave maser operating at 15.3 GHz has been used to 
test the noise temperature contribution of various waveguide com- 
ponents. An assembled system, consisting of the maser, a direc- 
tional coupler, a waveguide switch, a polarizer, and a feed horn, 
measured 23-K total system operating noise temperature. Previ- 
ously measured X-band data are shown for comparison. The maser 
was installed on the 64-m antenna at the Goldstone Deep Space 
Communication Complex. Maximum changes of 2.2-deg signal 
phase and 0.25-dB signal amplitude were observed during antenna 
motion tests for maser phase and gain stability. The excellent 
stability performance is attributed to the use of a superconducting 
maser magnet. 


COYNER, J. V., JR. 

C15 Redial Rib Antenna Surface Deviation Analysis Program 

J. V. Coyner, Jr. 

Technical Memorandum 33-518, December 15, 1971 

As described in this memorandum, a digital computer prograr l has 
been developed that analyzes any radial rib antenna that has ribs 
radiating from a center hub. The program has the capability to 
calculate the antenna surface contour (reversed pillowing effect), 
calculate the optimum rib shape that minimizes the RMS surface 
error, calculate the actual RMS surface error, compensate for rib 
deflection due to mesh tension and catenary cable tension, and 
determine the pattern from which the mesh gores are cut. 


CURRIE, D. G. 

C16 Preliminary Results of Laser Ranging to a Reflector on the 
Lunar Surface 

J. D. Mulholland, C. 0. Alley (University of Maryland), 

P. L. Bender (Joint Institute for Laboratory Astrophysics), 
D. G. Currie (University of Maryland), 

R. H. Dicke (Princeton University), J. E. Faller (Wesleyan 
Universit>. /( W. M. Kaula (University of California, Los 
Angeles), G. J. F. MacDonald (University of California, 
Santa Barbara), H. H. Piotkin (Goddard Space Flight 
Center), and D. T. Wilkinson (Princeton University) 
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Space Research XI, pp. 97-104, Akademie-Verlag, Berlin, 
1971 

For abstract, see Mulholland, J. D. 

DACHEL, P. 

D01 DSN Progress Report for July-August 1971: Superconducting 
Magnet for a Ku-Band Maser 

R. Berwin, E. Wiebe, and P. Dachel 

Technical Report 32-1526, Vol. V, pp. 109-114, 

October 15, 1971 

For abstract, see Berwin, R. 


DALLAS, S. S. 

D02 DSN Progress Report for July-August 1971: A Comparison of 
Cowell’s Method and a Variation-of-Parameters Method for 
the Computation of Precision Satellite Orbits 

S. S. Dallas and E. A. Rinderle 

Technical Report 32-1526, Vol. V, pp. 74-78, 

October 15, 1971 

A precision special perturbations program that uses either Cow- 
ell's method or a variation-of-parameters method to compute an 
elliptical orbit is analyzed to determine which mode is more 
efficient when computing satellite orbits. The results obtained 
indicate that the variation-of-parameters mode is significantly 
more efficient if the numerical integrator being used is optimized 
in that mode by varying the integration order and local error 
control and by using either a predictor or predictor-corrector 
algorithm. 


DAVIS, E. K. 

D03 DSN Progress Report for July-August 1971: Viking Mission 
Support 

E. K. Davis 

Technical Report 32-1526, Vol. V, pp. 24-28, 

October 15, 1971 

Since the redirection of the Viking Project in January 1970, the 
DSN interface organization has been heavily involved with project 
organizational elements in advanced planning, exchanges of techni- 
cal information, identification of requirements, capabilities, prob- 
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lems, and their resolution. This article is a general summary of the 
accomplishments in these areas cf long-range planning including 
pertinent open questions to be resolved in the detailed planning 
phase. 


DICKE, R. H. 

D04 Preliminary Results of Laser Ranging to a Reflector on the 
Lunar Surface 

J. D. Mulholland, C. 0. Alley (University of Maryland), 

P. L. Bender (Joint Institute for Laboratory Astrophysics), 
D. G. Currie (University of Maryland), 

R. H. Dicke (Princeton University), J. E. Faller (Wesleyan 
University), W. M. Kaula (University of California, Los 
Angeles), G. J. F. MacDonald (University of California, 
Santa Barbara), H. H. Plotkin (Goddard Space Flight 
Center), and D. T. Wilkinson (Princeton University) 

Space Research XI, pp. 97-104, Akademie-Verlag, Berlin, 
1971 

For abstract, see Mulholland, J. D. 


DON, J. P. 

D05 Application of Hybrid Propulsion Systems to Planetary 
Missions 

J. P. Don and R. L. Phen 

Technical Memorandum 33-483, November 1, 1971 

The feasibilit and application of hybrid rocket propulsion to 
outer-planet orbiter missions are assessed in this study, and guide- 
lines regarding future development are provided. A Jupiter orbiter 
mission was selected for evaluation because it is the earliest plane- 
tary mission that may require advanced chemical propulsion. Mis- 
sion and spacecraft characteristics that affect the selection and 
design of propulsion subsystems are presented. Alternative propul- 
sion subsystems, including space-storable bipropellant liquids, a 
solid/monopropellant vernier, and a hybrid, are compared on the 
basis of performance, reliability, and cost. The comparisons that 
assess performance, reliability, and cost independently do not yield 
a conclusive evaluation of each alternative propulsion subsystem's 
competitive position. This handicap was overcome by comparing 
the alternative propulsion subsystems with a cost-effectiveness 
model that combines the above three variables into a single param- 
eter. The results indicate that the hybrid and space-storable bipro- 
pellant mechanizations arc competitive. 
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DUMAS, L. N. 

D06 Temperature Control of the Mariner Mars 1971 Spacecraft 

L N. Dumas 

Technical Memorandum 33-515, January 1, 1972 

The Mariner Mars 1971 orbiter mission was a part of the ongoing 
program of unmanned planetary exploration. The spacecraft de- 
sign was based on that of Mariner Mars 1969, with changes as 
necessary to achieve mission objectives. The temperature control 
design for Mariner Mars 1971 is described in this memorandum, 
with emphasis on those areas in which significant changes were 
implemented. Developmental tasks are summarized and discussed, 
and initial flight data are presented. 


ECKSTEIN, R. E. 

E01 Program Listing for the Reliability Block Diagram Computation 
Program of JPL Technical Report 32-1543 

P. 0. Chelson and R. E. Eckstein 

Technical Memorandum 33-513, December 1, 1971 

For abstract, see Chelson, P. O. 


EISENMAN, A. 

E02 The Characterization of Facsimile Camera Systems for Lunar 
and Planetary Surface Exploration 

A. Eisenman 

JPL Quarterly Technical Review, Vol. 1, No. 3, pp. 1-16, 
October 1971 

For imagery of lunar and planetary surfaces from stationary or 
stopped remotely controlled vehicles, facsimile camera systems 
offer unique advantages: better geometric fidelity than television 
systems, more picture elements per frame, ability to take single 
frame panoramic pictures, a very high contrast ratio within a 
picture, a choice of one or more narrow-band spectral responses 
over a wide possible range, image transmission over a low bit-rate 
communications channel without storage, very low power, small 
size, low weight, and ruggedizability to meet space flight require- 
ments. Systems of this kind are under development for the 1975 
Viking Martian lander, and have been used by the Soviets on 
Lunas and Lunokhod. JPL has recently completed laboratory and 
field evaluation of an existing facsimile camera system. The results 
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of this test program are presented. The applicability of facsimile 
cameras to lunar and planetary rovers is demonstrated. 

ESTABROOK, F. B. 

E03 Hamiltonian Cosmology 

F. B. Estabrook and H. D. Wahlquist 

Phys. Lett., Vol. 35A, No. 6, pp. 453-454, July 12, 1971 

In this article, a simple polynomial Hamiltonian is given for type 
VIII and IX vacuum cosmologies. This Hamiltonian also describes 
solutions with time-like homogeneous 3-surfaces and thus suggests 
quantized versions that involve fluctuations of 3-space signature 
and topology. 


FALLER, J. E. 

F01 Preliminary Results of Laser Ranging to a Reflector on the 
Lunar Surface 

J. D. Mu : holland, C. 0. Alley (University of Maryland), 

P. L. Bender (Joint Institute for Laboratory Astrophysics), 
D. G. Currie (University of Maryland), 

R. H. Dicke (Princeton University), J. E. Faller (Wesleyan 
University), W. M. Kaula (University of California, Los 
Angeles), G. J. F. MacDonald (University of California, 
Santa Barbara), H. H. Plotkin (Goddard Space Flight 
Center), and D. T. Wilkinson (Princeton University) 

Space Research XI, pp. 97-104, Akademie-Verlag, Berlin, 
1971 

For abstract, see Mulholland, J. D. 


FANALE, F. P. 

F02 Lunar Fines and Terrestrial Rock Powders: Relative Surface 
Areas and Heats of Adsorption 

F. P. Fanale, D. B. Nash, and W. A. Cannon 

J. Geophys. Res., Vol. 76, No. 26, pp. 6459-6461, 

September 10, 1971 

Surface area measurements by Kr adsorption (BET method) indi- 
cate that Apollo 11 lunar fines and ground terrestrial mafic rock 
powders have similar effective surface areas that are a factor of 
10-100 higher than their geometrical or surficial surface areas. 
There is no evidence that a significant increase in surface rough- 
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ness for lunar fines has resulted from their peculiar history of 
exposure on the moon’s urface. Low apparent heats of adsorption 
for Kr adsorption on lunar fine material are consistent with the 
presence of glassy or glass-coated particles. 


FANSELOW, J. L. 

F03 DSN Progress Report for July-August 1971: The Goldstone 
Interferometer for Earth Physics 

J. L. Fanswlow, P. F. MacDoran, J. B. Thomas, 

J. G. Williams, C. Finnie, T. Sato, L. Skjerve (Philco-Ford 
Corporation), and D. Spitzmesser 

Technical Report 32-1526, Vol. V, pp. 45-57, 

October 15, 1971 

The first in a series of very long baseline interferometry feasibility 
demonstrations for applications to Earth Physics was conducted on 
January 29, 1971. In this demonstration two Goldstone tracking 
stations, the 26-m Echo station and the 64-m Mars station, were 
equipped with JPL hydrogen masei requency systems and oper- 
ated in electrically independent, although coordinated, observing 
modes. S-band (2.3 GHz) radio signals from 14 celestial radio 
sources were recorded at each station on digital magnetic tape. 
Later, these tapes were brought together for computer cross- 
correlation to produce an interferometric observable. Using the 
interferometer fringe frequency, a measurement of the two equa- 
torial baseline components was made. The discrepancy between 
these measurements and the inter-station geodetic survey was less 
than 30 cm in each of these two components. 


FINNIE, C. 

F04 DSN Progress Report for July-August 1971: The Goldstone 
Interferometer for Earth Physics 

J. L. Fanselow, P. F. MacDoran, J. B. Thomas, 

J. G. Williams, C. Finnie, T. Sato, L. Skjerve (Philco-Ford 
Corporation), and D. Spitzmesser 

Technical Report 32-1526, Vol. V, pp. 45-57, 

October 15, 1971 

For abstract, see Fanselow, J. L. 
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FLEISCHER, G. E. 


F05 Flexible Spacecraft Control System Design Procedures Utilizing 
Hybrid Coordinates 

P. W. Likins, E. L. Marsh, and G. E. Fleischer 
Technical Memorandum 33-493, September 15, 1971 
For abstract, see Likins, P. W. 


FLETCHER, B. C. 

F06 An Eclectic Integrated Circuit Reliability Model 

B. C. Fletcher 

Technical Memorandum 33-514, January 1, 1972 

With the advent of solid-state physics and the subsequent inte- 
grated device technology, the task of assessing device reliability 
has become increasingly more difficult. One problem is that the 
parameters of integrated circuits and large-scale-integration de- 
vices are readily accessible only from the periphery, resulting in a 
lack of data needed for precise model definitions ar.d constraints. 
Another problem is that existing modeling techniques can not, in 
general, be considered viable since there is no direct input for 
basic materials and process changes that tend to be inherent in 
technological change. Other problems witn existing n.odels origi- 
nate rrom the assumptions forming the basis of present model 
derivations. 

This memorandum presents an interim report on research relating 
to the development of an integrated circuit reliability model. 
Because of the complexities stated briefly above and described 
more fully in the memorandum, defects, modes, mechanisms 
(where possible), frequency of mechanism occurrence (where pos- 
sible), manufacturing/processing impact on failure, time depen- 
dency of mechanisms, screening influences (supplie and user), and 
device application had to be considered in the modeling rationale 
that is proposed. 


FLIEGEL, H. F. 

F07 DSN Progress Report for July-August 1971: A Worldwide 
Organization to Secure Earth-Related Parameters for Deep 
Space Missions 

H. F. Fliegel 
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Technical Report 32-1526, Vol. V, pp. 66-73, 
October 15, 1971 


A global express service to obtain timing and polar motion param- 
eters for deep space mission support has been organized through 
the Bureau International de l’Heure. The results are incorporated 
into a daily operation. This article outlines what the new sources 
of data are, what procedures are used to reduce the data, and what 
software is available to the user. 



j 

i 





GALULY, I. 

G01 On the Drag Experienced by a Spheroidal, Small Particle in 
Gravitational and Electrostatic Fields 

I. Gallily 

J. Colloid Interface Sci., Vol. 36, No. 3, pp. 325-339, 

July 1971 

The behavior of a suspended, rigid particle in a creeping motion 
depends generally on translation, rotation, and coupling tensors, 
which are geometric properties of the body. Previously, these 
second-rank tensors have been calculated only for a few simple 
forms. In this article, the drag on a spheroidal particle that moves 
in a creeping motion and performs simultaneously random rota- 
tions is calculated by a new method of averaging for cases of both 
gravitational and electrostatic fields. Appropriate expressions are 
written for the mobilities, which are computed subsequently for a 
range of spheroid sizes, axial ratios, and electrostatic field 
intensities. 


GEORGEVIC, R. M. 

G02 Mathematical Model of the Solar Radiation Force and Torques 
Acting on the Components of a Spacecraft 

R. M. Georgevic 

Technical Memorandum 33-494, October 1, 1971 

Solar radiation pressure exerts a mechanical force upon the surface 
of a spacecraft that intercepts the stream of photons coming from 
the Sun. For high-precision spacecraft attitude control and orbit 
determination, it is necessary to generate a precise mathematical 
model of the solar radiation force and the moment of that force; 
such a model must be more accurate than the currently used “flat 
surface” model, based on the radiation force on the effective cross- 
sectional area of the irradiated body. 


151 


‘i 


u 



l 

i- 


4 - 

i 



r 



In this memorandum, the general expressions for the solar radia- 
tion force and torques are derived in the vectorial form for any 
given reflecting surface, provided that the reflecting characteristics 
of the surface, as well as the value of the solar constant, are known. 
An appropriate choice of a spacecraft-fixed frame of reference 
leads to relatively simple expressions for the solar radiation forces 
and torques in terms of the functions of the Sun-spacecraft-Earth 
angle. 

The advantage of such a model over the standard flat-surface 
model is obvious, and it is very easy to find the expressions for the 
error of the standard model for any given reflecting surface. An- 
other advantage of the model is that it can be used for the effects 
of the air drag, solar wind, etc. 


GERBER, W. 

G03 Quenching of Solid-Propellant Rockets by Water Injection 

L. D. Strand and W. Gerber (Lockheed Aircraft Service Co.) 

J. Spacecraft Rockets, Vol. 8, No. 9, pp. 992-996, 
September 1971 

For abstract, see Strand, L. D. 


GOODWIN, P. S. 

G04 DSN Progress Report for July-August 1971: Helios Mission 
Support 

P. S. Goodwin 

Technical Report 32-1526, Vol. V, pp. 17-21, 

October 15, 1971 

Project Helios is the first NASA international deep space project, 
although there have been prior NASA international sub-orbital and 
Earth-orbiting cooperative space projects. Helios is a joint under- 
taking by the Federal Republic of West Germany and the United 
States of America, who divide the project responsibilities. Two 
unmanned scientific satellites are planned for heliocentric orbits: 
the first to be launched in mid-1974, and the second in late-1975. 
Prior volumes of this series describe the history and objectives of 
this program, the contemplated spacecraft configuration, and its 
telecommunications and telemetry systems. This article deals with 
the spacecraft’s command system, its requirements, and conceptual 
block diagrams. 
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G05 DSN Progress Report for September-October 1971: Helios 
Mission Support 

P. S. Goodwin 

Technical Report 32-1526, Vol. VI, pp. 25-32, 

December 15, 1971 

Project Helios, a joint endeavor between the United States and 
West Germany, will place two unmanned spacecraft into heliocen- 
tric orbits whose perihelion distance will come closer to the Sun 
than any previous or presently planned free-world deep space 
undertaking. The first spacecraft is expected to be launched in 
mid-1974 and the second in late-1975. Prior volumes of this series 
describe the history and objectives of this project, the contem- 
plated spacecraft configuration, and the spacecraft’s radio system. 
This article deals with the capabilities of the telecommunications 
link between the spacecraft and the DSN. 


GORDON, H. J. 

G06 Analysis of Mariner VII Pre-encounter Anomaly 

H. J. Gordon, S. K. Wong, and V. J. Ondrasik 

J. Spacecraft Rockets, Vol. 8, No. 9, pp. 931-937, 

September 1971 

The loss of signal at 127 h before the Mariner VII (Mariner 7) 
closest approach to Mars and the subsequent changes in the space- 
craft and its trajectory are described. Real-time orbi : i'ermina- 
tion activity and postencounter analysis are discussed. Great dif- 
ficulty was experienced in processing tracking data influenced by 
an unknown nongravitational force that could not be properly 
modeled. Results of simulations using known perturbations are 
presented. It is concluded that the battery case ruptured and 
vented into the interior of the spacecraft. This rupture caused 
corona discharges, and escaping gas produced the translational 
forces. 


GROTCH, S. L. 

G07 Computer Techniques for Identifying Low Resolution Mass 
Spectra 

S. L. Grotch 
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Anal. Chem., Vol. 43, No. 11, pp. 1362-1370, 

September 1971 

A number of computer programs have been developed for identi- 
fying low resolution mass spectra through search of an extensive 
library file. One set of programs used as a disagreement criterion 
the sum of the absolute values of the differences in peak height 
levels when peak height was encoded to 2, 8, or 10 4 levels at each 
nominal mass. Another program employed the maximum coinci- 
dence of the top N peaks. The programs were tested using 125 
unknowns and the .ecognition performances were compared. The 
maximum coincidence criterion was significantly poorer in recog- 
nition performance than the other techniques which increased in 
reliability as the number of levels increased. However, even the 
two-level system attained very high reliability. Since computer 
requirements and economic costs are likely to be minimal for this 
case, it might suffice for many applications. 


HADEK, V. 

HOI Electrical Properties of 7,7',8,8'-Tetracyanoquinodimethane 
Salts of lonene Polymers and Their Model Compounds 

V. Hadek, H. Noguchi, and A. Rembaum 

Macromolecules, Vol. 4, No. 4, pp. 494-499, 

July-August 1971 

Electrically conducting polymeric salts were prepared by the 
reaction of ionene polymers with LiTCNQ in the presence or 
absence of neutral TCNQ (tetracyanoquinodimethane). The spe- 
cific resistivity, the activation energy for conductivity, and the 
Seebeck coefficient were determined as a function of the number 
of CH 2 groups between positively charged nitrogens. The wide 
variations of electrical properties could not be correlated with the 
length of the polymethylene chain in the polymer. X-ray analysis 
of single crystals of model compounds revealed that the electrical 
properties depend mainly on crystal geometry. 


HANSON, R. J. 

H02 A Numerical Method for Solving Fredholm Integral Equations 
of the First Kind Using Singular Values 

R. J. Hanson 
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SIAM J. Numer. Anal., Vol. 8, No. 3, pp. 616-622, 

September 1971 

The integral equation in question is approximated by simple nu- 
merical quadrature formulas plus collocation. Each row of the 
resulting matrix equation for the unknown function values is 
weighted by the reciprocal of the standard deviation of the known 
function. A singular-value decomposition is used to obtain a solu- 
tion for the resulting linear system. By avoiding the use of the 
smallest singular values, an approximate solution is calculated that 
frequently solves the problem “close enough” and removes a great 
deal of the oscillations in the solution that are inherently present 
due to the ill-posed nature of the problem. Test cases and compu- 
tational results are presented. 

HARPER, L 

H03 DSN Progress Report for July-August 1971: Contributions to a 
Mathematical Theory of Complexity 

L. Harper (University of California, Riverside) and 
J. E. Savage (Brown University) 

Technical Report 32-1526, Vol. V, pp. 91-98, 

October 15, 1971 

This article is another in a series that attempts to define precisely 
and investigate the “computational complexity” of a general class 
of problems which includes many problems that occur in the DSN 
control center. Specific DSN control center questions that a theory 
of computational complexity will help define and answer include 
the problem of the optimum mix of core, disk, and drum storage in 
the control center, and the intelligent allocation of computational 
resources to flight projects of differing complexity in such a way 
that simultaneous real-time computing commitments can be made. 
This article shows that there exists such a theory which is capable 
of providing important information about the true complexity of 
several classes of non-trivial problems. 

H04 DSN Progress Report for July-August 1971: Some Results on 
the Matrix Multiplication Problem 

L. Harper (University of California, Riverside) and 
J. E. Savage (Brown University) 

Technical Report 32-1526, Vol. V, pp. 99-101, 

October 15, 1971 

Three results on the multiplication of two n X n matrices are 
presented. They contribute to our understanding of the complexity 
of matrix multiplication, and so of code decoding, tracking accu- 
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racy computation, antenna structural analysis and other DSN com- 
putational tasks. 


HARTLEY, R. B. 

H05 DSN Progress Report for July-August 1971: Apollo Mission 
Support 

R. B. Hartley 

Technical Report 32-1526, Vol. V, pp. 29-41, 

October 15, 1971 

The support provided by the DSN to the Manned Space Flight 
Network during the Apollo 15 mission is described. Support was 
provided from four 26-m (85-ft) DSN stations, the Goldstone 64-m 
(210-ft) antenna, the Ground Communications Facility, and the 
Space Flight Operations Facility. Premission and mission activities 
are discussed and a brief mission description is included. 


HARTOP, R. W. 

H06 DSN Progress Report for September-October 1971: 

Operational Time Sync Microwave Subsystem 

R. W. Hartop 

Technical Report 32-1526, Vol. VI, pp. 165-167, 

December 15, 1971 

The prototype Operational Time Sync Microwave Subsystem has 
been redesigned to accommodate the change in operational fre- 
quency and to increase its power-handling capability. The com- 
pleted system has been successfully tested to 115 kW continuous 
wave at 7150 MHz and is expected to operate at 150 kW continu- 
ous wave, the design goal. A description of the subsystem is 
presented. 


HEER, E. 

H07 Maximum Dynamic Response and Proof Testing 

J.N. Yang and E. Heer (NASA Headquarters) 

J. Eng. Mech. Div., Proc. ASCE, Vol. 97, No. EM4, 
pp. 1307-1313, August 1971 

For abstract, see Yang, J.-N. 
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HERMANN, A. M. 


H08 Electrical Conductivity of Elastomeric TCNQ Complexes Under 
Mechanical Stress 

A. M. Hermann (Tulane University), S. P. S. Yen, 

A. Rembaum, and R. F. Landel 

J. Polym. Sci. f Pt. B: Polym. Lett., Vol. 9, No. 8, 
pp. 627-633, August 1971 

A significant amount of research has been made in the area of 
macroscopically inhomogeneous filled rubbers: however, no study 
of a homogeneous rubbery semiconductor under mechanical stress 
has been reported previously. This article gives a preliminary 
report of a study on the electrical conductivity of rubbery poly- 
mers complexed with tetracyanoquinodimethane (TCNQ). Mea- 
surements on a rubbery polymer mixed intimately with LiTCNQ 
are also discussed. The carbon-filled rubbers previously investi- 
gated had decreased conductivities under load; the rubbery poly- 
meric-TCNQ polymers have either unchanged conductivities or, in 
some cases, increased conduct ities in the stretched state. This 
increase in conductivities ir *v be suggestive of stress-induced 
orientation. 


HINTZ, G. R. 

H09 A Viking Satellite Orbit Trim Strategy 

G. R. Hintz 

JPL Quarterly Technical Review, Vol. 1, No. 3, pp. 133-142, 
October 1971 

The Viking Project places a number of interesting and stringent 
requirements on the control of the satellite orbit to obtain recon- 
naissance and to prepare for lander release. To satisfy these re- 
quirements, different orbit trim maneuver strategies have been 
developed for two typical Viking missions. This article describes 
one of these strategies. In addition, a summary of recent numerical 
results is included to show that this strategy satisfies the mission 
requirements which have been identified. 


HOFMANN, A. H. 

H10 DSN Progress Report for July-August 1971: Pioneer F and G 
Mission Support Area 

A. H. Hofmann 
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Technical Report 32-1526, Vol. V, pp. 136-140, 
October 15, 1971 


With the advent of the third-generation computer systems and the 
increased complexity of interplanetary missions, the area required 
to support mission operations has exceeded available facilities 
within the SFOF. Consequently, the Pioneer F Jupiter flyby mis- 
sion support area will be located in the new Systems Development 
Laboratory, Building 264, at JPL. The mission support area, its 
relationship with the SFOF and other DSN facilities, and some 
anticipated operational problems caused by its remote location 
from the SFOF are described in this article. 


HOLMES, J. K. 

Hll A Note on Some Efficient Estimates of the Noise Variance for 
First-Order Reed-Muller Codes 

J. K. Holmes 

IEEE Trans. Inform. Theor., Vol. IT- 17, No. 5, pp. 628-630, 
September 1971 

The maximum-likelihood estimate based on order statistics for 
both a single-sample and an Af-sample estimate are derived for the 
first-order Reed-Muller code sent over a zero-mean white Gaussian 
noise channel. The estimation method is subject to an equipment- 
complexity constraint. The constraint imposed requires that the 
estimates be obtained from successive comparisons of the correla- 
tions in a serial system, without the use of arithmetic operations. 
The maximum-likelihood estimates obtained here are compared to 
other such estimates and also the unrestricted maximum-likelihood 
estimate. 


HUNTRESS, W. T., JR. 

H12 Ion Cyclotron Resonance Power Absorption: Collision 

Frequencies for C0 2 + , N 2 + , and H 3 + Ions in Their Parent 
Gases 

W. T. Huntress, Jr. 

J. Chem. Phys., Vol. 55, No. 5, pp. 2146-2155, 

September 1, 1971 

The complete solution for the equation of motion of an ion in the 
ICR cel) is shown to give results for the instantaneous power 
absorption in excellent agreement with experiment at all pressures. 
The instantaneous power absorption at resonance initially in- 
creases linearly with time, and at high pressures levels off to a 
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constant value at saturation where the energy gained by ions from 
the RF electric field is equal to the energy dissipated in collisions. 
An expression is also derived for the average kinetic energy of an 
ion at saturation in the steady-state limit. Pulsed ICR techniques 
are used to obtain the instantaneous power absorption curves for 
N 2 + , C0 2 + , and H 3 + ions in their parent gases as a function of 
pressure, from which are calculated the momentum transfer rate 
constants k, and the dependence of the rate constants on ion 
kinetic energy. At 293 °K, 


k(K 2 +) = k(C0 2 + ) 

= 0.67 X 10" 9 cm 3 -molecule _ 1 'Sec _1 
fc(H 3 + ) = 1.09 X 10" 9 cm 3 -molecule" 1 -sec -1 


For N 2 + ions in N 2 and C0 2 + ions in C0 2 , the rate constants are 
both significantly greater than that predicted by polarization the- 
ory and both rate constants increase significantly with increasing 
ion kinetic energy. This behavior is most likely a consequence of 
long-range resonant charge transfer outside the orbiting impact 
parameter. 


JACKSON, E. B. 

J01 DSN Progress Report for July-August 1971: DSN Research 
and Technology Support 

E. B. Jackson 

Technical Report 32-1526, Vol. V, pp. 120-121, 

October 15, 1971 

The major current activities of the Development Support Group, 
at both the Venus Deep Space Station and the Microwave Test 
Facility, are presented and accomplishments and progress for each 
are described. Activities include pulsar observations, planetary 
radar (including a general relativity experiment), 100-kW clock 
synchronization implementation, SDS 930 computer installation 
into the Mars Deep Space Station, precision antenna gain mea- 
surements, very long baseline interferometry, electromagnetic field 
survey at the Pioneer Deep Space Station, and new phase-lock 
receiver installation at the Microwave Test Facility. 

J02 DSN Progress Report for September-October 1971: DSN 
Research and Technology Support 

E. B. Jackson 
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Technical Report 32-1526, Vol. VI, pp. 147-148, 
December 15, 1971 


Accomplishments and progress during the major current activities 
of the Development Support Group at both the Venus Deep 
Space Station and the Microwave Test Facility are described, 
Activities include pulsar observations, planetary radar with very- 
high-resolution range measurement of the planet Mars, tricone 
support structure assembly, 100-kW clock synchronization imple- 
mentation, precision antenna gain measurement on the 26-m an- 
tenna, Block IV receiver/ exciter installation and testing, weak 
source observations of 13 sources, and observations of the planet 
Jupiter. 


JET PROPULSION LABORATORY 

J03 Proceedings of the Conference on Experimental Tests of 
Gravitation Theories, California Institute of Technology, 
Pasadena, California, November 11-13, 1970 

Jet Propulsion Laboratory 

Technical Memorandum 33-499, November 1, 1971 

The Conference on Experimental Tests of Gravitation Theories 
was the result of the feelings of a number of people that technol- 
ogy, particularly that spawned by previous space activities, had 
made it possible to think realistically in terms of a long-range 
cooperative effort in the testing of general relativity and other 
modern theories of gravity. The conference was sponsored by 
NASA, the European Space Research Organization, and JPL. 
These proceedings, edited by R. W. Davies, present 33 papers that 
explore various competitive philosophies on both the operational 
and the theoretical level. 

J04 Development and Testing of the Central Computer and 
Sequencer for the Mariner Mars 1971 Spacecraft 

Jet Propulsion Laboratory 

Technical Memorandum 33-501, October 15, 1971 

The central computer and sequencer subsystem (CC&S) has been 
an important part of the Mariner series of interplanetary space- 
craft since their inception. As with other spacecraft subsystems, 
the CC&S has increased in complexity and capability with each 
spacecraft project. This report describes the design, fabrication, 
and testing associated with the development of the Mariner Mars 
1971 CC&S subsystem. 
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J05 Development and Testing of the Pyrotechnic Subsystem for the 
Mariner Mars 1971 Spacecraft 

Jet Propulsion Laboratory 

Technical Memorandum 33-502, December 15, 1971 

The Mariner Mars 1971 pyrotechnic subsystem consists of the 
pyrotechnics switching assembly (PSA), the explosive squibs on the 
spacecraft, the spacecraft/Centaur release devices, and the pin- 
pullers. Numerous small changes from the Mariner Mars 1969 
design were effected in the PSA, but the basic capacitor-discharge 
approach remained the same. Two new squibs were developed; a 
completely new spacecraft/Centaur release device was designed, 
developed, and qualified; and minor design improvements were 
accomplished on the pinpullers for the Mariner Mars 1971 pro- 
gram. The basic support equipment hardware and approach were 
carried over from the Mariner Mars 1969 program. The purpose of 
this memorandum is to document the design, fabrication, and 
testing of the Mariner Mars 1971 pyrotechnic subsystem, with 
emphasis on changes from the Mariner Mars 1969 subsystem 
development. 

J06 Development and Testing of the S-Band Antenna Subsystem 
for the Mariner Mars 1971 Spacecraft 

Jet Propulsion Laboratory 

Technical Memorandum 33-503, November 1, 1971 

The Mariner Mars 1971 S-band antenna subsystem is used to 
transmit and receive S-band signals to and from the Deep Space 
Instrumentation Facility ground stations. The antenna subsystem 
consists of a low-gain antenna, a medium-gain antenna, a direc- 
tional coupler, a high-gain antenna, and all transmission lines 
required to interconnect the antennas to the spacecraft radio 
frequency subsystem. The low-gain antenna is used to transmit 
signals during cruise and receive signals throughout the mission. 
The medium-gain antenna is coupled to the low-gain antenna via 
the directional coupler and is used to transmit and receive signals 
during Mars orbit insertion. The high-gain antenna is used to 
transmit high-data-rate signals primarily during Mars orbit. 

J07 Development and Testing of the Television Instrument for the 
Mariner Mars 1971 Spacecraft 

Jet Propulsion Laboratory 

Technical Memorandum 33-505, November 1, 1971 

This memorandum describes the Mariner Mars 1971 television 
instrument, with emphasis on those aspects that are different from 
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the Mariner Mars 1969 television subsystem. The various modes of 
operation are described, and functional descriptions of the major 
elements in the system are summarized. An electronic description 
of the circuits that differ from those of Mariner Mars 1969 is also 
presented, along with a brief description of the calibration and test 
sequences. 

J08 A Reduced Star Catalog Containing 537 Named Stars 

Jet Propulsion Laboratory 

Technical Memorandum 33-507, November 15, 1971 

This document is the first of a series presenting data, compiled by 
JPL, that are to be included in the JPL Astronautical Star Catalog. 
Positional and color magnitude data for the 537 named stars, which 
are to be included in the catalog, are given. A brief translation of 
the star names and the source language of the names are also 
presented. 

J09 Development and Testing of the Data Automation Subsystem 
for the Mariner Mars 1971 Spacecraft 

Jet Propulsion Laboratory 

Technical Memorandum 33-508, November 15, 1971 

The data automation subsystem serves as the data handling subsys- 
tem for the science payload. It provide:- basic timing and com- 
mands to control the instrument’s modes and sequences. The 
instrument data are formatted 'or transmission to Earlh in real 
time via the telemetry channel or stored on the tape recorder for 
playback in non-real time. The design, fabrication, and testing of 
the Mariner Mars 1971 data automation subsystem are described in 
this memorandum. Major changes are discussed relative to the 
Mariner Mars 1969 subsystem design. 

J10 Development and Testing of the Infrared Interferometer 
Spectrometer for the Mariner Mars 1971 Spacecraft 

Jet Propulsion Laboratory 

Technical Memorandum 33-509, December 1, 1971 

The primary purpose of the Mariner Mars 1971 infrared spectros- 
copy experiment is to determine values of atmospheric parameters 
and information pertine: t to the solid surface of Mars. Both atmo- 
spheric and surface parameters can be derived by interpretation of 
thermal emission spectra, which will be observed from the orbiter 
as a function of time and location on the pianet. 


162 



The Mariner interferometer, which was designed, fabricated, and 
tested by Texas Instruments, Inc., under the direction of the NASA 
Goddard Space Flight Center, is similar to the interferometer 
flown on the meteorological research satellites Nimbus 3 and 4 in 
1969 and 1970. In this memorandum, the unique features of the 
Mariner interferometer are emphasized, with particular regard to 
the cesium iodide beamsplitter. The extended spectral range made 
possible with the cesium iodide beamsplitter is significant to the 
water vapor investigation. The major changes required to accom- 
modate the increased performance and to adapt to the Mariner 
spacecraft environment are described. 

Jll Mariner Mars 1971 Science Operational Support Equipment 
Final Report 

Jet Propulsion Laboratory 

Technical Memorandum 33-511, November 15, 1971 

The Mariner Mars 1971 science operational support equipment 
(SOSE) was developed to support the checkout of the proof test 
model and flight spacecraft. This memorandum discusses the test 
objectives for the SOSE and how these objectives were imple- 
mented. Attention is focused on the computer portion, since incor- 
poration of a computer in ground checkout equipment represents 
a major departure from the support equipment concepts previ- 
ously used at JPL. The major hardware elements and the SOSE 
operational performance during spacecraft testing are described. 

J12 u.^LG— A Program for Test Pattern Generation in 

Combinational Logical Circuits 

Jet Propulsion Laboratory 

Technical Memorandum c 3-510, November 15, 1971 

This memorandum is primarily a user’s manual for a computer 
program DALG that generates test patterns for detecting faults in 
combinational logic circuits containing up to 200 logical gates. The 
gates may be of logical types AND, OR, NAND, NOR, NOT, or 
Exclusive OR. The faults may be any one gate or input stuck at a 
fixed value (0 or 1). 

In addition to test pattern generation, DALG will also determine 
whether or not the given test pattern .vill detect given faults in a 
circuit. Sample problems are given along with input data sheets 
and piinted output to illustrate the capabilities of the program. 
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JOHNSON, D. 


J13 DSN Progress Report for July-August 1971: Combinational 
Complexity Measures as a Function of Fan-Out 

D. Johnson, J. E. Savage (Brown University), and 
L. R. Welch (University of Southern California) 

Technical Report 32-1526, Vol. V, pp. 79-81, 

October 15, 1971 

If C s (fi, . . . t fj) is the fan-out s combinational complexity of the 
functions fi,f 2 , . . . ,fi with respect to straight-line algorithms (or 
combinational machines) of fan-in r, then it is shown that 



1 





where N is the number of variables on which /j, . . . ,f L depend 
and d = C s (7), where I is the identity function in one variable. 
Thus, a well-designed combinational machine or algorithm will not 
have a fan-out which is more than several times its fan-in. 


JOHNSTON, A. R. 

J14 Dispersion of Electro-Optic Effect in BaTi0 3 
A. R. Johnston 

J. Appl. Phys., Vol. 42, No. 9, pp. 3501-3507, August 1971 

The dispersion of both the quadratic and linear electro-optic 
effects in flux-grown crystals of barium titanate (BaTi0 3 ) has been 
measured polarimetrically between 0.4 and 1.0 ft yielding (gn-gi 2 ) 
and 



A strong dispersion was found: (gn— gi 2 ) increases to double its 
long-wave-length limit at 0.4 ft. Both unclamped (low frequency) 
and clamped measurements were made. Birefringence and the 
principal indices were also determined. A two-oscillator Sellmeier 
model, in which one oscillator frequency is polarization depen- 
dent, was shown to represent closely all the data. The unclamped 
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polarization potential, which specifies the magnitude of the oscilla- 
tor frequency shift, was 3.9 eV m 4 C -2 , while its clamped value was 
one-third as large. Ultraviolet reflectivity calculated from the 
model assuming reasonable values of empirical damping agrees 
with the observations of Cardona and Gahwiller. 


JORDAN, J. F. 

J15 Guidance and Navigation for Solar Electric Interplanetary 
Missions 

K. H. Rourke and J. F. Jordan 

J. Spacecraft Rockets, Vol. 8, No. 9, pp. 920-926, 
September 1971 

For abstract, see Rourke, K. H. 


KALFAYAN, S. H. 

KOI Long-Term Aging of Elastomers: Chemical Stress Relaxation of 
Fluorosilicone Rubber and Other Studies 

S. H. Kalfayan, A. A. Mazzeo, and R. H. Silver 

JPL Quarterly Technical Review, Vol. 1, No. 3, pp. 38-47, 
October 1971 

Elastomers have varied aerospace applications, including: propel- 
lant binders, bladder materials for liquid propellant expulsion 
systems, and fuel tank sealants, particularly for high-speed aircraft. 
Predicting the long-term behavior of these materials is of para- 
mount importance. A comprehensive molecular theory for me- 
chanical properties has been developed at JPL. It has only been 
tested experimentally in cases where chemical degradation 
processes are excluded. Hence, a study is underway to ascertain 
the nature, extent, and rate of chemical changes that take place in 
some elastomers of interest. The results can then be incorporated 
in the theoretical framework. This article reports progress on the 
investigations of chemical changes that may take place in the 
fluorosilicone elastomer, LS 420, which is regarded as a fuel and 
high-temperature-resistant rubber. The kinetic analysis of the 
chemical stress relaxation and gel permeation chromatography 
studies comprise the major portion of this article. 


165 



KATOW, M. S. 


K02 DSN Progress Report for September-October 1971: S- and 
X-Band Feed System 

M. S. Katcw 

Technical Report 32-1526, Vol. VI, pp. 139-141, 

December 15, 1971 

The proposed S- and X-band feed system provides for simultaneous 
RF propagation from the 64-m antenna for both S- and X-band 
signals along the same boresight direction. The hardware for the 
tri-cone system consists of an ellipsoid reflector over the S-band 
horn and a dichroic reflector plate over the X-band cone. The 
ellipsoid reflector focuses the S-band signal in front of the dichroic 
plate. The dichroic plate is capable of transmitting an X-band 
signal through it and reflecting S-band. The dichroic plate, 
mounted in a position about 60 deg to the centerline of the X-band 
signal, then reflects the S-band signal coincident to the X-band 
signal. Preliminary hardware mounting schemes are outlined with 
probable operational requirements. 

K03 DSN Progress Report for September-October 1971: Overseas 
64-m RMS Program for SDS 920 

D. McCarty and M. S. Katow 

Technical Report 32 1526, Vol. VI, pp. 158-164, 

December 15, 1971 

For abstract, see McCarty, D. 


KAULA, W. M. 

K04 Preliminary Results of Laser Ranging to a Reflector on the 
Lunar Surface 

J. D. Mulholland, C. 0. Alley (University of Maryland), 

P. L. Bender (Joint Institute for Laboratory Astrophysics), 
D. G. Currie (University of Maryland), 

R. H. Dicke (Princeton University), J. E. Faller (Wesleyan 
University), W. M. Kaula (University of California, Los 
Angeles), G. J. F. MacDonald (University of California, 
Santa Barbara), H. H. Plotkin (Goddard Space Flight 
Center), and D. T. Wilkinson (Princeton University) 

Space Research XI, pp. 97-104, Akademie-Verlag, Berlin, 
1971 

For abstract, see Mulholland, J. D. 
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KHATIB, A. R. 
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K05 Dynamic Upper Atmospheric Force Model on Stabilized 

Vehicles for a High-Precision Trajectory Computer Program 

A. R. Khatib 

JPL Quarterly Technical Review, Vol. 1, No. 3, pp. 125-132, 
October 1971 

This article summarizes the results of research carried out at JPL 
for the design and implementation of dynamic upper atmosphere 
and lift and drag models into the advanced Double-Precision Tra- 
jectory Program (DPTRAJ). The upper atmosphere model draws 
heavily on the behavior of the Earth’s upper atmosphere which 
exhibits cyclic as well as irregular variations in density profile, 
temperature, pressure, and composition in unison with solar activi- 
ties as deduced from the more recent land-based and satellite 
observations. 

The lift and drag model is designed specifically for inertially 
stabilized vehicles of the Mariner class, with possible extension to 
gravity gradient stabilized vehicles of the GEOS class. The model 
considers operation in the free molecular flow regimes with large 
Knudsen numbers. The vehicle is considered a composite structure 
with basic components having well-defined shapes, each with its 
own surface characteristics in terms of temperature, reflectivity, 
and accommodation of free stream molecules. The model lakes 
into account both the calculation of precise aerodynamic force 
coefficients in terms of expansion of modified Bessel functions in 
speed ratios and angle of attack, and approximate force coefficients 
when the speed ratios approach infinity. Other considerations 
include specular and diffused reflectivity, shielding, and shadow 
effects. 


KINDER, W. J. 

K06 DSN Progress Report for September-October 1971: DSN 
Telemetry System Tests 

W I Kinder 2nd R. S. Basset 

Technical Report 32-1526, Vol. VI, pp. 10-12, 

December 15, 1971 

The overall DSN System test plan, as edited for the DSN Multi- 
mission Telemetry System, is briefly described. Specific results, 
with delivered Mariner Mars 1971 telemetry software, are pre- 
sented in relation to this test plan. Recommendations are included 
for future system demonstration, specifically as to system docu- 
mentation and training. 
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KRON, M. 


K07 DSN Progress Report for July-August 1971: Load Distribution 
on the Surface of Paraboloidal Reflector Antennas 

M. Kron 

Technical Report 32-1526, Vo! V, pp. 122-128, 

October 15, 1971 

Wind pressure coefficients have been measured using wind tunnel 
models of parabolic reflectors. The application of this data and its 
conversion to useful form for structural deflection analysis within 
the “NASTRAN” S f ructural Analysis Computer Program and ulti- 
mately Root Mean Square (RMS) program is described in this 
article. 


LACY, G. H. 

L01 Farthest South Soil Microbial and Ecological Investigations 

R. E. Cameron, G. H. Lacy, F. A. Morelli, and 
J. B. Marsh (University of California, Davis) 

Antarctic J. U.S., Vol. VI, No. 4, pp. 105-106, 
July-August 1971 

For abstract, see Cameron, R. E. 


LAESER, R. P. 

L02 DSN Progress Report for July-August 1971: Mariner Mars 
1971 Mission Support 

R. P. Laeser 

Technical Report 32-1526, Vol. V, pp. 22-23, 

October 15, 1971 

While continuing to provide the tracking and data acquisition 
function for Mariner 9 on its journey to Mars, the DSN is plaiming 
and practicing for the orbital operations. The final steps of imple- 
mentation and the test plans for orb*.ul operations are outlined in 
this article. 

L03 DSN Progress Report for September-October 1971: Mariner 
Mars 1971 Mission Support 

R. P. Laeser 
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Technical Report 32-1526, Vol. VI, pp. 33-36, 

December 15, 1971 

This article completes and updates the description of the planned 
DSN configuration for support of the Mariner 9 orbit insertion and 
orbital operations. Specifically covered are the S-band occultation 
experiment data handling, the planetary ranging configu-ition, and 
the simulation configuration. 
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LANDEL, R. F. 

L04 Electrical Conductivity of Elastomeric TCNQ Complexes Under 
Mechanical Stress 

A. M. Hermann (Tulane University), S. P. S. Yen, 

A. Rembaum, and R. F. Landel 

J. Polym. Sci., Pt. B: Polym. Lett., Vol. 9, No. 8, 
pp. 627-633, August 1971 

For abstract, see Hermann, A. M. 

L05 Properties of a Highly Crosslinked Elastomer 
R. F. Landel 

Polymer Networks : Structural and Mechanical Properties, pp. 
219-243, Plenum Publishing Corporation, New York, 1971 

With the objective of enhancing understanding of the behavior of 
elastomers, the properties of a highly crosslinked polyurethane 
elastomer, based on ricinoleic acid (castor oil) as the backbone, 
have been measured. The rubber, termed Galcit I, has been pro- 
posed as a standardized, highly birefringent, nearly elastic rubber 
for use in testing new apparatus on a material with known proper- 
ties and in checking or calibrating existing instruments. This article 
discusses the preparation and characterization of Galcit I. 


LAYLAND, J. W. 

L06 DSN Progress Report for Septem ber-October 1971: An 
Optimum Buffer Management Strategy for Sequential 
Decoding 

J. W. Layland 

Technical Report 32-1526, Vol. VI, pp. 106-111, 

December 15, 1971 

Sequential decoding has been found to be an efficient means of 
communicating at low undetected error rates from deep space 
probes, but a failure mechanism known as erasure or computa- 
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tional overflow remains a significant problem. The erasure of a 
block occurs when the decoder has not finished decoding that 
block at the time that it must be output. 

The erasure rate can be unacceptably high, even when the de- 
coder is spending over half of its time idly awaiting incoming data. 
By drawing upon analogies in computer time-sharing, this article 
develops a buffer management strategy that reduces the decoder 
idle time to a negligible level and, therefore, improves the erasure 
probability of a sequential decoder. For a decoder with speed 
advantage of 10 and buffer size of 10 blocks, operating at an 
erasure rate of 10~ 2 , use of the new buffer management strategy 
reduces the erasure rate to less than 10 -4 . 
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LAYMAN, W. E. 

L07 Development and Testing of the Beryllium Propulsion Support 
Structure for the Mariner Mars 1971 Spacecraft 

J. H. Stevens and W. E. Layman 

Technical Memorandum 33-517, January 1, 1972 

For abstract, see Stevens, J. H. 

LIKINS, P. W. 




LOS Finite Element Appendage Equations for Hybrid Coordinate 
Dynamic Analysis 

P. W. Likins 

Technical Report 32-1525, October 15, 1971 

The increasingly common practice of idealizing a spacecraft as a 
collection of interconnected rigid bodies to some of which are 
attached linearly elastic flexible appendages leads to equations of 
motion expressed in terms of a combination of: (1) discrete coordi- 
nates describing the arbitrary rotational motions of the rigid 
bodies, and ( 9, 1) distributed or modal coordinates describing the 
small, time-varying deformations of the appendages; such a formu- 
lation is said to employ a hybrid system of coordinates. In this 
report, the existing literature is extended to provide hybrid coordi- 
nate equations of motion for a finite element model of a flexible 
appendage attached to a rigid base undergoing unrestricted mo- 
tions, and some of the advantages of the finite element approach 
are noted. Transformations to the modal coordinates appropriate 
for the general case and various special cases are provided. 
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L09 Flexible Spacecraft Control System Design Procedures Utilizing 
Hybrid Coordinates 

P. W. Likins, E. L. Marsh, and G. E. Fleischer 

Technical Memorandum 33-493, September 15, 1971 

Procedures for the practical implementation of the hybrid coordi- 
nate methods of dynamic analysis of flexible spacecraft in applica- 
tion to vehicles of realistic complexity are briefly documented, 
with supporting examples. 


LINNES, K. W. 

L10 DSN Progress Report for July-August 1971: Radio Science 
Support 

K. W. Linnes 

Technical Report 32-1526, Vol. V, pp. 42-44, 

October 15, 1971 

Since 1967 radio scientists have used the DSN 26- and 64-m 
antenna stations to investigate pulsars, to study the effects of solar 
corona on radio signals, and to observe radio emissions of X-ray 
sources. More recently, very long baseline interferometry (VLBI) 
techniques have been used for high resolution studies of quasars. 
During the reporting period, several proposals were received for 
extension of VLBI observations which had reported the startling 
expansion of quasar 3C279. 

Lll DSN Progress Report for September-October 1971: Radio 
Science Support 

K. W. Linnes 

Technical Report 32-1526, Vol. VI, pp. 43-45, 

December 15, 1971 

Since 1967, radio scientists have used the DSN 26- and 64-m 
antenna stations to investigate pulsars, to study the effects of solar 
corona on radio signals, and to observe radio emissions from X-ray 
sources. More recently, very-long-baseline interferometry (VLBI) 
techniques have been used for high-resolrbon studies of quasars. 
During the reporting period, VLBI observations were made in 
support of investigations of quasars and the application of VLBI 
techniques to Earth physics problems. Support was also provided 
for preliminary investigation of the mapping of spiral galaxies. 
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L.ZZ r<?»i . v ogress Report for July-August 1971: Sequential Tests 
*or Exponential Distributions 

G. Lorden (California Institute of Technology) 

Technical Report 32-1526, Vol. V, pp. 82-90, 

October 15, 1971 

The problem is to test whether the frequency of random events 
(e.g., DSIF equipment failures) is at a nominally prescribed value. 
When the actual frequency is higher, a determination of this fact is 
to be made as quickly as possible. A test based on sequential 
maximum likelihood ratio methods is developed and approxima- 
tions of its performance characteristics are derived. Results of 
Monte Carlo sampling demonstrate that these approximations are 
accurate and that high statistical efficiency is attained over a broad 
range of possible higher frequencies. Some applications to reliabil- 
ity and inventory policies for the DSIF are indicated. 


MacDONALD, G. J. F. 

M01 Preliminary Results of Laser Ranging to a Reflector on the 
Lunar Surface 

J. D. Mulholland, C. 0. Alley (University of Maryland), 

P. L. Bender (Joint Institute for Laboratory Astrophysics), 
D. G. Currie (University of Maryland), 

R. H. Dicke (Princeton University), J. E. Faller (Wesleyan 
University), W. M. Kaula (University of California, Los 
Angeles), G. J. F. MacDonald (University of California, 
Santa Barbara), H. H. Plotkin (Goddard Space Flight 
Center), and D. T. Wilkinson (Princeton University) 

Space Research XI, pp. 97-104, Akademie-Verlag, Berlin, 
1971 

For abstract, see Mulholland, J. D. 


MacDORAN, P. F. 

M02 DSN Progress Report for July-August 1971: The Goldstone 
Interferometer for Earth Physics 

J. L. Fanselow, P. F. MacDoran, J. B. Thomas, 

J. G. Wihiams, C. Finnic, T. Sato, L. Skjerve (Philco Ford 
Corporation), and D. Spitzmesser 
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Technical Report 32-1526, Vol. V, pp. 45-57, 
October 15, 1971 

For abstract, see Fanselow, J. L. 
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McCarty, d. 

MC3 DSN Progress Report for September-October 1971: Overseas 
64-m RMS Program for SDS 920 

D. McCarty and M. S. Katow 

Technical Report 32 1526, Vol. VI, pp. 158-164, 

December 15, 1971 

With the completion of the 64-m antennas overseas and their 
performance testing, an important test measurem nt required is 
the reflector distortion RMS from gravity loading. In order to 
provide the paraboloid best-fitting capability for the available SDS 
920 computers at overseas sites, the RMS program was modified to 
suit the typewriter-and-two-tape-units input-output capabilities of 
the computers. The program computes the RMS after paraboloid 
best fit from field angle readings using typed inputs. The constant 
data, such as coordinates of targets, are supplied in a data tape 
with the binary program supplied in a second tape. 


McCLURE, J. P. 

M04 DSN Progress Report for September-October 1971: 

GCF 50- kbps Wideband Data Error Rate Test 

J. P. McClure 

Technical Report 32-1526, Vol. VI, pp. 149-157, 

December 15, 1971 

During June 1971, a 7-day wideband data error test was conducted 
between the Space Flight Operations Facility communications 
terminal and the NASCOM Madrid Switch Center. The test, 
which was run at 50 kbits/s, was conducted to determine both 
long-term and short-term error data for a wideband circuit compa- 
rable to those expected to be used to support the Mariner Venus- 
Mercury 1973 Project. Long-term end-to-end error rates of 6 X 
10" 5 or better were measured in both directions. The hourly and 
5-min error distributions indicate that the errors are grouped into 
bursts (as expected). Most of the time the error rate is substantially 
less than average. 
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McGINNESS, H. 


M05 DSN Progress Report for $ep*ember-October 1971- Movement 
of the Antenna Instrument Tower at DSS 14 

H. McGinness 

Technical Report 32-1526, Vol. VI, pp. 142-146, 

December 15, 1971 

The motions of the top of the instrument tower and its surround- 
ing windshield have been measured. A relationship between a 
static horizontal displacement and an angular displacement of the 
tower have been established through the use of optical apparatus. 
Displacements during excitation of the windshield have been de- 
termined by the use of accelerometers. The nature of the coupling 
between the windshield and tower is discussed. The conclusion 
reached is that the coupling is primarily an acoustical one. 


McLYMAN, W. T. 

M06 Magnetic Materials Selection for Static Invert*' and Converter 
Transformers 

W. T. McLyman 

Technical Memorandum 33-498, November 1, 1971 

A program was conducted to study magnetic materials for use in 
spacecraft transformers used in static inverters, converters, and 
transformer-rectifier supplies. A comparative investigation of dif- 
ferent magnetic alloys best suited for high-frequency and high- 
efficiency applications was conducted, together with an investiga- 
tion of each alloy’s inherent characteristics. The trade names ana 
magnetic alloys of the materials evaluated were: Orthono!: 50% 
Ni, 50% Fe; Sq. Permalloy: 79% Ni, 17% Fe, 4% Mo; 48 alloy: 48% 
Ni, 52% Fe; Supermalloy: 78% Ni, 17% Fe, 5% Mo; and Magnesil: 
3% Si, 97% Fe 

One characteristic of magnetic materials that is detrimental in 
transformer design is the residual flux density, which can be addi- 
tive on turn-on and cause the transformer to saturate. Investigation 
of this problem led to the design of a transformer with a very low 
residual flux. Tests were performed to determine the dc and ac 
magnetic properties at 2400 Hz using square-wave excitation. 
These tests were performed on uncut cores, which were then cut 
for comparison of the gapped and ungapped magnetic properties. 
When the data of many tra’ formers in many configurations were 
compiled, the optimum transformer was found to be that with the 
lowest residual flux and a small amount of air gap in *he magnetic 
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material. The data obtained from these tests are described, and the 
potential uses for the materials are discussed. 

MACIE, T. W. 

M07 Solar Electric Propulsion System Technology 

T. D. Masek and T. VV. Made 

Technical Memorandum 33-510, November 15, 1971 

For abstract, see Masek, T. D. 


MARGOLIS, J. S. 

M08 Recomputation of the Absorption Strengths of the Methane 
3v 3 J-Manifolds at 9050 cm 1 

J. T. Bergstralh (McDonald Observatory) and J. S. Margolis 

J. Quant Spectrosc. Radiat. Transfer, Vol. 11, No. 8, 
pp. 1285-1287, August 1971 

For abstiact, see Bergstralh, J. T. 


MARSH, E. L. 

M09 Flexible Spacecraft Control System Design Procedures Utilizing 
Hybrid Coordinates 

P. W. Likins, E. L. Marsh, and G. E. Fleischer 
Technical Memorandum 33-493, September 15, 1971 
For abstract, see Likins, P. W. 


MARSH, J. B. 

M10 Farthest South Soil Microbial and Ecological Investigations 

R. E. Cameron, G. H. Lacy, F. A. Morelli, and 
J. B. Marsh (University of California, Davis) 

Antarctic J. U.S., Vol. VI, No. 4, pp. 105-106, 
July-August 1971 


For abstract, see Cameron, R. E. 



MARTONCHIK, J. V. 


Mil Absorption by Venus in the 3-4-Micron Region 

R. Beer, R. H. Norton, and J. V. Martonchik (University of 
Texas) 

Astrophys. J., Vol. 168, No. 3, Pt. 2, pp. L121-L124, 
September 15, 19?1 

For abstract, see Beer, R. 

M'12 Astronomical Infrared Spectroscopy With a Connes-Type 
Interferometer: II. Mars, 2500-3500 cm- 1 

R. Beer, R. H. Norton, and J. V. Martonchik (University of 
Texas) 

Icarus: Int. J. Sol. Sys., Vol. 15, No. 1, pp. 1-10, 

August 1971 

For abstract, see Beer, R. 


MASEK, T. D. 

M13 Solar Electric Propulsion System Technology 

T. D. Masek and T. W. Macie 

Technical Memorandum 33-510, November 15, 1971 

As the number of possible applications for primary solar-powered 
electric propulsion grows, the burden of demonstrating this tech- 
nology grows in proportion. The solar electric propulsion system 
technology (SEPST) program at JPL is focusing on such a demon- 
stration. This memorandum reports the progress of the present 
JPL hardware program (SEPST III) and discusses certain propul- 
sion-system-spacecraft interaction problems being investigated. 

The basic solar electric propulsion system concept and elements 
are reviewed. Hardware is discussed only briefly, relying on de- 
tailed fabrication or assembly descriptions reported elsewhere. 
Emphasis is placed on recent performance data, which are pre- 
sented to show the relationship between spacecraft requirements 
and present technology. 


MATHUR, F. P. 

M14 Reliability Estimation Procedures and CARE: The Computer- 
Aided Reliability Estimation Program 

F. P. Mathur 
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JPL Quarterly Technical Review, Vol. 1, No. 3, pp. 17-26, 
October 1971 


Ultrareliable fault-tolerant onboard digital systems for spacecraft 
intended for long mission life exploration of the outer planets are 
under development. The design of systems involving self-repair 
and fault-tolerance leads to the companion problem of quantifying 
and evaluating the survival probability of the system for the 
mission under consideration and the constraints imposed upon the 
system. Methods have been developed to (1) model self-repair and 
fault-tolerant organizations; (2) compute survival probability, mean 
life, and many other reliability predictive functions with respect to 
various systems and mission parameters; (3) perform sensitivity 
analysis of the system with respect to mission parameters; and (4) 
quantitatively compare competitive fault-tolerant systems— various 
measures of comparison are offered. To automate the procedures 
of reliability mathematical modeling and evaluation, the CARE 
(computer-aided reliability estimation) program was developed. 
CARE is an interactive program residing on the UNIVAC 1108 
system, which makes the above calculations and facilitates report 
preparation by providing output in tabular form and graphical 
2-dimensional plots and 3-dimensional projections. The reliability 
estimation of fault-tolerant organization by means of the CARE 
program is described in this article. 


MAZZEO, A. A. 

M15 Long-Term Aging of Elastomers: Chemical Stress Relaxation of 
Fluorosilicone Rubber and Other Studies 

S. H. Kalfayan, A. A. Mazzeo, and R. H. Silver 

JPL Quarterly Technical Review, Vol. 1, No. 3, pp. 38-47, 
October 1971 

For abstract, see Kalfayan, S. H. 


MILLER, L. F. 

M16 DSN Progress Report for July-August 1971: Comparison of 
Faraday Rotation Measurements of the Ionosphere 

L. F. Miller and 6. D. Mulhall 

Technical Report 32-1526, Vol. V, pp. 58-65, 

October 15, 1971 

An evaluation of the mapping techniques employed to provide 
ionospheric charged particle calibration for post-flight analysis and 
for Mariner Mars 1971 tracking system analytical calibration oper- 
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ations was performed. Comparisons based on Faraday rotation data 
from geostationary satellites were made between various satellites 
as recorded at the Venus Deep Space Station and by Stanford 
Center for radar astronomy. 


M0RELL1, F. A. 

M17 Farthest South Soil Microbial and Ecological Investigations 

R. E. Cameron, G. H. Lacy, F. A. Morelli, and 
J. B. Marsh (University of California, Davis) 

Antarctic J. U.S., Vol. VI, No. 4, pp. 105-106, 
July-August 1971 

For abstract, see Cameron, R. E. 


MOYNIHAN, P. I. 

X 

. £ M18 TOPS Attitude Propulsion Subsystem Technology 

P. I. Moynihan 

JPL Quarterly Technical Review, Vol. 1, No. 3, pp. 48-56, 
October 1971 

This article summarizes the JPL Thermoelectric Outer-Planet 
Spacecraft (TOPS) Attitude Propulsion Subsystem effort through 
the end of fiscal year 1971. It includes the tradeoff rationale that 
went into the selection of anhydrous hydrazine as the propellant, 
followed by a brief description of three types of 0.445-N (100- 
mlbf) thrusters that were purchased for in-house evaluation. A 
discussion is also included of the 0.2224-N (50-mlbf) JPL-developed 
thrusters and their integration with a portable, completely en- 
closed, propulsion module that was designed and developed to 
support the TOPS single-axis attitude control tests in the JPL 
Celestarium. The article concludes with a synopsis of further work 
which will be accomplished prior to the onset of an outer-planet 
mission. 
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MUDGWAY, D. J. 

M19 DSN Progress Report for September-October 1971: Viking 
Mission Support 


D. J. Mudgway 
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Technical Report 32-1526, Vol. VI, pp. 37-42, 

December 15, 1971 

The Tracking and Data System Functional Specification and the 
NASA Support Plan have been completed for the Viking Project. 
A complex scheduling problem, created by the Viking Project 
request for mission design verification tests in late 1974 and early 
1975, has been solved by reworking early agreements on responsi- 
bility for software development. 

The Viking Project poses the problem of simultaneous multiple RF 
links to the DSN for the first time. Consequently, it has been 
necessary to introduce multiple-link requirements into the current 
DSN techniques for single-link RF compatibility testing. The effect 
of these new requirements, with particular reference to Viking, is 
discussed in this article. 
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MULHALL, B. D. 

M20 DSN Progress Report for July-August 1971: Comparison of 
Faraday Rotation Measurements of the Ionosphere 

L. F. Miller and B. D. Mulhall 

Technical Report 32-1526, Vol. V, pp. 58-65, 

October 15, 1971 

For abstract, see Miller, L. F. 


MULHOLLAND, J. D. 

M21 Preliminary Results of Laser Ranging to a Reflector on the 
Lunar Surface 

J. D. Mulholland, C. 0. Alley (University of Maryland), 

P. L. Bender (Joint Institute for Laboratory Astrophysics), 

D. G. Currie (University of Maryland), 

R. H. Dicke (Princeton University), J. E. Faller (Wesleyan 
University), W. M. Kaula (University of California, Los 
Angeles), G. J. F. MacDonald (University of California, 

Santa Barbara), H. H. Plotkin (Goddard Space Flight 
Center), and D. T. Wilkinson (Princeton University) 

Space Research XI, pp. 97-104, Akademie-Verlag, Berlin, 
1971 

The first Lunar Ranging Experiment (LURE) retroreflector array 
was placed on the lunar surface during the Apollo 11 mission. Prior 
to this event, a special high-precision lunar ephemeris (designated 
LE 16) was developed by means of a composite numeric /analytic 
process, so as to provide more accurate predictions for use at the 
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telescope. First returns were observed at the Lick Observatory on 
August 1, 1969 and at the McDonald Observatory shortly thereaf- 
ter. The observing program is to continue for several years. Pre- 
liminary use of these data consists of their comparison with the LE 
16 epbemeris, preparatory to a differential correction of the lunar 
elements. Present indications are that the ephemeris must undergo 
order-of-magnitude improvements before the full power of the 
laser data can be utilized. Other parameters in the predictive 
process are also capable of being corrected; an improvement in 
the coordinates of the Lick Observatory 120-in. telescope is al- 
ready indicated. 


NASH, D. B. 

N01 Lunar Fines and Terrestrial Rock Powders: Relative Surface 
Areas and Heats of Adsorption 

F. P. Fanale, D. B. Nash, and W. A. Cannon 

J. Geophys. Res., Vol. 76, No 26, pp. 6459-6461, 
September 10, 1971 

For abstract, see Fanale, F. P. 


NEUGEBAUER, M. 

N02 Computation of Solar Wind Parameters From the OGO-5 
Plasma Spectrometer Data Using Hermite Polynomials 

M. Neugebauer 

Technical Memorandum 33-519, December 15, 1971 

This memorandum presents the method used to calculate the 
velocity, temperature, and density of the solar wind plasma from 
spectra obtained by attitude-stabilized plasma detectors on the 
Earth satellite Orbiting Geophysical Observatory 5 (OGO-5). The 
method, which uses expansions in terms of Hermite polynomials, is 
very inexpensive to implement on an electronic computer com- 
pared to the least-squares and other iterative methods often used 
for similar problems in the past. 


NOGUCHI, H. 

N03 Electrical Properties of 7,7',8,8'-Tetracyanoquinodimethane 
Salts of lonene Polymers and Their Model Compounds 

V. Hadek, H. Noguchi, and A. Rembaum 
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Macromolecules, Vol. 4, No. 4, pp. 494-499, 
July-August 1971 

For abstract, see Hadek, V. 


NORTON, R. H. 

N04 Absorption by Venus in the 3-4-Micron Region 

R. Beer, R. H. Norton, and J. V. Martonchik (University of 
Texas) 

Astrophys. J., Vol. 168, No. 3, Pt. 2, pp. L121-L124, 
September 15, 1971 

For abstract, see Beer, R. 

N05 Astronomical Infrared Spectroscopy With a Connes-Type 
Interferometer: II. Mars, 2500-3500 cm 1 

R. Beer, R. H. Norton, and J. V. Martonchik (University of 
Texas) 

Icarus: Int. J. Sol. Sys., Vol. 15, No. 1, pp. 1-10, 

August 1971 

For abstract, see Beer, R. 


ODLYZKO, A. M. 

001 DSN Progress Report for September-October 1971: Data 
Storage and Data Compression 

A. M. Odlyzko 

Technical Report 32-1526, Vol. VI, pp. 112-117, 

December 15, 1971 

In this article, a sharp upper bound on the best possible data rate 
achievable is computed as a function of data storage capability in 
certain very general situations. The result shows that a dramatic 
increase in rate can be caused by a small increase in storage 
capability. 


ONDRASIK, V. J. 

002 DSN Progress Report for September-October 1971: The 

Repetition of Seasonal Variations in the Tropospheric Zenith 
Range Effect 

K. L. Thuleen and V. J. Ondrasik 
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Technical Report 32-1526, Vol. VI, pp. 83-98, 
December 15, 1971 

For abstract, see Thuleen, K. L. 

003 Analysis of Mariner VII Pre-encounter Anomaly 

H. J. Gordon, S. K. Wong, and V. J. Ondrasik 

J. Spacecraft Rockets, Vol. 8, No. 9, pp. 931-937, 
September 1971 

For abstract, see Gordon, H. J. 
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OTOSHI, T. Y. 

004 DSN Progress Report for July-August 1971: Antenna Noise 
Temperature Contributions Due to Ohmic and Leakage Losses 
of the DS$ 14 64-m Antenna Reflector Surface 

T. Y. Otoshi 

Technical Report 32-1526, Vol. V, pp. 115-119, 

October 15, 1971 

This article presents approximate formulas useful for computing 
antenna noise temperature contributions due to ohmic and leakage 
losses of a parabolic antenna reflector surface. The total noise 
temperature contributions due to ohmic and leakage losses for the 
DSS 14 64-m antenna were calculated to be 0.1, 0.3, and 0.6 K at 
2.295, 8.448, and 15.3 GHz, respectively. 

005 DSN Progress Report for September-October 1971: Further 
Studies of Microwave Transmission Through Perforated Flat 
Plates 

T. Y. Otoshi and K. Woo 

Technical Report 32-1526, Vol. VI, pp. 125-129, 

December 15, 1971 

This article presents approximate formulas useful for predicting 
transmission loss characteristics of a circular hole array in a metal- 
lic flat plate having finite thickness. The formulas apply to perpen- 
dicular and parallel polarizations of an obliquely incident plane 
wave. The approximate formulas are experimentally verified by 
free space measurements made on a sample of the mesh material 
used on the 64-m antenna at the Mars Deep Space Station. 



5 


s 



PHEN, R. L. 


P01 Application of Hybrid Propulsion Systems to Planetary 
Missions 

J. P. Don and R. L. Phen 

Technical Memorandum 33-483, November 1, 1971 
For abstract, see Don. J. P. 
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PLOTKIN, H. H. 

P02 Preliminary Results of Laser Ranging to a Reflector on the 
Lunar Surface 

J. D. Mulholland, C. 0. Alley (University of Maryland), 

P. L. Bender (Joint Institute for Laboratory Astrophysics), 
D. G. Currie (University of Maryland), 

R. H. Dicke (Princeton University), J.- E. Faller (Wesleyan 
University), W. M. Kaula (University of California, Los 
Angeles), G. J. F. MacDonald (University of California, 
Santa Barbara), H. H. Plotkin (Goddard Space Flight 
Center), and D. T. Wilkinson (Princeton University) 

Space Research XI, pp. 97-104, Akademie-Verlag, Berlin, 
1971 

For abstract, see Mulholland, J. D. 


POULSON, P. 

P03 computer Program for the Automated Attendance Accounting 
System 

P. Poulson and C. Rasmusson 

JPL Quarterly Technical Review, Vol. 1, No. 3, pp. 27-32, 
October 1971 

The Automated Attendance Accounting System (AAAS) was devel- 
oped under the auspices of the Space Technology Applications 
Office at JPL. The task is basically the adaptation of a small digital 
computer, coupled with specially developed pushbutton terminals 
located in school classrooms and offices for the purpose of taking 
daily attendance, maintaining complete attendance records, and 
producing partial and summary reports. Especially intended for 
high schools, the system will relieve both teachers and office 
personnel from the time-consuming and dreary task of recording 
and analyzing the myriad classroom attendance data collected 
throughout the semester. In addition, since many school district 
budgets are related to student attendance, the increase in account- 
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ing accuracy is expected to augment district income. A major 
component of this system is the real-time AAAS software system, 
which is described in this article. 


PRICE, T. W. 

P04 Long-Duration Firings of a Mariner Mars 1969 Catalytic 
Reactor 

T. W. Price 

JPL Quarterly Technical Review, Vol. 1, No. 3, pp. 57-66, 
October 1971 

Two long-duration tests were conducted with a surplus Mariner 
Mars 1969 monopropellant hydrazine reactor in an attempt to 
induce the “washout” phenomenon. The Mariner Mars 1969 reac- 
tor was chosen because it has a long development history and thus 
is well characterized. No “washout” occurred during either of the 
two 1000-s tests, although slow transients were observed in the 
reactor operation during what were nominally steady-state condi- 
tions. The 2000 s of operating time represents nearly an order of 
magnitude increase over the rated life of the engine. 


QUINN, R. 

QC' DSN Progress Report for July-August 1971: Tracking and 
Data Acquisition Elements Research: Low Noise Receivers: 
Microwave Maser Development 

R. C. Clauss and R. Quinn 

Technical Report 32-1526, Vol. V, pp. 102-108, 

October 15, 1971 

For abstract, see Clauss, R. C. 


RASMUSSQN, C. 

R01 Computer Program for the Automated Attendance Ac** junting 
System 

P. Poulson and C. Rasmusson 

JPL Quarterly Technical Review, Vol. 1, No. 3, pp. 27-32, 
October 1971 

For abstract, see Poulson, P. 
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REID, M. S. 


R02 DSN Progress Report for Septem ber- October 1971: Improved 
RF Calibration Techniques: System Operating noise 
Temperature Calibrations 

M. S. Reid 

Technical Report 32-1526, Vol. VI, pp. 130-138, 

December 15, 1971 

The system operating noise temperatures of the S-band research 
operational cone at the Venus Deep Space Station and the polar- 
ization diversity S-band cone at the Mars Deep Space Station are 
reported for the period June 1, 1971 through September 30, 1971. 
In addition, the performance of the multi-frequency X- and K-band 
(MXK) cone on the ground at the Venus Deep Space Station is 
reported for X-band operation, as well as for X-band operation on 
the 64-m antenna at the Mars Deep Space Station for the same 
period. Also presented are system operating noise temperature 
calibrations of the K-band system in the following configuration*: 
before installation in the MXK cone (approximately 23 K), installed 
in the cone, with the cone on the ground (approximately 25 K), 
and with the cone installed on the 64-m antenna at the Mars Deep 
Space Station (approximately 29 K). 


REILLY, H. F., JR. 

R03 DSN Progress Report for September-October 1971: Microwave 
Maser Development 

R. C. Clauss and H. F. Reilly, Jr. 

Technical Report 32-1526, Vol. VI, pp. 118-122, 

December 15, 1971 

For abstract, see Clauss, R. C. 


REMBAUM, A. 

R04 Onset of Superconductivity in Sodium and Potassium 
Intercalated Molybdenum Disulphide 

R. B. Somoano and A. Rembaum 

JPL Quarterly Technical Review, Vol. 1, rio. 3, pp. 33-37, 
October 1971 

For abstract, see Somoano, R. B. 
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R05 Electrical Conductivity of Elastomeric TCNQ Complexes Under 
Mechanical Stress 

A. M. Hermann (Tulane University), S. P. S. Yen, 

A. Rembaum, and R. F. Landel 

J. Polym. Sci.. Pt. B: Polym. Lett., Vol. 9, No. 8, 
pp. 627-633, August 1971 

For abstract, see Hermann, A. M. 

R06 Electrical Properties of 7,7',8,8'>Tetracyancquinodimethane 
Salts of lonene Polymers and Their Model Compounds 

V. Hadek, H. Noguchi, and A. Rembaum 

Macromolecules, Vol. 4, No. 4, pp. 494-499, 

July-August 1971 

For abstract, see Hadek, V. 

RC7 Superconductivity in Intercalated Molybdenum Disulfide 

R. B. Somoano and A. Rembaum 

Phys. Rev. Lett., Vol. 27, No. 7, pp. 402-404, 

August 16, 1971 

For abstract, see Somoano, R. B. 


RENZETTI, N. A. 

R08 DSN Progress Report for July-August 1971: DSN Functions 
and Facilities 

N. A. Renzetti 

Technical Report 32-1526, Vol. V, pp. 1-3, 

October 15, 1971 

The DSN, established by the NASA Office of Tracking and Data 
Acquisition and under the system management and technical direc- 
tion of JPL, is designed for two-way communications with unman- 
ned spacecraft traveling approximately 16,000 km (10,000 mi) from 
earth to planetary distances. The objectives, functions, and organi- 
zation of the DSN are summarized, and its three facilities— the 
Deep Space Instrumentation Facility, the Ground Communications 
Facility, and the Space Flight Operations Facility— are described. 

R09 DSN Progress Report for September-October 1971: DSN 
Functions and Facilities 

N. A. Renzetti 
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Technical Report 32-1526, Vol. VI, pp. 1-4, 

December 15, 1971 

The DSN, established by the NASA Office of Tracking and Data 
Acquisition and under the system management and technical direc- 
tion of JPL, is designed for two-way communications with unman- 
ned spacecraft traveling approximately 16,000 km (10,000 mi) from 
Earth to planetary distances. The objectives, functions, and organi- 
zation of the DSN are summarized, and its three facilities— the 
Deep Space Instrumentation Facility, the Ground Communications 
Facility, and the Space Flight Operations Facility— are described. 


REY, R. D. 

R10 DSN Progress Report for September-October 1971: Angle 
Tracking Analysis and Test Development 

R. D. Rey 

Technical Report 32-1526, Vol. VI, pp. 170-187, 

December 15, 1971 

The angle tracking systems are currently being analyzed, and tests 
are being developed to measure their performance. This article 
presents the progress made on the analysis and testing of the 
standard 26-m-diam antenna station automatic angle tracking sys- 
tem. The model is discussed, and certain important system con- 
stants are developed. Simulation runs of the model were per- 
formed, and comparisons are made with preliminary tests per- 
formed at the Echo Deep Space Station. The article also outlines 
the design of the test and software processing. 


RHO, J. H. 

Rll Grating Anomalies in Porphyrin Spectra 

J. H. Rho 

Geochim. Cosmochim. Acta, Vol. 35, No. 7, pp. 743-747, 

July 1971 

Experimental evidence is presented to show that porphyrin spec- 
tra are modified by grating anomalies to such an extent that 
interpretation is difficult. The anomalies appear in a definite wave- 
lei gth region which depends on the optical properties of individ- 
ual grating monochromators, and are associated with polarization 
effects. The presence of light-scattering material in the samples of 
any fluorescent compounds may induce the formation of grating 
anomalies in their fluorescence spectra. The anomalies may easily 
be mistaken for porphyrin peaks, especially in samples which 
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contain light-scattering material such as those from geological 
sources. 


R5NDERLE, E. A. 

R12 DSN Progress Report for July-August 1971: A Comparison of 
Cowell’s Method and a Variation-of-Parameters Method for 
the Computation of Precision Satellite Orbits 

S. S. Dallas and E. A. Rinderle 

Technical Report 32-1526, Vol. V, pp. 74-78, 

October 15, 1971 

For abstract, see Dallas, S. S. 


RODRIGUEZ, C. F. 

R13 Simulation of Mariner Mars 1971 Spacecraft 

N. E. Ausman, Jr., N. K. Simon, and C. F, Rodriguez 

JPL Quarterly Technical Review, Vol. 1, No. 3, pp. 67-78, 
October 1971 

For abstract, see Ausman, N. E., Jr. 

ROUKLOVE, P. 

R14 Parametric Testing of an Externally Configured Thermionic 
Converter 

K. Shimada and P. Rouklove 

Technical Memorandum 33-504, November 1, 1971 

For abstract, see Shimada, K. 


ROURKE, K. H. 

R15 Guidance and Navigation for Solar Electric Interplanetary 
Missions 

K. H. Rourke and J. F. Jordan 

J. Spacecraft Rockets, Vol. 8, No. 9, pp. 920-926, 

September 1971 

This article makes a practical analysis of closed-loop guidance of a 
solar electrically thrusted interplanetary spacecraft. A first-order 
guidance algorithm is used that allows easy interpretation yet 
rigorous treatment of trajectory terminal constraints. Interplane- 
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tary orbit determination is assumed to be provided by Earth-based 
radio (doppler) tracking. Guidance and orbit determination are 
evaluated with respect to given trajectories and then combined to 
yield information on the encounter accuracies and guidance cost- 
in terms of additional power consumption, control deviations, and 
payload penalty— expected from given injection and thruster per- 
formance dispersions. The results are presented in terms of two 
representative cases: a Jupiter flyby/orbiter mission and a comet 
rendezvous mission in the late 1970s. 


SATO, T. 

$01 DSN Progress Report for July-August 1971: The Goldstone 
Interferometer for Earth Physics 

J. L. Fanselow, P. F. MacDoran, J. B. Thomas, 

J. G. Williams, C. Finnie, T. Sato, L. Skjerve (Philco-Ford 
Corporation), and D. Spitzmesser 

Technical Report 32-1526, Vol. V, pp. 45-57, 

October 15, 1971 

For abstract, see Fanselow, J. L. 


SAVAGE, J. E. 

502 DSN Progress Report for July-August 1971: Combinational 
Complexity Measures as a Function of Fan-Out 

D. Johnson, J. E. Savage (Brown University), and 
L. R. Welch (University of Southern California) 

Technical Report 32-1526, Vol. V, pp. 79-81, 

October 15, 1971 

For abstract, see Johnson, D. 

503 DSN Progress Report for July-August 1971: Contributions to a 
Mathematical Theory of Complexity 

L. Harper (University of California, Riverside) and 
J. E. Savage (Brown University) 

Technical Report 32-1526, Vol. V, pp. 91-98, 

October 15, 1971 

For abstract, see Harper, L. 
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S04 DSN Progress Report for July-August 1971: Some Results on 
the Matrix Multiplication Problem 

L. Harper (University of California, Riverside) and 
J. E. Savage (Brown University) 

Technical Report 32-1526, Vol. V, pp. 99-101, 

October 15, 1971 

For abstract, see Harper, L. 






SAWYER, C. D. 

S05 DEXTER— A One-Dimensional Code for Calculating Thermionic 

Performance of Long Converters 

C. D. Sawyer 

Technical Report 32-1545, November 15, 1971 

This report describes a versatile code for computing the coupled 
thermionic electric-thermal performance of long thermionic con- 
verters in which the temperature and voltage variations cannot be 
neglected. The code is capable of accounting for a variety of 
external electrical connection schemes, coolant flow paths, and 
converter failures by partial shorting. Sample problem solutions 
are included, along with a user’s manual. 


SCHORN, R. A. J. 


S06 Comments on "The Venus Spectrum: New Evidence for Ice” 

R. A. J. Schorn and L. D. G. Young 

Icarus: Int. J. Sol. Sys. f Vol. 15. No. 1, pp. 103-109, 
August 1971 


In a recent article in Icarus , Plummer has attempted to show that 
high-altitude infrared spectra obtained by Kuiper’s group exhibit 
evidence for ice-crystal clouds on Venus. He also asserts that these 
data are consistent with ground-based spectroscopic work which 
indicates ~100 ppt fim of water vapor “above the clouds” of 
Venus. Such an interpretation of the high-altitude spectra is not 
required by the data and, in fact, raises more problems than it 
answers. The bulk of ground-based observations indicate an H 2 0 
abundance of much less than 100 ppt (im . 


In this article, no conclusions are made about the composition of 
the clouds of Venus. It is merely pointed out that the airborne and 
ground-based spectra offer no convincing evidence for an ice-cloud 
composition. 
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SHIMADA, K. 


507 Parametric Testing of an Externally Configured Thermionic 
Converter 

K. Shimada and P. Rouklove 

Technical Memorandum 33-504, November 1, 1971 

As described in this memorandum, a 25.4-cm-long externally con- 
figured converter was performance-tested at JPL by electrically 
heating the emitter to simulate reactor thermal power input. The 
measured maximum output power was limited by the maximum 
input power available from the electric RF induction heater. With 
maximum heater input power, the converter electric output was 
178 W (1.95 W/cm 2 ) at an emitter temperature of 1946 K. This 
electric output power was smaller than expected. The power 
output during acceptance testing at the contractor’s site was 5 
W/cm 2 at a 2000-K emitter temperature. The converter withstood 
46 controlled shutdowns and 13 abrupt shutdowns without shorting 
and without loss of ces ^m. 

A reactor-core-length (25.4-cm-long) cylindrical thermionic con- 
verter could be used as a full-length thermionic fuel element. 
Obtaining high output power and maintaining the emitter-to- 
collector gap without shorting are of major importance to the 
feasibility of a 25.4-cm-long reactor fuel element. The emitter of 
the converter is located externally to the collector to increase the 
fuel-volume fraction and to allow redundant collector cooling in a 
reactor configuration. 

508 Measurements of Plasma Parameters in a Simulated 
Thermionic Converter 

K. Shimada 

JPL Quarterly Technical Review, Vol. 1, No. 3, pp. 97-109, 
October 1971 

Cesium-filled thermionic energy converters are being considered as 
candidate electrical energy sources in future spacecraft requiring 
tens to hundreds of kilowatts of electrical power. The high operat- 
ing temperatures necessary for a large specific power and high 
efficiency inevitably impose stringent constraints on the converter 
fabrication to achieve the desired reliability of the power system. 
The converter physics for reducing operating temperatures and 
cesium plasma losses are being studied to achieve high reliability 
without sacrificing the power performance of the converters. Vari- 
ous cesium parameters which affect the converter performance 
are: (1) electron temperatures, (2) plasma ion densities, and (3) 
electric potential profiles. These were investigated using a Lang- 
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muir probe in a simulated converter. The parameters were mea- 
sured in different cesium discharge modes. 
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SHIRLEY, D. L. 

S09 An Approach to Automated Drug Identification 

D. L. Shirley 

J. Forensic Sci., Vol. 16, No. 3, pp. 359-375, July 1971 

The process of selecting one approach to the automation of drug 
screening tests is outlined. A systems analysis was performed which 
determined the requirements of forensic science laboratories for an 
automated drug identification system, and these requirements were 
compared with a set of model systems. A system that matched the 
most representative set of requirements was refined by the selec- 
tion of more detailed approaches. An approach selected for initial 
hardware development and critical technology testing involved: (1) 
wet chemical sample preparation; (2) gas chromatographic separa- 
tion, presumptive identification, and quantitation of drugs; (3) 
infrared spectrophotometric identification of drugs; and (4) com- 
puter control and data analysis. 


SIEGMETH, A. J. 

510 DSN Progress Report for July-August 1971: Pioneer Mission 
Support 

A. J. Siegmeth 

Technical Report 32-1526, Vol. V, pp. 4-16, 

October 15, 1971 

A description of the planned configuration and data flow methodol- 
ogy of the Mark III Deep Space Network System is given. This 
system will support the Pioneer F and G missions and the succes- 
sive projects of the NASA mission set of the 1970s. Block diagrams 
graphically illustrate the planned functions of the DSN Telemetry, 
Tracking, and Command Systems including their capabilities of 
being compatible with the forthcoming project requirements. The 
basic interfaces between subsystems of the three DSN facilities are 
defined. 

511 DSN Progress Report for September-October 1971: Pioneer 
Mission Support 

A. J. Siegmeth 
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Technical Report 32-1526, Vol. VI, pp. 13-24, 
December 15, 1971 


The DSN plans to use the Mark III system configuration for 
tracking and data acquisition support of the Pioneer F and G 
missions. As a continuation of the description of the network 
systems, the configurations of the Simulation, Monitoring, and 
Operations Control Systems are given. Block diagrams show the 
planned functions, data flow methodology, and interfaces between 
the subsystems and the three DSN facilities. 


SILVER, R. H. 

S12 Long-Term Aging of Elastomers: Chemical Stress Relaxation of 
Fluorosilicone Rubber and Other Studies 

S. H. Kalfayan, A. A. Mazzeo, and R. H. Silver 

JPL Quarterly Technical Review, Vol. 1, No. 3, pp. 38-47, 
October 1971 

For abstract, see Kalfayan, S. H. 


SIMON, N. K. 

S13 Simulation of Mariner Mars 1971 Spacecraft 

N. E. Ausman, Jr., N. K. Simon, and C. F. Rodriguez 

JPL Quarterly Technical Review, Vol. 1, No. 3, pp. 67-78, 
October 1971 

For abstract, see Ausman, N. E., Jr. 


SKJERVE, L. 

S14 DSN Progress Report for July-August 1971: The Goldstone 
Interferometer for Earth Physics 

J. L. Fanselow, P. F. MacDoran, J. B. Thomas, 

J. G. Williams, C. Finnie, T. Sato, L. Skjerve (Philco-Ford 
Corporation), and D. Spitzmesser 

Technical Report 32-1526, Vol. V, pp. 45-57, 

October 15, 1971 

For abstract, see Fanselow, J. L. 
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SOMOANO, R. B. 


515 Onset of Superconductivity in Sodium and Potassium 
Intercalated Molybdenum Disulphide 

R. B. Somoano and A. Rembaum 

JPL Quarterly Technical Review, Vol. 1, No. 3, pp. 33-37, 
October 1971 

Molybdenum disulfide in the form of natural crystals or powder has 
been intercalated at -65 to -70°C with sodium and potassium 
using the liquid ammonia technique. All intercalated samples were 
found to show a superconducting transition. A plot of the percent 
of diamagnetic throw versus temperature indicates the possible 
existence of two phases in the potassium intercalated molybdenum 
disulfide. The onset of superconductivity in potassium and sodium 
intercalated molybdenite powder was found to be =6.2 and =4.5 
K, respectively. The observed superconductivity is believed to be 
due to an increase in electron density as a result of intercalation. 

516 Superconductivity in Intercalated Molybdenum Disulfide 

R. B. Somoano and A. Rembaum 

Phys. Rev. Lett, Vol. 27, No. 7, pp. 402-404, 

August 16, 1971 

In an attempt to investigate superconductors which are two- 
dimensional in nature, molybdenum disulfide has been intercalated 
with sodium and potassium. Both natural crystals of molybdenite 
and synthetic crystals were used, and measurements made on the 
intercalated products indicate a superconducting transition tem- 
perature of ~1.3°K for sodium and ~4.5°K for potassium. 


SPITZMESSER, D. 

S17 DSN Progress Report for July-August 1971: The Goldstone 
Interferometer for Earth Physics 

J. L. Fanselow, P. F. MacDoran, J. B. Thomas, 

J. G. Williams, C. Finnie, T. Sato, L. Skjerve (Philco-Ford 
Corporation), and D Spitzmesser 

Technical Report 32-1526, Vol. V, pp. 45-57, 

October 15, 1971 

For abstract, see Fanselow, J. L. 
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STEVENS, J. H. 


S18 Development and Testing of the Beryllium Propulsion Support 
Structure for the Mariner Mars 1971 Spacecraft 

J. H. Stevens and W. E. Layman 

Technical Memorandum 33-517, January 1, 1972 

In November 1971, the Mariner 9 spacecraft will be injected into 
Martian orbit by a 574-kg (1265-lbm) propulsion system. Support 
for that system is provided by an 8.9-kg (19.5-lbm) truss assembly 
consisting of beryllium tubes adhesively bonded to magnesium end 
fittings. Beryllium was selected for the tubular struts in the truss 
because of its exceptionally high stiffness-to-weight ratio. Adhesive 
bonding, rather than riveting, was utilized to join the struts to the 
end fittings because of the low toughness (high notch sensitivity) of 
beryllium. Magnesium, used in the end fittings, resulted in a 50% 
weight savings over aluminum, since geometric factors in the 
fitting design resulted in low stress areas where magnesium’s lower 
density was a benefit. This memorandum describes the design, 
development, testing, and fabrication procedures and problems 
associated with the development of the Mariner 9 propulsion 
support truss structure. 


STRAND, L. D. 

S19 Low-Acceleration Solid-Propellant Rocket Ignition Study 
L. D. Strand 

Technical Memorandum 33-506, December 1, .1971 

A study was conducted to develop a solid-propellant rocket igniter 
system that would build up thrust at a controlled rate of less than 
0.2 g/s. The system consisted of a long-burning, regressive-burning, 
control led-flow igniter and an inhibited progressive-burning sur- 
face in the main rocket motor. The igniter performed the dual role 
of igniting, under vacuum back-pressure and low L* conditions, the 
nonrestricted portion of the propellant and providing the mass 
addition necessary to sustain combustion until the propellant burn- 
ing area had increased sufficiently to provide a stable motor- 
chamber pressure. Two series of tests were conducted with exist- 
ing small test motor hardware to: (1) demonstrate the feasibility of 
the concept, (2) determine the important parameters governing 
the system, and (3) obtain design guidelines for future scaled-up 
motor tests. A quasi-steady-state mass balance tor the ignition 
system was written and programmed for use as a motor design 
tool. 
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S20 Quenching of Solid-Propellant Rockets by Water Injection 

L. D. Strand and W. Gerber (Lockheed Aircraft Service Co.) 

J. Spacecraft Rockets , Vol. 8, No. 9, pp. 992-996, 

September 1971 

In other programs, command termination of solid-propellant mo- 
tors by water quench has been demonstrated in motors having 
propellant weights up to 8600 kg. However, most attempts to 
correlate extinguishment data have given questionable results. 
Three mechanisms have been proposed: (1) rapid cooling of gases 
causes a dP/dt sufficient for extinction; the water then wets the 
surface, cooling it to prevent reignition; (2) cooling of the gases 
below some threshold value lowers heat feedback to the propellant 
below that necessary for self-supported combustion; and (3) a 
water film covers the entire surface of the propellant, cooling it 
below the temperature required for burning. 

The study described in this article was initiated to improve under- 
standing of the quench mechanism and to determine the optimum 
method of water injection. The slab-burning window motor used is 
capable of accepting several different types of water injectors: 
head-end injectors, multiple injectors impinging normal to the 
propellant surface, and sheet injectors which lay a thin sheet of 
water onto the propellant surface. The spray form varied from a 
fine mist to a solid stream. High-speed movies were taken during 
water injection and quenching. 


THOMAS, J. B. 

T01 DSN Progress Report for July-August 1971: The Goldstone 
Interferometer for Earth Physics 

J. L. Fanselow, P. F. MacDoran, J. B. Thomas, 

J. G. Williams, C. Finnie, T. Sato, L. Skjerve (Philco*Ford 
Corporation), and D. Spitzmesser 

Technical Report 32-1526, Vol. V, pp. 45-57, 

October 15, 1971 

For abstract, see Fanselow, J. L. 


THORMAN, H. C. 

T02 DSN Progress Report for Sep - nber-October 1971: DSN 
Simulation System 
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H. C. Thorman 



Technical Report 32-1526, Vol. VI, pp. 5-9, 
December 15, 1971 


The DSN Simulation System provides real-time insertion of simu- 
lated Tracking, Telemetry, Command, Monitor, and Operations 
Control Systems data into the DSN. Data flows originating from 
the system are used extensively in testing and training activities to 
prepare the DSN and its users for coverage of planned missions. 
This article describes the upgrading of the Simulation System that 
was accomplished to provide support of DSN development, test- 
ing, and training activities in 1970 and 1971. 


THULEEN, K. L. 

T03 DSN Progress Report for September-October 1971: The 

Repetition of Seasonal Variations in the Tropospheric Zenith 
Range Effect 

K. L. Thuleen and V. J. Ondrasik 

Technical Report 32-1526, Vol. VI, pp. 83-98. 

December 15, 1971 

Using radiosonde balloon data taken from sites close to the DSN 
tracking stations, the tropospheric zenith range effect A p z has 
been computed throughout 1967 and 1968. The behavior of A p z 
has definite seasonal trends that are similar in both years. With the 
modification of the tropospheric model, which is used to calibrate 
radio tracking data to include these seasonal trends, the naviga- 
tional errors, produced by inaccuracies in representing the zenith 
range effect, may possibly be reduced by as much as 40%. 


TRAJMAR, S. 

T04 Differential and Integral Cross Sections for the Electron-Impact 
Excitation of the a'A g and b J 2 g + States of 0 2 

S. Trajmar, D. C. Cartwright (The Aerospace Corporation), 
and W. Williams 

Phys. Rev., Pt. A: Gen. Phys., Vol. 4, No. 4, pp. 1482-1492, 
October 1971 

Electron-impact energy-loss spectra of 0 2 have been analyzed for 
incident electron energies from 4 to 45 eV, scattering angles from 
10 to 90 deg, and energy losses from 0 to 5 eV. The inelastic 
processes observed were the excitation of the a ] Ag and b l 2, g J 
electronic states and vibrational exc s tation in some cases to v" “ 
13. The excitation cross sections at each energy were made abso- 
lute by normalizing the s>un of the integral cross sections (all 
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inelastic, ionization, and elastic) to measured electron-0 2 total 
cross sections. The differential cross sections for the and 

& 1 S g + states show nearly isotropic behavior, as expected for opti- 
cally spin-forbidden transitions. The elastic differential cross sec- 
tions are strongly forward peaked at higher energies, but become 
only slightly forward peaked at the lower energies. The integral 
cross sections for the excitation of the a 1 A g and ^ 1 2 g + states reach 
their maxima near 7 eV and are more than an order of magnitude 
larger than previous estimates. The integral elastic cross section 
reaches its maximum at around 10 eV. 


TUSTIN, D. G. 

T05 DSN Progress Report for September-October 1971: Network 
Allocation Schedules 

D. G. Tustin 

Technical Report 32-1526, Vol. VI, pp. 168-169, 

December 15, 1971 

This article reviews the reasons and needs for the Network Alloca- 
tion Schedules and briefly describes the make-up of these sched- 
ules and how they are used. The major emphasis is placed on the 
implementation of these schedules, including new special-purpose 
software. This software makes use of an existing file management 
program and IBM 360 utility programs. 


VON ROOS, 0. H. 

V01 DSN Progress Report for September-October 1971: Analysis of 
Dual-Frequency Calibration for Spacecraft VLBI 

0. H. von Roos 

Technical Report 32-1526, Vol. VI, pp. 46-56, 

December 15, 1971 

In this article, a feasibility study is undertaken and a detailed 
analysis is made of a wide-band very-long-baseline interferometer 
(VLBI) for the purpose of ranging and tracking a spacecraft. The 
system works on two frequencies (S- and X-band). By a new 
correlation technique, it is shown that it is possible to extract 
information on the total electron content with a rather high degree 
of accuracy, an accuracy certainly impossible to achieve with 
tracking modes currently in use. The total electron content and its 
time variation are valuable quantities on their own; they give 
important information on the solar wind. It is also shown that, at 
the same time, the declination and right ascension of a spacecraft 
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can potentially be determined much more accurately than by 
existing procedures. 

V02 DSN Progress Report for September-October 1971: 

Tropospheric and Ionospheric Range Corrections for an 
Arbitrary Inhomogeneous Atmosphere (First Order Theory) 

0. H. von Roos 

Technical Report 32-1526, Vol. VI, pp. 99-105, 

December 15, 1971 

In this article, a simple and concise expression is presented for the 
range correction for an atmosphere that possesses arbitrary radial, 
lateral, and azimuthal gradients of the index of refraction. The 
validity of this expression hinges only on the assumption that the 
index of refraction is close to unity, an assumption that is well- 
satisfied for the Earth's atmosphere. Furthermore, it is shown that 
the range corrections for a simple model of the Earth’s tropo- 
sphere, including typical lateral variations, are in close agreement 
with existing computer solutions. 


WAHLQUIST, H. D. 

W01 Hamiltonian Cosmology 

F. B. Estabrook and H. D. Wahlquist 

Phys. Lett., Vol. 35A, No. 6, pp. 453-454, July 12, 1971 

For abstract, see Estabrook, F. B. 


WELCH, L. R. 

W02 DSN Progress Report for July-August 1971: Combinational 
Complexity Measures as a Function of Fan-Out 

D. Johnson, J. E. Savage (Brown University), and 
L. R. Welch (University of Southern California) 

Technical Report 32-1526, Vol. V, pp. 79-81, 

October 15, 1971 

For abstract, see Johnson, D. 


WIEBE, E. 

W03 DSN Progress Report for July-August 1971: Superconducting 
Magnet for a Ku-Band Maser 

R. Berwin, E. Wiebe, and P. Dachel 
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Technical Report 32-1526, Vol. V, pp. 109-114, 
October 15, 1971 

For abstract, see Berwin, R. 


WILKINSON, D. T. 

W04 Preliminary Results of Laser Ranging to a Reflector on the 
Lunar Surface 

J. D. Mulholland, C. 0. Alley (University of Maryland), 

P. L. Bender (Joint Institute for Laboratory Astrophysics), 
D. G. Currie (University of Maryland), 

R. H. Dicke (Princeton University' J. E. Faller (Wesleyan 
University), W. M. Kaula (University of California, Los 
Angeles), G. J. F. MacDonald (University of California, 
Santa Barbara), H. H. Plotkin (Goddard Space Flight 
Center), and D. T. Wilkinson (Princeton University) 

Space Research XI, pp. 97-104, Akademie-Verlag, Berlin, 
1971 

For abstract, see Mulholland, J. D. 


WILLIAMS, J. G. 

W05 DSN Progress Report for July-August 1971: The Goldstone 
Interferometer for Earth Physics 

J. L. Fanselow, P. F. MacDoran, J. B. Thomas, 

J. G. Williams, C. Finnie, T. Sato, L. Skjerve (Philco-Ford 
Corporation), and D. Spitzmesser 

Technical Report 32-1526, Vol. V, pp. 45-57, 

October 15. 1971 

For abstract, see Fanselow, J. L. 


WILLIAMS, W. 

W06 Differential and Integral Cross Sections for the Electron-Impact 
Excitation of the a^ and b 1 ^* States of 0 2 

S. Trajmar, D. C. Cartwright (The Aerospace Corporation), 
and W. Williams 

Phys. Rev., Pt. A: Gen. Phys., Vol. 4, No. 4, pp. 1482 -1492, 
October 1971 

For abstract, see Trajmar, S. 
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WINKELSTEIN, R. 

W07 Minicomputer-Controlled Programmed <'* i.or 
R. Winkelstein 

JPL Quarterly Technical Review, Vol. 1, No. 3, pp. 79-87, 
October 1971 

The programmed oscillator is a telecommunications receiver or 
transmitter subsystem which compensates for the known doppler 
frequency effect produced by the relative motion between a 
spacecraft and a tracking station. Two such programmed oscilla- 
tors have been constructed, each using a low-cost minicomputer 
for the calculation and control functions, and each contained in a 
single rack of equipment. They are capable of operation in a 
phase-tracking mode as well as a frequency-tracking mode. When 
given an ephemeris suitable for the planet Venus, these units 
maintained phase coherence of better than 5 deg rms at 2388 
MHz. 


WONG, $. K. 

W08 Analysis of Mariner VII Pre-encounter Anomaly 

H. J. Gordon, S. K. Wong, and V. J. Ondrasik 

J. Spacecraft Rockets, Vol. 8, No. 9, pp. 931-937, 
September 1971 

For abstract, see Gordon, H. J. 


WOO, K. 

W09 DSN Progress Report for September-October 1971: Further 
Studies of Microwave Transmission Through Perforated Flat 
Plates 

T. Y. Otoshi and K. Woo 

Technical Report 32-1526, Vol. VI, pp. 125-129, 

December 15, 1971 

For abstract, see Otoshi, T. Y. 


WOO, R. 

W10 A Multiple-Beam Spherical Reflector Antenna 

R. Woo 
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JPL Quarterly Technical Review, Vol. 1, No. 3, pp. 88-96, 
October 1971 


A spherical reflector with multiple feeds is an attractive possibility 
for application in future communications satellite systems. Data 
are presented which show that spherical reflectors possessing rela- 
tively high gain (40 dB) and very small phase path error (<A/32) 
are feasible. A design of a spherical reflector utilizing corrugated 
horn feeds is considered. Radiation patterns are computed using 
the physical-optics technique. The designed antenna is approxi- 
mately 60A in diameter. Calculations performed for this antenna 
with three beams indicate that each beam has a gain of about 42 
dB, a beam width of 1.4 deg, and sidelobes that can be expected to 
be at least 28 dB down. These results indicate that the feature of 
low sidelobes makes the spherical reflector ^ promising candidate 
for a multiple-beam communications satellite antenna. 


YANG, J.-N. 

Y01 Maximum Dynamic Response and Proof Testing 

J.-N. Yang and E. Heer (NASA Headquarters) 

J. Eng. Mech. Div., Proc. ASCE, Vol. 97, No. EM4, 
pp. 1307-1313, August 1971 

Recent flight data show that the major spacecraft excitations 
during any one flight are not only highly transient, but are associ- 
ated with considerable statistical variation from one flight to an- 
other. However, the flight data are usually insufficient for the 
statistical characterization of random excitation inputs, so other 
means of characterization are being attempted. In the time-do- 
main analysis presented in this article, an approach is taken that 
yields an upper bound of the maximum dynamic response as well 
as a dynamic proof-testing excitation that, in turn, will produc? the 
expected upper bound of the maximum response. Although a 
structure with a single input is considered in this article, extension 
of the technique to a structure with multi-excitations appears to be 
possible. 


YASUI, R. K. 

Y02 Effects of Storage Temperatures on Silicon Solar Cell Contacts 
P. A. Berman and R. K. Yasui 
Technical Report 32-1541, October 15, 1971 
For abstract, see Berman, P. A. 
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Y03 


Supporting Data Package for TR 32-1541, Effects of Storage 
Temperatures on Silicon Solar Cell Contacts 

P. A. Berman and R. K. Yasui 

Technical Memorandum 33-497, October 15, 1971 

For abstract, see Berman, P. A. 



YEN, S. P. S. 

Y04 Electrical Conductivity of Elastomeric TCNQ Complexes Under 
Mechanical Stress 

A. M. Hermann (Tulane University), S. P. S. Yen, 

A. Rembaum, and R. F. Landel 

J. Polym. Sci., Pt. B: Polym. Lett., Vol. 9, No. 8, 
pp. 627-633, August 1971 

For abstract, see Hermann, A. M. 


YOUNG, A. T. 

Y05 Interpretation of Interplanetary Scintillations 

A. T. Young 

Astrophys. J., Vol. 168, No. 3, Pt. 1, pp. 543-562, 

September 15, 1971 

In this article, radio observations of interplanetary scintillations are 
interpreted by means of a theory previously used for the detailed 
interpretation of optical observations of atmospheric scintillations. 
The theory allows a number of restrictive assumptions to be 
removed, and thus gives a more realistic picture of the interplane- 
tary medium and a more accurate representation of the observa- 
tions. Observers are requested to report power spectra of the 
logarithm of the intensity, to avoid artifacts due to “spikes.” 


YOUNG, L. D. G. 

Y06 Comments on “The Venus Spectrum: New Evidence for Ice” 

R. A. J. Schorn and L. D. G. Young 

Icarus: Int. J. Sol. Sys., Vol. 15, No. 1, pp. 103-109, 
August 1971 

For abstract, see Schorn, R. A. J. 
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Y07 Calculation of the Partition Function for 14 N 2 16 0 
L. D. G. Young 

J. Quant Spectrosc. Radiat. Transfer, Vol. 11, No. 8, 
pp. 1265-1270, August 1971 

In calculating the transmission of the Earth’s atmosphere in the 
infrared, it is necessary to include the effects of such minor atmo- 
spheric constituents as H 2 0, C0 2 , CH 4 , N 2 0, CO, and O 3 . Before 
the line intensities for a given band can be accurately calculated, 
the partition function must be known. The most common isotope 
of nitrous oxide, 14 N 2 !e O, has been the subject of many high- 
resolution laboratory studies, and the best values of the molecular 
constants for this isotope have been tabulated elsewhere. In this 
article, internal and vibrational partition functions are tabulated 
for 14 N 2 ,6 0 for the temperature range 200-350°K at 10°K inter- 
vals. Rotational partition functions are also computed for the five 
lowest vibrational states. 


ZANDELL, C. 

Z01 DSN Progress Report for July-August 1971: CPS Sustaining 
Engineering 

C. Zandell 

Technical Report 32-1526, Vol. V, pp. 129-131, 

October 15, 1971 

The Sustaining Engineer Program for the Central Processing Sys- 
tem (CPS) of the Mark IIIA Space Flight Operations Facility is 
based on optimizing flight support capabilities. This is being 
achieved by testing hardware responses during simulated critical 
conditions, by supporting software development, and by monitor- 
ing system performance. This article defines the major hardware/ 
software problem areas and reports the results of studies made to 
resolve these problem areas. 
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